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used  to  monitor  tropopause  folding  events  under  a  variety  of  background  conditions  and 
a  preliminary  climatology  of  the  occurrence  of  folds  was  developed.  (5)The  spectrum 
of  vertical  motions  at  Flatland  was  analyzed  and  was  found  to  closely  match  predictions 
for  a  spectrum  of  Doppler-shifted  gravity  waves.  (6)Temporal  mean  vertical  velocities 
measured  by  the  radar  reflect  the  large-scale  flow  when  weather  systems  are  changing 
slowly  enough  that  the  radar  sample  is  representati ve  of  a  large  area..  ( 7 ) Fi nal 1 y ,  it 
was  found  that  the  sampling  strategy  for  vertical  motions  must  be  carefully  designed 
to  accommodate  the  frequency  spectrum. 
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ABSTRACT  OF  OBJECTIVES  AND  ACCOMPLISHMENTS: 

This  study  was  aimed  at  understanding  the  mesoscale  motions 
of  the  atmosphere,  their  interactions  with  other  scales  of 
flow,  and  applications  of  MST  radars  to  meteorological 
problems.  The  primary  observational  data  source  was 
vertical  and  horizontal  wind  measurements  from  MST  radars, 
also  called  clear-air  Doppler  profilers,  plus  some  analyses 
of  balloon  and  aircraft  data.  Investigations  included  (1) 
case  studies  of  the  interactions  between  motions  of  various 
scales,  including  study  of  the  influence  of  underlying 
terrain  on  the  interactions  of  meso-  and  larger-scale 
processes,  (2)  study  of  the  frequency  spectra  of  vertical 
and  horizontal  winds  as  a  function  of  background  weather 
conditions  such  as  stability  and  wind  speed,  and  (3) 
interpretation  of  the  observed  mesoscale  variations  in  the 
framework  of  current  theories  of  gravity  waves  and 
stratified  turbulence. 

It  was  concluded  from  analyses  of  aircraft  data  of  wind  and 
temperature  that  the  observed  scalar  spectra  are  reasonably 
consistent  with  a  simple  model  in  which  fluctuations  are 
produced  by  quasi  two-dimensional  motions  acting  in  a  plane 
inclined  slightly  to  isentropic  surfaces  and  to  the 
isopleths  of  constant  mixing  ratios  of  trace  species.  It 
was  also  found  that  the  variances  are  up  to  six  times  larger 
over  mountainous  terrain  than  over  oceans  or  plains.  The 
results  showed  only  partial  agreement  with  predictions  from 
the  theories  of  stratified  turbulence  and  of  a  spectrum  of 
gravity  waves  (see  paper  C.2),  indicating  that  further  study 
was  needed  to  understand  the  physics  of  the  motions 
involved.  MST  data  from  Colorado  and  Alaska  were  used  to 
examine  the  response  of  mesoscale  horizontal  motions  to 
changes  in  weather  conditions.  It  was  found  that  there  is 
an  enhancement  of  amplitude  at  high  frequencies,  with 
magnitude  of  about  a  factor  of  four  relative  to  values 
extrapolated  from  the  low  frequency  spectra,  regardless  of 
season,  altitude,  or  weather  conditions,  and  which  likely 
represents  high  frequency  gravity  waves  launched  by  flow 
over  rough  terrain.  Comparisons  of  the  results  from 
Colorado  with  theoretical  predictions,  including  the  effects 
of  Doppler-shifting,  again  showed  that  the  observed  spectra 


do  not  follow  either  the  wave  or  turbulence  models  closely 
in  all  respects. 

The  refractivity  turbulence  structure  constant,  C  2 ,  was 
measured  with  several  methods  and  instruments,  and 
comparisons  among  the  results  obtained  during  night  hours 
showed  good  agreement.  During  the  daylight  hours  the 
thermosonde  results  were  consistently  larger  than  those  from 
other  techniques.  An  important  result  of  this  study  was 
that  the  MST  radar  data  at  Flatland,  in  Illinois,  could  be 
modeled  using  the  same  model  parameters  developed  from  radar 
data  collected  in  Colorado,  except  that  the  distribution  of 
wind  shears  had  to  be  reduced  about  a  factor  of  two, 
suggesting  that  fine  structure  wind  shear  is  smaller  over 
the  plains  than  over  mountains,  and,  consequently,  that  eddy 
dissipation  rates  are  smaller  there. 

Several  results  were  obtained  using  the  Flatland  MST  radar: 
(1) there  is  enhanced  gravity  wave  activity  around  fronts, 
jet  streams,  and  convective  areas.  (2) The  radar  data  can  be 
used  to  study  tropopause  folding  events.  In  addition  to 
case  studies,  a  preliminary  climatology  of  folding  events 
was  developed.  (3) The  spectrum  of  vertical  motions  at 
Flatland  was  found  to  closely  match  that  predicted  for  a 
spectrum  of  gravity  waves  under  the  effects  of  Doppler 
shifting  by  the  horizontal  winds.  (4) Finally,  it  was  found 
that  the  radar  measurements  of  vertical  velocity  reflect  the 
large-scale  motions  when  storms  are  not  changing  rapidly  and 
conditions  are  statistically  stationary.  It  was  also  found 
that  care  must  be  used  to  ensure  that  adequate  temporal 
sampling  is  made  when  deriving  time  averages  of  vertical 
velocity;  this  is  due  to  the  relatively  large  energy  of 
vertical  motions  at  high  frequencies. 
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A-  RESEARCH  OBJECTIVES  (STATEMENT  OF  WORK) : 

1.  Study  the  joint  space  and  time  variability  of  small- 
scale  turbulence,  gravity  wave  activity,  and  background  flow 
using  data  from  the  Colorado  Profiler  Network  and  from  the 
Sunset  VHF  radar. 

2.  Draw  conclusions  regarding  the  influence  of  the 
underlying  terrain  on  the  interactions  of  the  processes 
occurring  at  different  scales. 

3.  Study  the  variations  of  vertical  velocity  and  small- 
scale  turbulence  during  summer  to  define  the  diurnal 
variations  and  determine  the  influence  of  background 
stability  on  gravity  wave  activity. 

4.  Compare  the  frequency  spectra  of  horizontal  and  vertical 
winds  during  different  seasons,  times  of  day,  locations, 
altitudes,  and  background  weather  conditions  to  draw 
conclusions  regarding  the  interpretation  of  the  mesoscale 
wind  spectra  in  terms  of  gravity  waves,  quasi  two- 
dimensional  turbulence  or  a  mix  of  processes. 

B.  STATUS  OF  THE  RESEARCH  EFFORT: 

Data  from  the  Platteville  MST  radar  in  the  Colorado 
Profiler  Network  were  used  to  study  changes  in  the  spectrum 
of  horizontal  winds  as  a  function  of  mean  wind  speed, 
altitude,  season,  and  stability.  It  was  found  that  the  high 
frequency  portion  of  the  spectrum  is  enhanced  relative  to  an 
extrapolation  of  the  low  frequency  regime.  This  high 
frequency  enhancement  is  likely  due  to  gravity  waves 
launched  by  flow  over  rough  terrain.  This  was  reported  in 
the  papers  listed  below  as  C.13  and  D.ll. 

Spectra  of  winds  and  temperatures  from  aircraft  were 
sorted  according  to  underlying  terrain,  and  it  was  found 
that  there  is  a  significant  enhancement  of  the  variance  over 
mountainous  regions.  This  was  reported  in  papers  C.2  and 
D.2. 


The  small-scale  eddies  that  lead  to  radar  reflections 
at  VHF  are  linked  to  the  eddy  dissipation  rate  and  to  the 
radar  refractivity  turbulence  structure  constant,  Cn2 . 
Measurements  of  C2  had  been  made  in  the  past  at  Sunset  and 
at  Denver,  and  haa  been  used  to  calibrate  a  theoretical 
model  of  Cn2  as  a  function  of  large-scale  winds  and 
stability.  *  This  model  was  applied  to  VHF  radar  data  from 
the  Flatland  radar  in  Illinois  and  it  was  found  that  the 
distribution  of  wind  shears  had  to  be  reduced  by  about  a 
factor  of  two  from  the  values  used  in  Colorado.  This 
indicates  the  fine  structure  wind  shear  over  the  plains  is 
smaller  than  over  the  mountains  for  the  same  large-scale 
wind  structure,  and  implies  that  eddy  dissipation  rates  are 
also  smaller.  This  was  reported  in  papers  C.7  and  D.9. 

Further,  an  effort  was  made  to  obtain  estimates  of  the 
vertical  flux  of  horizontal  momentum  using  data  from  the 
Sunset  radar  in  the  so-called  five-beam  configuration. 
Unfortunately,  the  data  were  spaced  in  time  such  that 
motions  with  time  scales  less  than  about  40  minutes  were  not 
resolved.  From  our  own  work  and  that  of  others  it  has 
become  clear  that  a  significant  amount  of  the  gravity  wave 
energy  is  at  periods  less  than  an  hour  or  so,  especially 
over  the  mountains  (see  paper  C.13),  which  makes  results 
based  on  coarse  time  resolution  of  limited  value.  Further, 
since  Sunset  is  near  the  continental  divide,  conditions 
change  rapidly  and  a  statistically  stationary  period  is 
difficult  to  define.  As  a  result  of  these  limitations,  it 
was  decided  to  delay  expanding  our  initial  analysis, 
reported  in  paper  D.l,  until  a  more  appropriate  data  set 
becomes  available. 

Simultaneous  measurements  of  horizontal  and  vertical 
velocity  at  Platteville  were  spectrally  analyzed  and  the 
results  compared  with  predictions  from  the  model  of  Doppler- 
shifted  gravity  waves  presented  by  Fritts  and  VanZandt  (J. 
Geophys.  Res.,  92:9723-9732,  1987).  Large  discrepancies 
were  noted,  but  as  Platteville  is  near  the  mountains  they 
could  be  due  to  another  process  such  as  mountain  wave 
effects.  This  was  reported  in  paper  D.3  and  expanded  upon 
in  paper  D.4 

Analyses  of  aircraft  data  were  contrasted  with 
predictions  from  the  theory  of  stratified  turbulence 
(reported  in  paper  C.l)  and  it  was  found  that  the  scalar 
spectra  are  reasonably  consistent  with  a  simple  model  in 
which  fluctuations  are  produced  by  quasi  two-dimensional 
motions  acting  in  a  plane  inclined  slightly  to  isentropic 
surfaces  and  to  the  isopleths  of  constant  mixing  ratios  of 
trace  species.  Similar  concepts  of  tilting  of  atmospheric 
layers  were  used  to  derive  a  simple  model  for  the  enhanced 
frequency  spectrum  of  vertical  velocities  and  are  reported 
in  papers  C.5,  C.ll,  and  D.7. 


Vertical  velocity  measurements  from  the  Flatland  VHF 
radar  in  Illinois  were  spectrally  analyzed  and  the  results 
compared  with  the  cheory  of  Doppler-shifted  gravity  waves  as 
presented  by  Fritts  and  VanZandt.  Excellent  agreement  was 
found,  on  a  statistical  basis,  suggesting  that  over  the 
plains  the  vertical  velocity  fluctuations  are  primarily  due 
to  propagating  gravity  waves.  This  was  reported  in  papers 
C . 14  and  D . 10 . 

The  Flatland  radar  data  were  also  used  to  study  the 
correlation  between  periods  of  increased  gravity  wave 
activity  and  the  background  weather  conditions.  As  reported 
in  papers  C.10  and  D.8,  there  is  evidence  that  every  weather 
system,  such  as  a  front  or  jet  streak  passage,  produces 
enhanced  gravity  wave  activity  (as  detected  by  the  radar) 
and  that  every  instance  of  enhanced  gravity  wave  activity 
detected  by  the  radar  could  be  traced  to  a  nearby  weather 
system. 

On  a  larger  scale,  the  Flatland  radar  was  used  to 
monitor  variations  in  tropopause  height  during  several  cases 
when  tropopause  folds  developed,  as  determined  from  analyses 
of  isentropic  potential  vorticity.  This  work,  reported  in 
papers  C.3  and  C.9,  demonstrates  the  usefulness  of  MST 
radars  for  monitoring  the  tropopause  and  its  folds  which  are 
believed  to  play  a  key  role  in  stratosphere-troposphere 
exchange  processes. 

Other  large-scale  processes  considered  during  this 
study  were  the  divergence  and  vertical  velocity.  These 
studies  were  reported  in  papers  C.4,  C.15,  D.5,  and  D.6. 

They  were  proof-of-concept  efforts  intended  to  show  the 
feasability  of  this  analysis  approach,  and  demonstrated  the 
great  versatility  of  MST  radar  data  in  application  to  many 
meteorological  analysis  and  forecasting  problems.  The  data 
are  relatively  new  to  meteorologists,  and  care  must  be  taken 
to  interpret  them  properly  (as  discussed  in  papers  C.6  and 
C.12).  Other  applications  of  the  data  to  broad  problems 
were  discussed  in  a  major  review  session,  reported  in  paper 
C.8. 
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As  reflected  by  the  coauthors  of  the  papers  listed 
above,  much  of  this  work  was  done  in  collaboration  with 
other  scientists.  Joint  studies  with  NOAA's  Aeronomy 
Laboratory  and  several  university  groups  were  made,  with 
contact  provided  during  personal  visits,  seminars  or 
workc-  ;  ,  and  at  conferences  attended  as  reflected  in  the 

confe  .  .  ,.je  papers  listed. 

In  u  member,  1987,  Dr.  G.D.  Nastrom  participated  in  the 
8th  iJi.iual  EOSAEL/TWI  Conference  (Las  Cruces,  NM)  .  The 
joint  studies  of  refractivity  turbulence  involving  groups 


from  the  Air  Force  Geophysics  Laboratory,  the  Army 
Atmospheric  Sciences  Laboratory,  and  the  Flatland  MST  radar 
were  an  outgrowth  of  discussions  at  this  conference. 

Finally,  peer  reviews  rendered  on  journal  papers  and  on 
proposals  for  new  work  represent  an  interaction  with  the 
scientific  community  in  a  broad  sense.  During  the  period  of 
this  contract  Dr.  Nastrom  provided  reviews  on  22  papers  for 
the  Journal  of  Atmospheric  Sciences,  Journal  of  Atmospheric 
and  Terrestrial  Physics,  Geophysical  Research  Letters, 
Journal  of  Geophysical  Research,  Radio  Science,  and  Journal 
of  Climate  and  Applied  Meteorology.  Also,  reviews  were 
given  for  8  proposals  to  the  Air  Force  Office  of  Scientific 
Research  and  the  National  Science  Foundation. 
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Spectrum  of  Atmospheric  Vertical  Displacements  and 
Spectrum  of  Conservative  Scalar  Passive  Additives 
Due  to  Quasi-Horizontal  Atmospheric  Motions 


K.  S.  Gage 


Laborator  \ .  \j;n>nd/  ami  Atmospheru  Administration  Environmental  Research  Laboratories. 

Boulder.  Colorado 


G  D.  Nastrom 


Meteor{>!og\  Research.  C ontrol  Data  Corporation.  Minneapolis .  Minnesota 

This  paper  is  concerned  with  the  relationship  between  the  spectrum  of  vertical  displacements  associ¬ 
ated  with  the  mesoscale  spectrum  of  quasi-horizontal  motions  in  the  atmosphere  and  the  mesoscale 
spectra  of  scalar  passive  additives  This  work  is  based  on  observations  of  mesoscale  atmospheric  wave 
number  specira  of  wind,  temperature,  and  certain  trace  species  sampled  in  the  altitude  range  9-14  km  bv 
commercial  aircraft  during  the  Global  Atmospheric  Sampling  Program  (GASP)  The  principal  result  is 
that  the  scalar  spectra  are  reasonable  consistent  with  a  simple  model  in  which  scalar  fluctuations  are 
produced  bv  quasi  two-dimensional  motions  acling  in  a  plane  inclined  slightly  to  isentropic  surfaces  and 
to  the  isopleths  of  constant  mixing  ratio  of  trace  species. 


1  Introduction 

Following  the  work  of  Reed  and  German  [1965]  and  many 
others,  it  has  been  widely  recognized  that  much  of  the  large- 
scale  vertical  transport  in  the  atmosphere  is  due  to  quasi- 
horizontal  motions  that  transport  mass  to  great  distances 
along  sloping  isentropic  surfaces.  Large-scale  motions  in  the 
atmosphere  are  quasi-horizontal  and  are  much  more  effective 
in  diffusing  mass  in  the  horizontal  than  in  the  vertical  direc¬ 
tion  Yet  since  this  diffusion  takes  place  along  sloping  surfaces, 
there  is  a  nonnegligible  vertical  component  of  mass  transport. 
The  diffusivity  along  the  mixing  surfaces  so  greatly  exceeds  the 
diffusi \ i ty  normal  to  the  surfaces  of  constant  concentration 
that  even  a  very  small  slope  of  the  mixing  surfaces  is  enough 
to  produce  a  major  contribution  to  vertical  transport. 

In  this  paper  we  are  concerned  primarily  with  the  mesoscale 
contribution  to  the  process  of  vertical  transport  due  to  the 
spectrum  of  nearly  horizontal  motions  that  has  been  studied 
using  wind  and  temperature  data  collected  during  routine 
commercial  airline  flights  during  the  Global  Atmospheric 
Sampling  Program  IGASP)  [Nastrom  et  al..  1984;  Nastrom 
and  Gage.  1985]  The  interpretation  of  these  mesoscale  spectra 
offered  by  Gage  and  Nastrom  [1986]  points  toward  the  exis¬ 
tence  of  a  mesoscale  spectrum  of  quasi  two-dimensional  tur¬ 
bulence  which  should  contribute  very  effectively  to  horizontal 
diffusivity.  Analysis  of  the  temperature  spectra  suggests  that 
they  are  also  the  result  of  quasi  two-dimensional  turbulence 
acting  on  the  background  temperature  gradient  The  quasi 
two-dimensional  turbulence  acting  on  a  background  gradient 
of  the  mixing  ratio  of  a  passive  scalar  quantity  should,  in  a 
similar  way.  lead  to  a  mesoscale  spectrum  of  scaiar  fluctu¬ 
ations 

In  this  paper  we  develop  a  physical  model  linking  the  ob- 
*erved  GASP  spectra  of  wind,  temperature,  and  trace  species 
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In  a  companion  paper  by  Nastrom  et  al.  [this  issue]  we  con¬ 
sider  the  consistency  of  the  observed  spectra  with  the  physical 
model  proposed  here. 

1  Observed  Spectra  of  Atmospheric  Motions 
and  Temperature  From  the  GASP  Data  Set 

A  detailed  climatology  of  mesoscale  horizontal  wave 
number  spectra  of  wind  at  tropopause  altitude  has  been  pre¬ 
sented  by  Nastrom  et  al.  [1984]  and  Nastrom  and  Gage 
[1985],  These  spectra  were  deduced  from  an  analysis  of  data 
from  over  6900  routine  commercial  airline  flights  which  com¬ 
prised  GASP  The  wind  spectra  are  shown  in  summary  form 
in  Figure  1.  The  spectra  show  two  regimes  of  distinctly  differ¬ 
ent  slope.  The  mesoscale  regime  extends  from  about  500  km 
or  so  on  the  large-scale  end  down  to  a  few  kilometers.  Here 
the  slope  is  —5  3.  as  expected  for  quasi  two-dimensional  stra¬ 
tified  turbulence  [Gage.  1979:  Lilly.  1983]  At  large  scales  the 
spectra  steepen  toward  the  —3  slope  expected  for  the 
enstrophy-cascading  range  of  quasi  two-dimensional  turbu¬ 
lence.  which  is  usually  referred  to  as  geostrophic  turbulence 
[C/iurney.  1971],  A  schematic  representation  of  these  two  iner¬ 
tial  ranges  shown  in  Figure  2  [ Larsen  et  al..  1982]  fits  the 
GASP  spectra  quite  well. 

The  wave  number  spectrum  of  potential  temperature  ob¬ 
tained  from  analysis  of  the  GASP  data  is  shown  in  summary 
form  in  Figure  1  and  for  the  longest  flights  in  Figure  3  It  is 
obvious  at  first  glance  that  the  temperature  spectrum  is  basi¬ 
cally  the  same  shape  as  the  velocity  spectra.  Note  that  this 
result  differs  from  model  calculations  of  quasi  two- 
dimensional  turbulence,  which  predicts  the  slope  of  a  passive 
scalar  quantity  to  be  —  1  in  the  enstrophy-cascading  range 
and  +1  3  in  the  energy-cascading  range  [Lesieur  and  Herring. 
1985]  Reasons  for  this  will  be  discussed  in  section  4. 

Evidently,  the  magnitude  and  shape  of  the  potential  temper¬ 
ature  spectrum  are  determined  by  the  same  dynamics  that 
govern  the  velocity  spectra  Indeed,  as  discussed  by  Gage  and 
Nastrom  [1986],  the  potential  temperature  spectrum  appears 
to  be  in  equilibrium  with  the  velocity  spectra,  as  evidenced  by 
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Fig  I  Summary  comparison  of  GASP  wave  number  spectra  of  potential  temperature  and  zonal  and  meridional 
velocity  composited  from  three  groups  of  flight  segments  of  different  lengths.  The  three  types  of  symbols  show  results  from 
each  group  Straight  lines  indicate  slopes  of  -  3  and  -  5  3.  The  mendtonal  wand  spectra  arc  shifted  one  decade  to  the 
right,  and  the  potential  lemperature  spectra  are  shifted  two  decades  to  the  right  [after  Sastrom  and  Gage.  1985]. 


the  relationship  of  the  potential  energy  spectrum  to  the  kinetic 
energy  spectrum  As  anticipated  by  Charney  [1971]  for  geo- 
strophic  turbulence,  the  spectral  energy  is  partitioned  approxi¬ 
mately  equally  among  the  two  components  of  kinetic  energy 
and  potential  energy.  This  result  is  reproduced  in  Figure  4. 
The  potential  temperature  spectrum  is  related  to  the  potential 
energy  spectrum  by 


where  N  is  the  Brunt-Vaisala  frequency  and  g  is  gravitational 
acceleration  Since  N1  is  larger  in  the  stratosphere  than  in  the 
troposphere,  the  amplitude  of  the  potential  temperature  spec¬ 
trum  is  larger  in  the  stratosphere  than  the  troposphere.  These 
relationships  appear  to  hold  throughout  the  atmospheric 
mesoscale  as  well  as  in  the  geostrophic  regime. 

Note  that  the  equipartitioning  referred  to  above  is  evidently 
intrinsic  to  the  quasi-horizontal  motions  themselves.  This  does 
not  imply  an  equipartitioning  between  the  quasi-horizontal 
motions  and  wave  motions,  as  discussed  by  Gage  and  Nastrom 
[1986],  In  this  context  the  low  Froude  number  regime  of  stra¬ 
tified  turbulence  [Riley  et  a/.,  1981;  Lilly.  1983]  implies  a 
regime  of  decoupled  modes  of  interna!  waves  and  quasi¬ 
horizontal  motions.  Indeed,  if  the  stratification  is  sufficiently 
stable,  the  kinetic  energy  associated  with  the  quasi-horizontal 


mode  of  motion  can  be  expected  to  dominate  the  kinetic  and 
potential  energy  associated  with  the  internal  wave  mode. 

In  section  3  we  consider  the  consequences  of  a  spectrum  of 
quasi-horizontal  eddies  flowing  along  sloping  "mixing  sur- 


Fig.  2  Schematic  representation  of  spectral  ranges  in  quasi  two- 
dimensional  turbulence  [after  Larsen  el  al .  1982] 
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Fig  3  Poiemi.il  temperature  spectra  from  flights  at  least  4800 
km  long  Stratospheric  and  tropospheric  spectra  are  plotted  separate¬ 
ly.  The  number  of  cases  averaged  in  each  group  is  given  in  parenthe¬ 
ses  [after  Gam  and  Nastrom.  1986] 


faces."  as  anticipated  in  the  analysis  of  Reed  and  German 
[1965].  Specifically,  we  shall  consider  what  happens  when 
these  "mixing  surfaces"  are  inclined  relative  to  isentropic  sur¬ 
faces  or  isopleths  of  constant  mixing  ratio  of  trace  species  such 
as  ozone 


3.  Spectrum  of  Vertical  Displacement  Deduced 
From  the  GASP  Temperature  Spec  tri  m 

Fluctuations  in  any  conservative  passive  scalar  quantity  will 
be  produced  whert*vcr  there  is  a  component  of  motion  parallel 
to  the  gradient  of  the  scalar  quantity.  In  the  atmosphere,  and 
especially  in  the  stratosphere,  gradients  of  conservative  passive 
scalar  quantities  usually  are  very  much  larger  in  the  vertical 
than  in  the  horizontal  directions.  This  is  just  another  way  of 
saying  the  atmosphere  is  horizontally  stratified,  and  to  first 
order,  gradients  of  potential  temperature,  ozone,  etc  .  are  in 
the  vertical  direction.  The  exception  being  in  the  vicinity  of 
fronts,  strong  convection,  etc. 

For  the  purposes  of  this  paper  we  shall  define  "vertical 
displacements"  relative  to  isentropic  surfaces  Thus  quasi- 
hortzontal  motions  which  are  coincident  with  isentropic  sur¬ 
faces  would  not  produce  "vertical  displacements”  even  if  the 
isentropic  surfaces  were  inclined  to  the  horizontal  Our  con¬ 
cern  is  rather  with  the  small  displacements  taking  place  or¬ 
thogonal  to  the  isentropic  surfaces  or  isopleths  of  constant 
mixing  ratio  From  this  perspective  the  vertical  displacement 
spectrum  G>..  is  the  projection  of  the  quasi-horizontal  displace¬ 
ment  spectrum  in  the  mixing  surface  <t>, ,  in  the  direction 
normal  to  the  isentropic  surface  Thus  the  vertical  displace¬ 
ment  spectrum  is  given  by 

d>;;  =  (sin  (3i 

where  5  is  the  angle  that  the  isentropic  surface  makes  with  the 
mixing  surface. 

With  vertical  displacement  defined  in  this  way  the  vertical 
displacement  spectrum  can  be  deduced  from  the  potential 
temperature  spectrum,  provided  the  background  potential 
temperature  gradient  is  known.  If  the  background  potential 
temperature  gradient  is  constant,  the  relationship  between  the 
potential  temperature  spectrum  and  the  vertical  displacement 
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Fig  4  Kinetic  and  potential  energy  spectra  from  flights  at  least  4800  km  long  (a)  tropospheric  energy  and  (b)  strato¬ 
spheric  energy  spectra  [after  G age  and  Nastrom,  1986] 
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where  ii  represent  ihe  angle  of  ihe  isopleth*  of  convlam 
mixing  ratio  relame  to  the  mixing  surface*  Then  the  fluctu¬ 
ations  of  conservative  passive  scalar  quantity  will  satisfy 

The  ratio  of  slopes  sin  O'  sin  A  can  be  evaluated  as  follow* 
The  variance  of  each  spectrum  can  be  obtained  by  integrating 
aero**  the  wave  number  spectrum 

A:  =  j  d>_  i Ik 
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f  rom  ( 3 1  we  find 


Fig  5  Vertical  displacement  spectra  in  the  troposphere  and  _  rr„„ 

stratosphere,  deduced  from  the  potential  temperalure  spectrum  and  ~  • 

assumed  background  gradients  of  potential  temperature 

from  l " )  we  find 


spectrum,  as  defined  above,  is  given  by 


<D;; 


(cfl  r:\-  ,\- 


and  from  161  we  find 
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where  we  have  made  use  of  ( 1 1  The  vertical  displacement 
spectrum  defined  in  this  way  is  easily  deduced  from  either  the 
potential  temperature  spectrum  of  Figure  3  or  the  potential 
energy  spectrum  of  Figure  4  The  result  is  shown  in  Figure  5 
Note  that  the  magnitude  of  the  displacement  spectrum  is  re¬ 
duced  in  the  stratosphere  relative  to  the  troposphere  Evident¬ 
ly,  stable  stratification  causes  stratospheric  motions  to  deviate 
less  from  isentroptc  surfaces  than  tropospheric  motions. 

Consider  now  any  conservative  passive  scalar  quantity  / 
that  has  isopleths  of  constant  mixing  ratio  coincident  with 
isentropic  surfaces  It  would  then  follow  that 

(^)*  *4) 

provided  the  gradient  of  ■/  is  constant  within  the  region  where 
/  is  measured  Under  these  circumstances  we  would  anticipate 


and  this  relationship  can  be  used  to  check  for  consistency 
between  passive  scalar  spectra  and  potential  temperature  spec¬ 
tra.  The  results  of  this  calculation  are  given  by  Sasirom  e t  al 
[this  issue")  and  show  reasonable  agreement  among  spectral 
amplitudes 

More  generally,  we  expect  that  the  slopes  of  isopleths  of 
constant  mixing  ratio  will  not  coincide  with  isentropic  sur¬ 
faces  In  this  event  the  vertical  displacements  defined  relative 
to  the  isopleths  of  constant  mixing  ratio  will  difer  from  the 
displacements  relative  to  the  isentropic  surfaces  so  that 


(A  ):  / sin  h'  \ : 
4. A* I  \sin  6  J 

From  these  relationships  we  deduce 


The  quantities  A  and  A'  are  analogous  to  Ehhalt's  [ Ehluilr  t‘i 
al..  1983]  equivalent  vertical  displacement. 

4  Discussion 

In  this  paper  we  have  developed  a  simple  model  relating  the 
spectrum  of  conservative  passive  scalars  to  the  spectrum  of 
vertical  displacements  of  quasi-horizontal  motions.  The  model 
provides  a  straightforward  interpretation  of  the  spectra  of 
temperature  and  trace  species  observed  during  GASP 

For  the  development  presented  here  we  have  focused  on  the 
mesoscale  atmospheric  motions  and  assumed  the  background 
mean  concentrations  and  associated  isopleths  of  constant 
mixing  ratio  to  be  given  In  reality,  of  course,  isentropic  sur¬ 
faces  and  mixing  ratio  isopleths  must  evolve  out  of  the  large- 
scale  dynamics,  as  described  by  Mahlman  el  al.  [1984]  and  by 
Mahlman  [1985]  According  to  this  picture,  the  slopes  of  con¬ 
stant  mixing  ratios  of  trace  species  are  determined  by  a  bal¬ 
ance  between  the  meridional  circulation,  which  tends  to  in¬ 
crease  their  slope,  and  the  quasi-horizontal  eddy  motions, 
which  tend  to  decrease  their  slope  Radiation  has  the  effect  of 
changing  the  potential  temperature  of  parcels  which  are  other¬ 
wise  assumed  to  remain  on  isentropic  surfaces  Chemical  reac¬ 
tions  can  alter  this  picture  for  reactive  species  and  these 
species  can  have  slopes  differing  from  nonreactive  species 
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The  dynamical  model  developed  here  attributes  the  fluctu¬ 
ations  of  scalar  quantities  solely  to  the  vertical  displacements 
experienced  by  parcels  of  fluid  in  quasi-horizontal  motion 
The  motion  field  in  the  presence  of  a  mean  gradient  of  the 
conservative  scalar  quantity  is  from  this  perspective  respon¬ 
sible  for  the  observed  scalar  fluctuations.  This  picture  of  pas¬ 
sive  scalar  transport  implies  a  spectrum  of  scalar  fluctuations 
closely  tied  to  the  spectrum  of  fluid  motions  There  is  no 
reason  why  a  spectrum  of  scalar  fluctuations  determined  in 
this  way  should  bear  any  particular  relationship  to  the  scalar 
spectrum  predicted  from  turbulence  model  simulations  [Les- 
ieur  and  Herring,  1985]  The  latter  spectra  are  produced  when 
scalar  variance  in  an  otherwise  homogeneous  fluid  is  added  to 
the  turbulence  model  at  specified  wave  numbers.  The  scalar 
quantity  then  participates  in  the  turbulence  cascades  of  the 
model,  and  a  scalar  spectrum  is  produced  For  example,  when 
scalar  variance  is  added  at  small  wave  numbers  to  a  turbu¬ 
lence  model  that  simulates  the  enstrophy  cascade  of  geo- 
strophic  turbulence,  a  scalar  variance  spectrum  with  spectral 
slope  —  1  is  produced  These  model  simulations  are  unlikely  to 
simulate  the  spectra  observed  in  GASP,  since  they  do  not 
account  for  the  gradients  in  mean  scalar  concentration. 

Throughout  this  work  we  have  assumed  that  the  motions 
represented  by  the  GASP  velocity  spectra  are  quasi-horizontal 
and  turbulencelike.  Dewan  [1979]  and  VanZandt  [1982]  sug¬ 
gested  that  such  spectra  could  be  due  to  internal  waves  famil¬ 
iar  to  oceanographers  as  the  Garrett-Munk  spectrum  [Garrett 
and  Munk.  1972,  1975;  Others,  1983],  and  the  two  hypotheses 
are  contrasted  next. 

Our  main  reasons  for  asserting  that  the  GASP  spectra  are 
not  internal  wave  spectra  are  given  in  detail  by  Gage  and 
Nastrom  [1986]  The  principal  problem  with  the  Dewan- 
VanZandt  hypothesis  can  be  stated  very  succinctly:  There  is 
too  much  energy  in  the  spectrum  of  horizontal  atmospheric 
motions  compared  to  the  spectrum  of  atmospheric  vertical 
motions  to  be  consistent  with  the  idea  that  both  spectra  are 
manifestations  of  a  common  spectrum  of  internal  waves.  The 
background  vertical  velocity  spectrum  in  the  atmosphere  does 
appear  to  have  a  striking  resemblance  to  the  Garrett-Munk 
spectrum  in  the  ocean,  but  the  horizontal  wind  spectral  ampli¬ 
tude  consistent  with  this  is  very  much  smaller  than  the  ampli¬ 
tudes  observed.  A  preliminary  climatology  of  vertical  velocity 
spectra  taken  at  various  geographical  locations  by  vertically 
looking  clear-air  Doppler  radar  under  light-wind  conditions 
shows  a  flat  spectrum  of  fairly  universal  amplitude  with  a 
sharp  cutoff  at  frequencies  higher  than  the  Brunt-Vaisala  fre¬ 
quency  [ Ecklund  et  al,  1986],  They  show  that  the  variance  of 
this  spectrum  is 

Jmf*  A 

<*»„»(/)  =  —  *  1 -0-2.0  x  10'2  m2  s'2  (13) 
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where  the  frequency  spectrum  of  vertical  motions  is  approxi¬ 
mated  by 

/;</£/* 


Otherwise, 

<t>  =0 
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where  /.( =  I/t<)  is  the  inertial  frequency  and  /*<  s  !/r,)  is  the 
buoyancy  frequency  N/2n.  The  variance  of  horizontal  velocity 
for  a  spectrum  of  internal  waves  can  be  estimated  as 
r/w 

U7  *  J  <*>„(/)  =  Amr,  (14) 


where  the  frequency  spectrum  of  zonal  horizontal  motions  is 
approximated  by 

/.</</, 


Otherwise. 


<t>  =  0 


The  constants  Aww  and  can  be  related  by  the  polarization 
relation  in  the  form 


£,  w2 


where  to  =  2 nf.  £„  =  I®...  EH  ~  | and  we  have 

assumed  =  <J>tl  The  value  of  to  for  which  <!>„„  =  <J>H„  gives 
the  desired  relationship  between  Auu  and  AWH„  namely. 


(16) 


The  ratio  of  variances  of  horizontal  and  vertical  motions  is 
then  found  to  be 


I  li 
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If  we  take  rB  =  10  min.  ?  =  2.0  *  10“2  m2  s'2,  and  r,  =  18 
hours,  corresponding  to  approximately  42  latitude,  we  find 
that  ti1  *  0.72  m2  s'2.  This  variance  associated  with  the  inter¬ 
nal  wave  spectrum  is  about  an  order  of  magnitude  less  than 
the  variance  contained  in  the  mesoscale  spectrum  of  quasi¬ 
horizontal  eddies  considered  in  this  paper  and  evident  in  the 
GASP  spectra  [cf..  Gage  and  Nastrom,  1985.  1986] 


5.  Conclusions 

In  this  paper  we  have  developed  a  simple  physical  model  to 
relate  the  spectra  of  conservative  passive  scalar  tracers  to  the 
spectra  of  vertical  displacements  associated  with  quasi¬ 
horizontal  atmospheric  motions.  The  vertical  displacement 
spectra  were  deduced  from  the  GASP  spectra  of  potential 
temperature.  The  consistency  between  the  spectral  shape  and 
amplitude  of  observed  GASP  spectra  with  this  simple  model 
supports  the  hypothesis  that  mesoscale  quasi-horizontal  mo¬ 
tions  play  an  important  role  in  vertical  as  well  as  horizontal 
transport  for  the  atmospheric  mesoscale  much  as  they  do  in 
large-scale  planetary  wave  motions.  The  existence  of  an  entire 
spectrum  of  quasi  two-dimensional  motions  has  important  im¬ 
plications  for  the  horizontal  and  vertical  transport  of  trace 
species  in  the  atmosphere.  At  the  very  least,  local  mixing  by 
this  process  should  complement  the  large-scale  transport  as¬ 
sociated  with  planetary  waves  and  geostrophic  turbulence. 

Acknowledgment  One  of  us  (G  D.N.)  was  partially  supported  by 
the  Air  Force  Office  of  Scientific  Research 

References 

Charnev.  J  G  ,  Geostrophic  turbulence.  J.  Atmos  St  i .  2k.  1087-1095. 
197).' 

Dewan.  E  M  .  Stratospheric  spectra  resembling  turbulence.  Science, 
204.  832-835.  1979. 

Ecklund.  W  L.,  K  S  Gage.  G.  D  Nastrom.  and  B  B  Balsley  A 
preliminary  climatology  of  the  atmospheric  vertical  velocity  spec¬ 
trum  observed  by  clear-air  radar,  J  dim.  Appi  Meteorol..  22.  885 
892.  1986 

Ehhalt.  D.  H  ,  E  P  Roth,  and  U.  Schmidt.  On  the  temporal  variance 
of  stratospheric  trace  gas  concentrations.  J  Atmos  Chem  .  I.  27-51. 
1983 

Gage.  K  S  .  Evidence  for  a  L~!3  law  inertial  range  in  mesoscale 
two-dimensional  turbulence.  J  Atmos  Set .  }6.  1950-1954.  1979 


n.:ib 


Gagl  and  Nastrom  Spectrum  From  Quasi-Horizontal  Atmospheric  Motions 


Gage.  K  S  .  and  G  D  Nastrom.  On  the  spectrum  of  atmospheric 
velocity  fluctuations  seer.  by  ST  MST  radar  and  their  interpreta¬ 
tion.  Radio  Sc r.  20.  1339  1347.  1985 
Gage.  K  S  .  and  G  D  Nastrom,  Theoretical  interpretation  of  atmo¬ 
spheric  wavenumber  spectra  of  wind  and  temperature  observed  b> 
commercial  aircraft  during  GASP.  J.  Atmos.  See.  A3.  729-740.  1986 
Garrett.  C  .  and  Munh.  Space-time  scales  of  internal  waves.  Geo- 
ph\s  Astrophys.  Fluid  D\n..  2.  225-264.  1972 
Garrett.  C.  and  W  Munk.  Space-time  scales  of  internal  waves  A 
progress  report.  J  Geophys  Res..  #0.291  297.  1975. 

Larsen.  M  F.  M  C  Kelley  and  K  S  Gage.  Turbulence  spectra  in 
the  upper  troposphere  and  lower  stratosphere  between  2  hours  and 
40  days.  J  Atmos  Set  .  39.  1035-1041.  1982 
Lesieur.  M  .  and  J  Herring,  Diffusion  of  a  passive  scalar  in  two- 
dimensional  turbulence.  J  Fluid  Mech..  161.  77-95.  1985 
Lilly  D  K  .  Stratified  turbulence  and  the  mesoscale  variability  of  the 
atmosphere.  J  Armos.  See.  40.  '49-761,  1983 
Mahlman.  J  D  .  Mechanistic  interpretation  of  stratospheric  tracer 
transport.  Adi  Geophys..  2# A.  301-323,  1985 
Mahlman.  J  D.  G  Andrews.  H  V  Dutsch.  D  L  Hartmann.  T 
Matsuno,  and  R  J  Murgatrovd.  Transport  of  trace  constituents  in 
the  stratosphere,  in  Dynamos  of  the  Middle  Atmosphere,  edited  by  J. 
R  Holton  and  T  Matsuno.  pp  387-416.  Terra  Reidel.  Tokyo. 
1984 

Nastrom.  G  D  .  and  K  S.  Gage.  A  climatology  of  atmospheric  wave- 
number  spectra  observed  by  commercial  aircraft.  J  Atmos.  Set .  42. 
950-960.  1985 


Nastrom.  G  D.  K  S  Gage,  and  W  H  Jasperson.  Kinelu  energy 
spectrum  of  large-  and  mesoscale  atmospheric  processes,  .\uiur, 
310.  36  38.  1984 

Nastrom,  G  D  .  W  H  Jasperson,  and  K  S  Gage.  Horizontal  spectra 
of  atmospheric  tracers  measured  during  the  Global  Atmospheric 
Sampling  Program.  J  Geophys.  Res.,  this  issue 

Olbers.  D  J  .  Models  of  the  oceanic  internal  wave  held.  Ret  Geoph\ s 
2/.  1567  1606,  1983 

Reed.  R  J  .  and  K  E  German,  A  contribution  10  the  problem  of 
large-scale  mixing.  Mon.  Weather  Ret  ..  93.  313  321.  1965 

Riley.  J  J  .  R  W  Metcalfe,  and  M.  A  Weissman.  Direct  numerical 
simulations  of  homogeneous  turbulence  in  densny-siralitied  fluids, 
in  Proceedings  of  AIP  Conference  on  \onhnear  Properties  ot  Inter¬ 
nal  Wares,  edited  by  Bruce  J  West,  pp  79-  1 12.  American  Institute 
-of  Physics.  New  York.  1981 

VanZandt.  T  E.,  A  universal  spectrum  of  buoyancy  waves  in  ihe 
atmosphere.  Geophys  Res  Lett .  9,  575-578,  1982 


K  S  Gage.  Aeronomy  Laboratory.  National  Oceanic  and  Atmo¬ 
spheric  Administration.  Environmental  Research  Laboratories.  325 
Broadway.  Boulder.  CO  80303 

G  D  Nastrom.  Meteorology  Research.  Control  Data  Corporation. 
Box  1249-B,  Minneapolis.  MN  55420. 

(Received  November  28.  1985. 
revised  July  28.  1986. 
accepted  July  29.  1986  1 


I 


Paper  C.2 

Reprinted  from  Journal  of  the  Atmospheric  Sciences.  Vol  44.  No  20.  15  October  1987 

American  Meteormopcal  Sonny 


An  Investigation  of  Terrain  Effects  on  the  Mesoscale  Spectrum  of  Atmospheric  Motions 

G.  D.  Nastrom 

Meteorology  Research  Control  Data.  Minneapolis.  MS  55440 

D.  C.  Fritts 

Geophysical  Institute.  University  of  Alaska.  Fairbanks.  AK  V9~?5 


K..  S.  Gage 

Agronomy  Laboratory.  SO.A.A/ERL.  Boulder.  CO  SOM)} 

(Manuscript  received  3  February  1987.  in  final  form  13  May  1987) 

ABSTRACT 

Wind  and  temperature  data  collected  on  commercial  aircraft  dunng  the  Global  Atmospheric  Sampling  Program 
(GASP)  are  used  to  investigate  the  effects  of  underlying  terrain  on  mesoscale  variability,  and  the  observational 
results  are  interpreted  within  the  theories  of  gravity  wave  motions  and  quasi-two-dimensional  turbulence.  The 
data  show  the  variances  are  up  to  six  times  larger  over  mountainous  terrain  than  over  oceans  or  plains,  with 
the  most  striking  differences  at  horizontal  scales  from  4  to  80  km  Results  were  subdiv  ided  between  the  stratosphere 
and  troposphere,  and  between  high-  and  low -background  wind  speed  cases,  and  show  basically  the  same  response 
to  topography  in  all  cases  The  linear  theory  of  grav  ity  waves  is  found  to  predict  correctly  the  scaling  of  wave 
amplitude  with  background  stability  in  the  case  of  low-background  wind  speeds,  while  the  two-dimensional 
turbulence  theory  correctly  predicts  the  shape  of  the  variance  spectrum  and  the  observ  ed  amplitude  scales  w  ith 
f  as  required  by  the  theory  .  In  other  cases  the  theoretical  predictions  are  less  satisfactory  Some  possible  causes 
of  the  discrepancies  and  likely  methods  to  resolve  them  are  discussed 


I.  Introduction 

It  has  become  increasingly  apparent  in  recent  years 
that  mesoscale  variations  of  the  atmosphere  have  an 
important  impact  on  a  wide  variety  of  meteorological 
problems.  Part  of  this  awareness  comes  from  the  use 
of  improved  measurement  systems  that  permit  us  more 
adequately  to  sample  mesoscale  motions,  and  part  of 
it  comes  from  theoretical  advances  that  highlight  the 
relevance  of  the  mesoscale  in  meteorological  processes 
(e  g..  Atkinson.  1981;  Lilly  and  Gal-Chen,  1983).  De¬ 
spite  this  progress,  many  observational  and  theoretical 
problems  remain.  For  example,  the  relationship  of  me¬ 
soscale  variability  near  the  tropopause  with  the  un¬ 
derlying  topography  is  of  interest  for  several  reasons. 
One  prominent  reason  is  that  mountain  drag  plays  a 
significant  role  in  the  general  circulation  of  the  at¬ 
mosphere  (White,  1949).  Palmer  et  al.  (1986)  have  re¬ 
cently  pointed  out  the  possible  need  to  model  more 
closely  mountain  wave  drag  in  order  to  achieve  more 
realistic  simulations  of  the  tropospheric  and  lower 
stratospheric  circulation.  Dedicated  field  programs  such 
as  ALPEX  can  help  us  to  understand  mountain  drag 
(Davies  and  Phillips,  1985),  and  observations  from  less 
closely  focused  data  collection  efforts  also  can  provide 
information  on  the  coupling  of  mesoscale  variance  and 
underlying  topography.  In  a  study  of  the  climatology 


of  the  variance  power  spectrum  of  winds  and  temper¬ 
atures  at  airline  cruise  altitudes  based  on  data  from  the 
Global  Atmospheric  Sampling  Program  (GASP).  Nas¬ 
trom  and  Gage  (1985.  hereafter  called  NG)  found, 
among  other  things,  that  the  variance  associated  with 
motions  for  scales  less  than  about  400  km  is  somewhat 
greater  over  land  than  over  oceans.  However,  no  details 
were  given  in  that  paper  of  the  changes  in  the  shape 
or  amplitude  of  the  spectrum  as  a  function  of  topog¬ 
raphy.  or  of  any  dependence  of  the  relationship  on 
synoptic  weather  conditions.  One  purpose  of  the  pres¬ 
ent  paper  is  to  describe  more  thoroughly  the  changes 
in  the  mesoscale  spectra  as  functions  of  topography 
and  background  conditions  by  analysis  of  data  from 
GASP. 

The  interpretation  of  observed  mesoscale  spectra  has 
been  approached  from  two  fundamentally  different 
theoretical  bases  in  recent  years:  waves  and  turbulence. 
Waves  are  dispersive,  with  well-defined  relationships 
between  frequency  and  wavenumber  and  among  the 
amplitudes  of  atmospheric  variables,  while  turbulence 
is  a  strongly  nonlinear  process  with  strong  coupling 
between  different  scales.  There  is  no  question  that  well- 
defined  lee  waves  are  sometimes  present  over  moun¬ 
tainous  terrain,  and  it  seems  likely  that  small-amplitude 
bouyancy  waves  are  present  to  some  extent  under  most 
atmospheric  conditions  (Ecklund  et  al.,  1986).  It  has 
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been  suggested  that  a  background  spectrum  of  gravity 
waves  gives  rise  to  the  observed  spectrum  of  mesoscale 
variability  (Dewan.  1979:  VanZandt.  1982).  Recent 
studies  of  the  vertical  wavenumber  (m)  spectrum  of 
horizontal  velocities  show  that  the  observed  spectra 
appear  consistent  with  the  notion  of  a  saturated  gravity 
wave  spectrum  (Dewan  and  Good.  1986:  Smith  et  a!.. 
1985.  1 987 ).  On  the  other  hand,  certain  aspects  of  the 
observed  spectra  also  appear  consistent  with  the  con¬ 
cept  of  a  spectrum  of  quasi-two-dimensional  stratified 
turbulence  as  put  forward  by  Gage  (1979)  and  Lilly 
( 1983).  Aspects  of  agreement  are  the  ratio  of  spectral 
amplitudes  of  vertical  velocity  and  temperature  to  hor¬ 
izontal  velocity  and  the  apparent  success  of  the  Tay  lor 
hy  pothesis  for  transforming  frequency  and  wavenum¬ 
ber  spectra  (Gage  and  Nastrom.  1985.  1986).  The  point 
is  that  two  very  different  interpretations  can  be  offered 
for  the  observed  mesoscale  variations.  The  second  pur¬ 
pose  of  this  paper  is  to  examine  and  compare  the  im¬ 
plications  of  both  theories  in  light  of  the  present  ob¬ 
servational  results. 

The  data  and  method  of  data  analysis  are  described 
in  section  2.  Results  of  the  analysis  procedure  are  dis¬ 
cussed  in  section  3.  Sections  4  and  5  examine  the  in¬ 
terpretation  of  these  results  in  the  context  of  gravity 
wave  and  two-dimensional  turbulence  theory,  respec¬ 
tively.  Our  conclusions  are  presented  in  section  6. 

2.  Data  and  method 

All  aircraft  data  used  here  are  from  the  observations 
taken  during  G  ASP  and  are  a  subset  of  the  data  ana¬ 
lyzed  in  NG.  The  data  collection  phase  of  GASP  was 
during  1 975— 1979.  Instruments  for  the  measurement 
of  trace  gases  such  as  ozone  w  ere  placed  aboard  B747 
aircraft  operating  in  routine  commercial  service.  These 
and  other  variables  such  as  aircraft  location,  flight  level, 
pressure,  wind,  and  temperature  data  were  automati¬ 
cally  recorded  on  a  cassette  tape  during  flight  above 
6.1  km.  and  the  full  cassettes  were  returned  to  NASA 
for  quality  control  and  archiving.  Tropopause  pressures 
were  interpolated  in  time  and  space  from  gridded  NMC 
analyses  to  each  GASP  data  location  during  archiving. 
On  about  93  of  the  flights,  the  system  was  preset  to 
record  data  at  4  s  intervals  ( 1  km  intervals)  at  all  times 
above  6. 1  km  altitude.  Normally,  on  other  flights,  data 
were  recorded  at  5  min  intervals  (75  km  intervals  at 
250  m  s'1  airspeed).  However,  when  moderate  or 
greater  turbulence  was  encountered  (defined  as  when¬ 
ever  the  vertical  accelerometer  reading  fell  outside  the 
range  0.8- 1 .2  g)  the  system  recorded  data  at  4  s  interv  als 
from  the  beginning  of  the  turbulence  until  I  min  after 
the  last  “limit”  exceedence.  A  few  of  these  episodes 
lasted  for  up  to  several  minutes.  Other  details,  including 
a  map  of  the  geographic  distribution  of  GASP  data, 
are  given  in  NG  and  references  therein.  All  GASP  data 
are  available  from  the  National  Climatic  Center,  Ashe¬ 
ville. North  Carolina. 


One  goal  of  this  study  is  to  examine  mesoscale  vari¬ 
ability  over  regions  of  the  globe  with  vary  ing  surface 
topography  .  As  a  first  step  toward  that  goal,  we  have 
limited  attention  to  the  three  zones  shown  in  Fig  1  in 
the  latitude  belt  from  30°  to  50°N.  Zone  1  is  entirely 
over  the  Pacific  Ocean,  zone  2  includes  the  rough  ter¬ 
rain  of  the  western  United  States,  and  zone  3  lies  mostly 
over  the  central  plains.  The  Appalachian  Mountains 
extend  into  the  southeast  comer  of  zone  3.  but  since 
nearly  all  GASP  flights  to  the  Eastern  United  States 
terminated  in  New  York  or  Boston  the  flight  paths 
rarely  crossed  the  mountains. 

The  first  step  of  this  analysis  was  to  search  for  a 
geographic  bias  in  the  frequency  of  occurrence  of  tur¬ 
bulence  “limit"  episodes  and  of  the  mean  vertical  ac¬ 
celerometer  range  (the  range  over  the  previous  4  s  in¬ 
tervals  was  recorded  with  each  observation).  Most  case’s 
when  the  vertical  accelerometer  reading  fell  outside  the 
range  0.8  to  1 .2  g  are  within  a  100  km  or  so  of  airports: 
this  is  perhaps  not  surprising  as  it  is  in  the  approach 
region  that  pilots  have  the  least  freedom  to  choose  their 
flight  path  to  avoid  turbulence.  Away  from  airports, 
there  is  no  obvious  difference  in  the  frequency  of 
"limit"  episodes  among  the  three  zones.  This  result  is 
likely  an  indication  of  the  skill  of  aircrews  in  avoiding 
turbulence  rather  than  evidence  that  turbulence  levels 
are  the  same  over  oceans,  plains  and  mountains,  es¬ 
pecially  since  research  flights  (e.g..  Lilly  et  al..  1974) 
found  more  turbulence  over  rough  terrain.  We  con¬ 
clude  that  our  data  from  commercial  airliners  are  not 
useful  for  study  ing  the  frequency  of  turbulence  occur¬ 
rence. 

Variance  power  spectra  of  mesoscale  variations  of 
wind  and  temperature  were  computed  by  the  same 
methods  used  in  NG.  Briefly,  data  from  the  high-fre¬ 
quency  recording  ( I  km  data  interval)  flights  were  di¬ 
vided  into  segments.  A  segment  was  retained  for  anal¬ 
ysis  only  if  it  fell  entirely  within  a  zone,  did  not  cross 
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the  tropopause.  the  flight  level  did  not  van,  more  than 
50  m.  and  data  gaps  were  not  longer  than  a  prescribed 
limit.  The  analvsis  was  made  using  segments  64  or  256 
km  long.  On  the  64  km  (256  km)  segments,  the  data 
were  linearly  interpolated  to  1.25  km  (2.58  km)  inter¬ 
vals.  with  no  interpolation  allowed  over  gaps  greater 
than  3.5  km  (7.0  km).  The  mean  and  a  linear  trend 
were  removed  from  the  interpolated  data,  and  the 
spectrum  was  estimated  with  a  fast  Fourier  transform. 
The  results  were  stratified  by  geographic  zone.  b>  tro¬ 
posphere  and  stratosphere,  and  according  to  whether 
the  average  wind  speed  along  the  segment  was  greater 
or  less  than  25  m  s~  '. 

3.  Results 

Before  comparing  the  spectra  among  different  data 
stratification  groups,  it  is  interesting  to  compare  the 
total  variance  under  the  spectra  to  gain  a  first  impres¬ 
sion  of  topographic  effects.  Figure  2  show  s  the  variance 
values  in  a  bar  graph  format  for  zonal  (it)  and  merid¬ 
ional  (v)  wind  and  potential  temperature  (0)  for  each 
geographic  zone  and  separated  between  troposphere 
and  stratosphere  and  between  low-  and  high-back¬ 
ground  wind  speeds.  The  shaded  portion  of  each  bar 
represents  the  average  variance  on  the  64  km  segments, 
and  the  entire  bar  represents  the  average  variance  on 
256  km  segments.  The  number  of  values  averaged  in 
each  case  is  listed  in  Table  I  along  with  the  corre¬ 
sponding  mean  wind  speed,  and  the  average  variances 
are  listed  in  Table  2.  Comparing  zones  within  each 
data  group,  we  note  that  for  the  64  km  segments 
(shaded),  in  every  case  the  variance  in  zone  2  is  largest. 


and  zone  1  smallest.  For  zonal  w  ind,  the  ratio  of  zone 
2  to  zone  1  variance  ranges  from  4.3  in  the  stratosphere, 
high  wind  speed  cases,  to  2.5  in  the  troposphere,  low 
wind  speed  cases.  The  pattern  of  ratios  for  t  and  11  is 
similar.  Further,  the  same  pattern  among  zones  is 
found  for  the  variances  on  256  km  segments  except  in 
the  troposphere,  low  wind  cases,  for  r  and  rt. 

Comparison  of  variances  in  the  same  zones  shows 
that  for  u  and  t  from  the  64  km  segments  for  high 
wind  speeds,  there  is  little  difference  between  the  tro¬ 
posphere  and  stratosphere.  This  can  be  seen  in  Fig.  I. 
and  also  in  Table  2.  In  fact,  for  u  and  v  in  zones  2  and 
3.  the  only  noticeable  difference  is  between  the  tro- 


Txble  I.  Number  of  spectra  averaged  b>  zone  and  data  group 


Group 

Troposphenc  zone 

1  2  3 

Stratosphenc  zone 

1  2  3 

a  Case*  with  mean 

«  md  speed  less 

than  2 

5  m  v  1 

64-km  segments 

136 

130 

71 

49 

78 

120 

256-km  segments 

30 

25 

18 

12 

18 

25 

1200-km  segments 

75 

82 

46 

58 

76 

40 

Mean  wind  speed 

(m  s  ') 

14 

18 

19 

20 

18 

18 

b  Cases  with 

mean  m 

md  speed  greater  than 

25  ms' 

64-km  segments 

132 

162 

127 

42 

149 

166 

256-km  segments 

29 

35 

31 

9 

35 

37 

1200-km  segments 

68 

98 

48 

32 

85 

85 

Mean  ,nd  speed 

(m  s'1) 

37 

40 

44 

36 

36 

37 
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Table  2.  Average  variance  under  the  spectrum  tm*  s"  ’  and  k b>  geographic  zone 


Zone  1  (ocean  1 

Zone  2 
(mountains! 

Zone  3  (plains i 

Level 

Component 

LW 

HW 

LW 

HW 

LW 

HW 

Stratosphere 

u 

0  28 

64-km  segments 

0  34 

1  ()V 

1  45 

0.58 

0.52 

Troposphere 

u 

or 

0  34 

0  6V 

1.40 

0  34 

0  53 

Stratosphere 

l 

0  40 

0.46 

III 

1.20 

0  55 

0  56 

Troposphere 

r 

0.30 

0  32 

0  54 

1.00 

0  33 

0  52 

Stratosphere 

H 

0.2V 

0  31 

0  60 

1  05 

0  55 

0  41 

Troposphere 

H 

0.14 

0  IV 

0  23 

0  5  V 

UP 

0  2h 

Stratosphere 

u 

1.03 

256-km  segments 

LOO 

3. VI 

3.38 

1  P6 

-S 

Troposphere 

u 

1  P 

1  15 

1.61 

3.30 

1.15 

1  42 

Stratosphere 

l 

1  42 

13’ 

3.34 

3.25 

1.4! 

1  4  3 

Troposphere 

t 

0  V' 

1.0' 

1.22 

1P1 

I  fill 

Stratosphere 

H 

OP  3 

0  60 

1  60 

1  86 

110 

1  1 1 

Troposphere 

H 

0  60 

0.51 

0  46 

O.Vv 

0.34 

0  '2 

*  LW  Low  wind  speed  cases 
+  HW  High  wind  speed  cases 


pospheric  low  wind  speed  cases,  and  the  other  three 
data  groups.  In  zone  I.  the  differences  among  all  four 
data  groups  are  small.  For  temperature,  the  most  strik¬ 
ing  differences  among  data  groups  are  between  the  tro¬ 
posphere  and  stratosphere,  with  the  stratospheric  values 
a  factor  of  2.0  larger  on  the  average.  In  all  zones,  the 
smallest  temperature  results  are  found  in  the  tropo¬ 
spheric  low  wind  speed  cases.  In  zone  2.  the  result  for 
temperature  in  the  stratospheric  high  wind  speed  cases 
is  much  larger  than  that  for  other  data  groups  The 
comparisons  for  256  km  segment  results  show  similar, 
but  not  identical,  patterns  to  the  64  km  results,  but  the 
small  number  of  256  km  segments  give  reduced  sta¬ 
tistical  significance. 

There  are  sufficient  segments  of  the  64  km  results 
to  check  the  statistical  signficance  of  the  differences 
among  the  means,  and  to  make  other  statistical  anal¬ 
yses.  For  example,  the  frequency  distributions  of  vari¬ 
ances  of  zonal  wind  speed  along  segments  from  each 
geographic  zone  for  the  tropospheric  high  wind  cases 
are  compared  in  Fig.  3.  The  upper  panel  of  Fig.  3  shows 
the  curves  in  standard  histogram  format,  while  the 
lower  panel  shows  them  in  cumulative  probability  co¬ 
ordinates.  The  curves  appear  approximately  lognormal, 
consistent  with  the  findings  of  NG.  The  Student's  t- 
test  indicates  that  the  mean  values  for  zone  1  and  zone 
3  are  different  from  the  mean  for  zone  2  at  the  99 
percent  significance  level,  out  that  the  means  for  zones 
1  and  3  are  not  significantly  different  in  this  example. 
Results  for  other  data  groups  are  similar,  and  we  con¬ 
clude  that  there  is  significantly  more  variance  at  me- 
soscales  over  the  mountains  than  over  oceans  or  plains. 

In  order  to  assess  the  horizontal  wavelength  (A)  de¬ 
pendence  of  the  differences  among  variances,  it  is  nec¬ 
essary  to  examine  the  spectra.  Figures  4  and  5  show 
the  spectra  for  u  for  low  and  high  wind  conditions, 
respectively.  As  these  illustrate  the  main  points,  the 


spectra  of  r  and  0  are  not  shown.  The  coordinates  along 
the  ordinate  apply  to  the  upper  three  curves  which  are 
for  the  stratosphere,  and  the  three  troposphenc  curves 
are  shifted  down  one  decade.  Dashed  lines.  arbitranK 
chosen  with  slope  -5/3  and  offset  by  one  decade,  have 
been  added  to  aid  comparison.  At  wavelengths  less  than 
64  km.  the  curves  are  based  on  analyses  of  the  64  km 
segment  data.  At  longer  wavelengths,  they  are  based 
on  the  256  km  segments,  with  some  fairing  near  64 
km.  The  zone  number  is  entered  several  times  along 
each  curve.  The  plotted  values  are  the  averages  over  d 
samples,  with  d given  in  Table  1:  no  other  smoothing 
has  been  applied.  Error  bars  are  entered  at  two  wave¬ 
lengths  to  illustrate  the  uncertainty  of  the  means  and 
extend  2a(dy'/:  above  and  below  the  mean,  where  a 
is  the  standard  deviation. 

Figure  4.  for  low  wind  conditions,  shows  that  in  the 
troposphere  there  are  small  differences  among  zones 
at  A  <  4  km  and  A  >  80  km.  At  intermediate  A.  the 
curve  for  zone  2.  over  the  mountains,  reflects  signifi¬ 
cantly  increased  energy.  Zone  1,  over  the  ocean,  has 
the  lowest  energy  at  intermediate  A.  In  the  stratosphere, 
zone  1  has  the  lowest  energy  at  all  A  and  zone  2  the 
highest  energy  at  all  A  >  4  km.  For  high  wind  conditions 
(Fig.  5).  the  energy  in  zone  2  appears  high,  and  in  zone 
1  is  low.  at  all  A  >  4  km  in  both  the  troposphere  and 
stratosphere.  We  have  checked  the  continuation  of 
these  patterns  to  longer  wavelengths  using  the  1 200- 
km  segment  analyses  based  on  normal  recording-rate 
data  used  in  NG.  The  number  of  segments  is  given  in 
Table  1 .  These  results  (not  shown)  indicate  the  patterns 
seen  at  long  A  in  Figs.  4  and  5  persist  to  A  —  500  km. 
At  A  >  500  km.  the  flow  enters  the  enstrophy  cascading 
regime  associated  with  geostrophic  turbulence  follow¬ 
ing  a  A-3  spectral  shape. 

The  variance  and  spectral  analysis  results  presented 
here  show  there  are  significant  differences  in  mesoscale 
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variability  related  to  surface  topographs.  The  differ¬ 
ences  range  up  to  a  factor  of  6  in  power,  and  are  largest 
in  the  scale  ra  ige  10  <  \  <  80  km.  A  similar  factor 
was  noted  by  Lilly  et  al.  ( 1974)  between  variances  over 
ocean  and  mountains  at  scales  less  than  610  m.  They 
interpreted  their  observations  w  ithin  the  framework  of 
three-dimensional  turbulence  theory.  The  present  re¬ 
sults  are  examined  using  the  theories  of  internal  gravity 
waves  and  quasi-two-dimensional  stratified  turbulence 
in  the  following  sections. 

4.  Internal  gravity  wave  interpretation 

We  begin  our  discussion  in  this  section  by  postulat¬ 
ing  that  the  motion  and  temperature  spectra  presented 


previously  are  due  to  internal  gravity  waves  and  that 
these  waves  comprise  a  spectrum  in  which  the  ampli¬ 
tudes  of  the  smallest  vertical  scales  are  limited  by  sat¬ 
uration  processes  (VanZandt,  1982;  Dewan  and  Good. 
19b6;  Smith  et  ai..  1987).  The  effects  of — and  processes 
contributing  to — gravity  wave  saturation  have  been 
reviewed  by  Fritts ( 1984)and  Frittsand  Rastogi  ( 1985) 
and  are  now  believed  to  contribute  to  the  large-scale 
circulation  and  thermal  structure  throughout  the  at¬ 
mosphere. 

We  assume  a  saturated  vertical  wavenumber  spec¬ 
trum  of  the  form  (Smith  et  al..  1987) 


E(m) 


1 

I  + 


( 1 1 


such  that  the  maximum  energy  weighted  by  wave- 
number  (mEim))  occurs  at  a  wavenumber  m  -  mt. 
21  '  w  here  mt  is  the  characteristic  vertical  wavenumber 
Typical  values  of  inferred  in  the  upper  tropo¬ 

sphere  and  lower  stratosphere  are  —2-2.5  km  (Fritts 
and  Chou.  1987;  Fntts  et  al..  1987).  We  then  assume 
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Flo  4  Average  spectra  of  zonal  wind  speed  for  cases  of  low  -back¬ 
ground  winds  Geographic  zone  number  is  entered  along  each  curve 
Values  on  ordinate  apply  to  stratosphenc  curves,  tropospheric  curves 
are  shifted  down  one  decade  (See  text  ) 
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that  this  spectrum  responds  to  increases  in  .V  following 
Vv'KB  scaling  and  relaxation  to  the  new  saturated  am¬ 
plitudes  as  anticipated  by  VanZandt  and  Friits  ( 1982). 
This  causes  the  saturated  spectrum  to  evolve  as 

w*  ~  A'2  \  (2) 

Thus  a  characteristic  velocity,  assuming  hydrostatic 
two-dimensional  wave  motions,  is 

X  2|,.V 

u  ~~(c-u)  =  —  ~ — —  —  A  .  (3) 

m 

Consistent  with  (1),  we  have  assumed  here  that  the 
dominant  \ertical  scale  is  near  saturation  amplitude. 
Also,  from  the  linear  thermodynamic  equation, 

T 

(4) 

implying  that  temperature  fluctuations  should  be  more 
sensitive  than  wind  fluctuations  to  changes  in  stratifi¬ 
cation.  Note  that  this  saturation  theory  leads  to  a  dif¬ 
ferent  dependence  on  X  than  the  usual  WKB  scaling 
of  the  wave  spectrum  fie.,  u'  ~  A’1'2). 

Comparing  these  scalings  with  our  observations,  we 


find  the  best  agreement  to  occur  for  the  average  ratios 
of  stratospheric  to  tropospheric  wind  variances  under 
low  wind  conditions.  These  are  1.62  and  1.64  for  the 
64  km  and  256  km  segments,  respectively,  and  com¬ 
pare  very  well  with  the  value  anticipated  from  (3)  for 
a  change  in  X  of  2  (22/;i  =  1.59).  The  ratios  are  some¬ 
what  smaller  for  high  wind  conditions  (near  1.2).  per- 
hans  reflecting  a  tendency  for  higher  wind  speeds  to 
allow  for  greater  reflection  or  dissipation  of  wave  energy 
near  the  tropopause.  depending  on  the  direction  of 
wave  propagation. 

Like  the  velocity  variances,  the  ratios  of  stratospheric 
to  tropospheric  temperature  variances  are  most  con¬ 
sistent  with  the  foregoing  scalings  for  low  wind  speeds. 
The  64  and  256  km  segment  averages  are  2.67  and 
2.45.  with  the  largest  individual  values  (from  2.61  to 
3.48)  occurring  in  zones  2  and  3  and  the  smallest  values 
(1.22  and  2.07)  in  zone  I.  These  are  approximately 
half  of  that  expected  (2267  =  6.36).  suggesting  that  an¬ 
other  process,  such  as  convection,  rather  than  wave 
activity  might  have  contributed  appreciably  to  tem¬ 
perature  variance  in  the  troposphere. 

We  now  consider  the  implications  of  the  differences 
in  wind  and  temperature  variance  in  the  three  zones 
for  the  gravity  wave  spectrum  and  for  wave  propagation 
and  effects  in  the  lower  and  middle  atmosphere.  As 
noted  in  section  3.  the  variances,  particularly  for  the 
64  km  data  segments,  exhibit  statistically  significant 
differences  between  zones.  Variances  are  much  larger, 
on  average,  in  zone  2  than  in  zones  1  and  3.  and  those 
in  zone  3  are  somewhat  larger  than  in  zone  I.  The 
zonal  velocity  spectra  presented  in  Figs.  4  and  5  reveal 
a  much  more  energetic  motion  field  at  horizontal  scales 
from  —  4  to  80  km  in  zone  2.  suggesting  that  these 
scales  may  be  preferentially  excited  over  rough  terrain. 
The  larger  characteristic  velocities  imply  a  smaller  value 
of  mt  (and  larger  characteristic  vertical  wavelength) 
from  (3)  at  a  given  horizontal  scale,  corresponding  to 
wave  motions  with  larger  intrinsic  frequencies.  Thus, 
on  the  basis  of  the  velocity  variance  data,  we  expect 
that  the  wave  spectrum  in  zone  2.  relative  to  that  in 
zones  1  and  3.  is  composed  of  wave  motions  with  larger 
vertical  velocities,  steeper  propagation  angles,  and  a 
greater  potential  for  large  vertical  fluxes  of  energy  and 
horizontal  momentum.  mV  and  r  V. 

To  quantify  these  observations,  we  assume  that  the 
velocity  variance  is  the  best  indication  of  the  relative 
level  of  wave  activity  and  note  that  a  representative 
value  in  zone  2  is  twice  that  in  zone  I  or  3.  Tnen  a 
characteristic  velocity  or  intrinsic  frequency  is  also 
larger  by  — 1.4  in  zone  2.  This  causes  the  energy  and 
momentum  fluxes. 

c^E-^Ac-u?  (5) 

and 

mV  -  Vc-  u)2 
A 


(6) 
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to  increase  in  /one  2  by  ^4  and  2.8.  respectively .  rel¬ 
ative  to  those  in  /ones  1  and  3.  assuming  that  the  wave 
spectrum  in  each  zone  has  the  same  degree  of  aniso¬ 
tropy  If  the  majority  of  the  enhanced  wave  activity  in 
/one  2  can  be  attributed  to  topography,  however,  the 
potential  momentum  flux  and  applied  drag  in  this  re¬ 
gion  may  be  even  larger  due  to  a  greater  anisotropv  of 
the  wave  held  anu  because  the  greater  variances  in  /one 
2.  if  applicable  to  a  smaller  mountainous  area,  miplv 
even  greater  enhancements  of  the  energy  and  momen¬ 
tum  fluxes  in  these  regions 

These  results  are  significant  because  they  suggest  that 
rough  terrain  may  account  for  a  disproportionutelv 
large  traction  of  the  total  momentum  flux  that  aflects 
the  atmosphere  at  greater  heights  due  »o  the  larger  wave 
amplitudes  and  intrinsic  frequencies  that  arise  in  these 
regions  The  results  also  support  the  suggestion  by  Pal¬ 
mer  et  al  I  I4SM  that  major  topographs  might  be  re¬ 
sponsible  for  the  drag  that  is  needed  to  achieve  realistic 
mean  umds  in  numerical  simulations  of  the  lower  at¬ 
mosphere  A  related  point  is  that  due  to  the  higher 
intrinsic  frequencies  and  essentially  vertical  group  ve¬ 
locities  that  appear  to  be  associated  with  topograph¬ 
ically  generated  waves,  the  drag  and  induced  diffusion 
due  to  these  wave  motions  should  occur  in  the  middle 
atmosphere  primarily  in  those  regions  that  overlie  ma¬ 
jor  topography . 

The  interpretation  of  the  GASP  data  offered  in  this 
section  was  based  on  the  assumption  that  the  observed 
motion  held  is  due  primarily  to  internal  gravity  waves. 
This  interpretation  will  be  compared  and  contrasted 
with  that  offered  in  terms  of  quasi-two-dimensional 
turbulence  in  section  6. 

5.  Wave  breaking,  small-scale  turbulence  and  meso- 
scale  stratified  turbulence 

In  section  4  it  was  assumed  that  the  GASP  spectra 
presented  in  section  3  were  due  solely  to  a  spectrum 
of  internal  waves.  In  this  section  we  consider  the  more 
general  condition  when  the  spectra  are  comprised  of  a 
combination  of  waves  and  turbulence.  Specifically,  we 
will  consider  what  happens  when  enhanced  wave  ac¬ 
tivity  over  mountainous  terrain  leads  to  small-scale 
turbulence  generation.  Under  these  circumstances  it  is 
anticipated  that  some  of  the  small-scale  turbulent  ki¬ 
netic  energy  participates  in  an  inverse  cascade  to  larger 
scales  as  suggested  by  Gage  ( 1979)  and  Lilly  ( 1983)  for 
quasi-two-dimensional  stratified  turbulence.  According 
to  this  scenario  a  field  of  motion  comprised  of  waves 
and  turbulence  can  evolve  into  a  field  of  quasi-hon- 
zontal  eddy  motions  as  w  ave  energy  radiates  out  of  the 
fluid  volume.  The  residual  field  of  quasi-horizontal  ed¬ 
dies  is  the  vorticity -bearing  mode  of  the  motion  field. 

The  possibility  of  a  vorticity -bearing  mode  of  quasi- 
horizonta!  motions  in  the  atmospheric  mesoscale  fol¬ 
lows  directly  from  the  RMW  scaling  of  the  equations 
of  motion  (Riley  et  ah.  1981:  Lilly.  1983).  Briefly,  at 


low  Froude  numbers  pertinent  to  the  stably  stratified 
atmosphere,  the  equations  of  motion  decouple  into  a 
set  of  equations  governing  internal  wave  motions  and 
a  set  of  equations  governing  stratified  turbulence 
Stratified  turbulence  as  described  by  Lilly  ( 1983 1  is  the 
vorticity -bearing  mode  of  motion.  The  possibility  that 
the  vonical  mode  of  motion  exists  in  the  ocean  has 
been  raised  by  Muller  ( 1984)  and  Muller  et  al  i  198m 
as  a  way  to  account  for  departures  from  the  Garret- 
Munk  internal  wave  spectrum  observed  during  IWI  \ 
(Muller  et  ah.  IP’S). 

The  basic  idea  of  stratified  turbulence  has  been  dis¬ 
cussed  by  Gage  <|9~9i  and  Lilly  (1983).  Briefly,  n 
inverse  cascade  is  hypothesized  to  redistribute  turbulent 
kinetic  energy  throughout  the  atmospheric  mesoscale 
from  a  small-scale  source  Lilly  (1983)  has  described 
the  source  as  upscale  leakage  from  three-dimensional 
turbulence.  He  found  that  only  a  few  percent  of  the 
energy  present  in  three-dimensional  turbulence  is 
needed  to  account  for  the  observed  level  of  (he  me¬ 
soscale  atmospheric  horizontal  velocity  spectrum 

The  range  of  horizontal  scales  pertinent  to  a  tran¬ 
sition  between  microscale  three-dimensional  turbu¬ 
lence  and  mesoscale  stratified  turbulence  can  be  defined 
as  follows.  The  smallest  scale  for  this  transition  /one 
can  be  identified  with  the  buoyancy  length  scale 


w  here  c  is  the  dissipation  rate  of  turbulent  kinetic  en¬ 
ergy  and  A  the  Brunt-Yaisala  frequency.  In  the  at¬ 
mosphere.  the  buoyancy -length  scale  can  vary  from  10 
m.  or  less,  in  very  stable  regions  of  the  stratosphere  to 
300  m.  or  more,  in  turbulent  regions  of  the  troposphere 
For  scales  L  <  Lb  the  dynamics  is  fairly  well  de¬ 
scribed  by  the  familiar  Kolmogoroff  turbulence  theory 
Consistent  with  this  theory,  the  kinetic  energy  spectrum 
is  described  by 

£(A )  =  .4  ,r'  k  .  k  >  k#  =  -  (8) 

'■H 

where  k  is  horizontal  wavenumber  and  .  f ,  is  a  univervil 
constant. 

The  energy  source  for  the  microscale  turbulence  in 
the  stable  atmosphere  is  though;  to  be  breaking  internal 
waves,  although  other  sources  are  possible.  Presumably . 
the  energy  content  of  stratified  turbulence  should  be 
largest  in  regions  of  enhanced  wave  breaking 

The  stratified  turbulence  has  a  spectrum  of  the  form 

£<A)  =  .4:|^j  k  k  <  A„  (9) 

where  is  the  largest  wavenumber  in  the  inertial  range 
of  stratified  turbulence.  dEjdi  the  rate  of  energy  inser¬ 
tion  into  the  inverse  cascade  and  .4;  the  universal  con¬ 
stant  for  stratified  turbulence.  ,4;  —  9.4,  (Lilly.  1983) 
The  scale  of  energy  insertion  k,  =  2v/L,  is  presumably 
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related  to  the  vertical  scale  of  wave  breaking  and  is 
necessarily  in  the  range 

ke  >  k,  >  ku .  (10) 

The  spectral  forms  given  in  Eqs.  (8)  and  (9)  both 
represent  idealized  equilibrium  model  spectra  for  in¬ 
ertial  ranges  in  which  the  energy  transfer  up  or  down 
the  spectrum  is  independent  of  wavenumber.  At  in¬ 
termediate  scales,  i.e..  those  that  occupy  the  range  of 
wavenumber  space  adjacent  to  k,.  the  energy  cascade 
rate  must  be  dependent  on  k. 

Weinstock  (1980)  considered  a  spectral  model  in 
which  the  energy  cascade  rate  is  wavenumber  depen¬ 
dent.  He  found  a  tendency  for  a  spectral  hump  to  form 
at  length  scales  smaller  than  the  source  scale  and  for  a 
gap  to  form  at  length  scales  associated  with  the  source 
scale.  At  length  scales  larger  than  the  scale  of  the  source, 
there  is  thus  a  wavenumber  domain  of  steepened  (more 
negative)  spectral  slope.  Similar  issues  have  been  dis¬ 
cussed  by  Lilly  (1983)  in  the  context  of  an  intermittent 
source  of  turbulent  kinetic  energy  associated  with  con¬ 
vection.  On  the  other  hand,  while  Lilly's  model  shares 
some  of  the  features  of  the  Weinstock  model,  such  as 
the  inverse  cascade,  there  is  no  clear  connection  be¬ 
tween  the  Weinstock  and  Lilly  models. 

The  physical  mechanisms  associated  with  wave 
breaking  in  the  atmosphere  have  been  discussed  at 
length  in  the  literature.  (See.  for  example.  Bretherton. 
1969:  Fritts  and  Rastogi.  1985;  Dewan  and  Good. 
1986.)  Usually  ,  wave  breaking  is  treated  as  a  manifes¬ 
tation  of  a  local  instability,  such  as  convective  over¬ 
turning  or  Kelvin-Helmholtz  instability.  The  turbu¬ 
lence  that  results  from  wave-breaking  episodes  is  often 
referred  to  as  wave-induced  turbulence  since  it  tends 
to  occur  when  a  local  instability  is  produced  in  certain 
phases  of  an  atmospheric  wave.  The  consequence  of 
wave  breaking  and  turbulence  generation  on  atmo¬ 
spheric  spectra  of  stratihed  turbulence  is  presumably 
to  create  a  transition  region  spanning  a  decade  or  so 
in  wavenumber  space  at  scales  larger  than  the  buoy¬ 
ancy-length  scale.  Presumably,  the  most  intense  wave- 
breaking  episodes  will  be  associated  with  instabilities 
on  the  largest  vertical  wave  scales  and  with  the  largest 
turbulence  production.  Consequently,  under  condi¬ 
tions  of  enhanced  wave-breaking  activity  we  miy  an¬ 
ticipate.  at  least  qualitatively,  steepening  of  spectral 
slopes  at  larger  scales,  and  an  increased  energy  level  of 
mesoscale  spectra  since  dE/dt  should  be  proportional 
to  c. 

An  examination  of  the  G  ASP  spectra  presented  in 
Fig.  4  and  Fig.  5  reveals  the  presence  of  enhanced  spec¬ 
tral  energy  levels  and  steepened  spectral  slopes  as  pre¬ 
viously  discussed.  To  illustrate  these  features.  Fig.  6 
contains  the  GASP  spectra  over  mountainous  terrain 
under  high  wind  conditions.  In  terms  of  the  previous 
discussion,  the  vertical  arrows  indicate  the  spatial 
scale  Lu  «  (2 r/ku)  associated  with  the  transition  from 
a  AT 5/3  spectral  slope  at  low  wavenumber  to  a  steepened 


Zonal  Wind  Spectra 
High  Wind  Speed  Cases 


Fig  6  As  in  Fig  5.  excepi  only  for  geographic  rone  2 
and  wnh  spectral  breaks  noted 


spectral  slope  at  high  wavenumber.  The  scale  of  the 
source  region  lies  to  the  right  of  the  vertical  arrows  in 
Fig.  6  probably  in  the  1-3  km  range  and  not  ade¬ 
quately  resolved  in  the  GASP  spectra.  In  Fig.  6  the 
transition  scale  occurs  at  about  15  km  in  the  strato¬ 
sphere  and  about  23  km  in  the  troposphere.  In  the 
A' 5/3  region  at  larger  scales  the  energy  levels  are  about 
the  same  in  the  troposphere  and  stratosphere.  These 
energy  levels  are  more  than  half  an  order  of  magnitude 
larger  than  the  more  ty  pical  spectral  amplitude  evident 
in  the  G  ASP  spectra  (Nastrom  et  al..  1984;  Nastrom 
and  Gage.  1985)  as  shown  in  Fig.  7. 

Under  the  assumptions  made  earlier  in  this  section, 
regional  differences  in  eddy  dissipation  rates  appear  to 
reflect  regional  differences  in  wave-breaking  intensity  . 
If  the  energy  insertion  rate  dE/dl  also  reflects  these 
regional  differences,  we  would  anticipate  that  the  spec¬ 
tral  amplitude  in  the  k~ 5/3  range  of  the  mesoscale  spec¬ 
trum  would  scale  as  r2/\  At  least  in  the  stratosphere, 
eddy  dissipation  rates  pertinent  to  the  GASP  data  are 
probably  about  1  cm2  s" 3  over  the  mountainous  region 
and  0  1  cm2  s"3  over  the  ocean  (cf.  Lilly  et  al..  1974). 
We  anticipate,  accordingly,  about  a  factor  of  5  increase 
of  mesoscale  spectral  amplitude  over  the  mountainous 
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Fig  1  C ompanson  of  GASP  zonal  wind  spectrum  o\er  moun¬ 
tainous  terrain  with  the  climatological  mean  GASP  spectrum  (Nas- 
trom  and  Gage.  1985 1  and  the  spectrum  of  Lill>  and  Lester  ( 1974) 


terrain  compared  to  mesoscale  spectral  amplitude  over 
the  oceans.  The  observed  ratio,  from  Figs.  4  and  5,  is 
about  3  or  4  on  the  average,  and  ranges  up  to  6.  As 
might  have  been  anticipated,  these  effects  seem  to  be 
most  pronounced  with  strong  winds.  Comparable  cli¬ 
matological  studies  of  turbulence  intensity  in  the  tro¬ 
posphere  have  yet  to  be  made. 

The  higher  amplitudes  of  GASP  spectra  over  the 
plains  compared  to  the  ocean  may  be  the  result  of  east¬ 
ward  advection  of  decaying  stratified  turbulence  or  of 
the  injection  of  turbulent  energy  by  enhanced  convec¬ 
tive  activity.  Since  stratified  turbulence  should  be  most 
persistent  in  the  stratosphere,  relatively  high  spectral 
levels  may  be  anticipated  in  the  stratosphere  over  the 
plains  even  under  light  wind  conditions.  Indeed,  under 
such  conditions  stratospheric  spectra  are  relatively 
energetic  compared  to  the  tropospheric  spectra  ob¬ 
served  over  the  plains. 

In  comparing  the  results  presented  here  with  earlier 
work  on  mountain  waves,  it  is  important  to  recognize 
that  the  GASP  flights  avoid  areas  of  intense  wave  ac¬ 
tivity  and  turbulence.  Thus,  spectra  observed  under 
very  active  conditions  should  be  more  energetic  than 
the  GASP  spectra.  The  spectra  reported  by  Lilly  and 
Lester  (1974)  for  example,  under  active  conditions 
appear  to  be  about  an  order  of  magnitude  more  intense 
than  those  reported  here  (Fig.  7).  Lilly  and  Lester  (1974) 
estimate  the  eddy  dissipation  rate  to  be  about  30  cm2 
s'  5  for  these  distributed  conditions.  An  order  of  mag¬ 
nitude  increase  of  mesoscale  spectral  amplitude  above 
the  GASP  spectra!  amplitude  is  consistent  with  the 
foregoing  arguments. 

In  this  section  we  have  considered  the  consequences 
of  an  evolving  field  of  motion  comprised  of  waves  and 
turbulence.  We  have  found  that  many  of  the  features 
of  the  GASP  spectra  reported  in  section  3  are  consistent 
with  the  idea  that  these  spectra  have  evolved  from  more 
active  conditions  in  which  the  wave  breaking  and  tur¬ 


bulence  generation  leads  to  an  inverse  cascade  and  a 
persistent  field  of  stratified  turbulence.  We  have  hy¬ 
pothesized  that  the  dynamics  of  the  disturbed  flow  over 
mountainous  terrain  may  serve  as  a  significant  regional 
source  for  the  GASP  spectra  reported  in  Nastrom  and 
Gage  (1985). 

6.  Conclusions 

Wind  and  temperature  measurements  taken  at  air¬ 
line  cruise  altitudes  during  GASP  have  been  analyzed 
over  three  geographic  regions:  oceans,  mountains,  and 
plains.  The  resuits  show  that  underlying  terrain  has  a 
significant  effect  on  the  atmospheric  variability  at  scales 
from  about  4  to  80  km.  w  ith  variances  over  mountains 
up  to  six  times  larger  than  those  over  oceans.  The  en¬ 
hanced  energy  over  mountains  is  found  in  both  wind 
and  temperature  results  and  is  persistent,  being  found 
in  the  troposphere  and  stratosphere  during  both  low 
and  high  background  wind  speed  conditions.  At  longer 
wavelengths.  80  to  500  km.  the  variances  are  also  larger 
over  mountains  than  in  other  regions  in  the  strato¬ 
sphere  and.  in  most  cases,  in  the  troposphere. 

These  results  complement  past  observational  studies 
such  as  Lilly  et  al.  (1974)  and  Pinus(  1979).  The  GASP 
data  are  unique  in  that  they  are  sufficiently  copious  to 
support  comparisons  of  mean  values  for  various  back¬ 
ground  conditions,  such  as  high  and  low-background 
wind  speed  and  cover  larger  spatial  scales  than  other 
datasets.  Also,  the  large  number  of  spectra  available 
has  permitted  us  to  establish  the  lognormal  distribution 
of  variances  regardless  of  underlying  terrain.  Finally, 
the  capability  to  subdivide  results  has  allowed  some 
comparisons  with  predictions  derived  from  the  gravity 
wave  and  stratified  turbulence  theories  of  mesoscale 
variability. 

The  linear  theory  of  gravity  wave  motions  was  used 
to  predict  the  scaling  of  wave  amplitude  with  back¬ 
ground  stability  both  for  wind  and  temperature.  Ex¬ 
cellent  agreement  of  the  predicted  scaling  was  found 
for  the  wind  energy  in  the  case  of  low-background  wind 
speed.  The  comparison  of  the  predicted  scalings  was 
less  satisfactory  for  wind  energy  in  the  case  of  high- 
background  wind  speed,  and  for  temperature  for  all 
background  wind  speed*.  It  was  noted  that  other  pro¬ 
cesses  such  as  wave  reflections,  wave  dissipation,  con¬ 
vection.  etc.,  which  were  not  accounted  for  here  may 
be  important  in  reconciling  the  wave  theory  with  the 
observations. 

The  theory  of  stratified  turbulence  was  used  to  pre¬ 
dict  an  excess  of  spectral  energy  over  mountainous  re¬ 
gions  in  the  mesoscale  AT 5/3  region  and  a  steepened 
spectral  slope  at  higher  wavenumbers.  The  observed 
spectral  amplitudes  in  the  mesoscale  AT 5/3  region  vary 
among  geographic  regions  approximately  in  accord 
with  an  c2/3  scaling.  Also,  the  notion  of  advection  of 
turbulence  appears  to  be  supported  by  the  enhanced 
spectra!  amplitudes  found  downwind  of  the  mountains 
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although  no  quantitative  predictions  were  tested,  and 
other  source  mechanisms,  such  as  convection,  are  pos¬ 
sible  in  such  areas. 

In  summary,  the  theories  of  internal  gravity  waves 
and  stratified  turbulence  each  account  for  some,  but 
not  all.  aspects  of  the  observations  as  analyzed  in  this 
brief  studs.  While  there  is  no  question  that  wave  and 
turbulence  processes  both  contribute  to  the  observed 
mesoscale  spectra,  we  have  made  only  a  small  effort 
to  synthesize  the  two  theories,  ft  appears  that  a  more 
detailed  investigation  would  have  mem 

Resolution  of  the  question  of  interpretation  is  of 
great  practical  significance  as  well  as  of  theoretical  in¬ 
terest.  If  the  linear  gravity  wave  theory  can  be  correctly 
applied  to  the  GASP  data,  results  such  as  those  given 
here  could  be  extended  to  infer  important  quantities 
such  as  momentum  and  energy  transfer.  In  addition 
to  the  linear  theory  progress  is  being  made  in  devel¬ 
oping  nonlinear  wave  theory  (e.g..  Fritts.  1984.  Wein- 
stock.  1985:  Muller  et  al..  1986).  The  potential  im¬ 
portance  of  the  vortical  mode  and  stratified  turbulence 
has  onlv  recently  been  appreciated.  Here.  too.  progress 
is  being  made  (e  g..  Metais  and  Herring.  1985:  Herring 
and  Metais.  1987).  A  major  challenge  for  the  future  is 
to  understand  how  each  of  these  processes  contributes 
to  the  observ  „d  mesoscale  spectra  and  how  they  interact 
with  one  another. 
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1.  Introduction 


Attempts  to  observe  the  vertical  motion  field  at  synoptic 
and  larger  scales  with  MST  radars  have  been  frustrated  by  the 
presence  of  geophysical'  "noise"  in  the  data,  consisting  of 
contributions  from  small  scales,  which  usually  overwhelms  any 
large-scale  signal  in  the  data.  Two  recent  such  studies  are 
Nastrom,  et  al.  (1985),  and  Larsen,  et  al.  (1988),  where  more 
detailed  discussion  of  the  background  and  references  to  other 
work  may  be  found.  Briefly,  all  past  studies  were  performed  in 
regions  where  the  noise  produced  by  orographic  effects  could  not 
usually  be  ignored,  and  the  possibility  existed  that  a  station 
well-removed  from  such  effects  might  be  able  to  observe  the  large 
scale  components  in  the  wind  field  under  more  general  conditions. 
Ongoing  studies  at  the  Flatland  radar  in  central  Illinois,  USA, 
are  testing  this  hypothesis,  and  an  emerging  result  is  that  the 
vertical  motion  field  is  still  dominated  by  scales  smaller  than 
synoptic  even  though  orographic  effects  are  evidently  not 
present.  This  effect  is  clearly  reflected  in  the  power  spectrum 
of  observed  vertical  motions  (Figure  1)  as  there  is  much  more 
energy  at  periods  less  than  an  hour  or  so  than  at  the  longer 
periods  commonly  associated  with  synoptic-scale  activity.  While 
the  high  frequency  variations  in  vertical  velocity  are  apparently 
due  to  a  field  of  propagating  gravity  waves  which  are  Doppler 
shifted  by  the  ambient  horizontal  wind  f ield_(VanZandt,  et  al., 
1989)  ,  we  expect  that  much  of  the  energy  at  longer  periods 
represents  contributions  from  large-scale  motion  systems  such  as 
baroclinic  storms.  However,  it  is  not  obvious  that  temporal 
noise  at  a  single  station  can  in  general  be  filtered  to  provide 
mean  values  consistent  with  a  large  spatial  average  (i.e., 
application  of  the  Taylor  hypothesis) .  Consequently,  the 
observation  of  synoptic  scale  vertical  motions,  believed  to  be  of 
the  order  of  a  few  cm  s“  ,  may  not  be  straightforward  even  at 
Flatland.  The  removal  of  orographic  effects  from  the  signal  has 
simplified  the  problem  to  the  extent  that  the  nature  of  the  other 
noise  sources  may  now  be  more  clearly  identified.  In  this  study 
we  explore  some  preliminary  aspects  of  these  problems:  first,  we 
show  that  periods  of  enhanced  variance  of  the  vertical  velocity 
during  the  interval  studied  are  associated  with  the  passage  of 


fronts  at  the  surface,  and  idenvify  a  frontal  signature  that  is 
sometimes  present  in  the  lower  troposphere.  Second,  ve 
illustrate  the  difficulties  in  verification  of  the  Taylor 
hypothesis  by  comparing  radar  observed  vertical  motions  with 
those  deduced  from  the  horizontal  wind  field  derived  from  the 
twice  daily  NWS  rawinsonde  network,  and  discuss  the  problems 
associated  with  such  a  comparison. 

2.  Analysis 

Time  series  of  vertical  velocity  over  Flat land  during  March 
1987  are  shown  in  Figure  2.  As  at  other  locations  (e.g. , 

Ecklund,  et  al.,  1982),  the  presence  of  quiet  and  active  periods 
is  clearly  visible  ,  though  the  variance  in  all  cases  is  less  at 
Flatland  than  that  at  mountain  sites  (Green,  et  al.,  1988). 

Along  the  bottom  of  Figure  2  letters  have  been  entered  to 
indicate  when  cold  (c) ,  warm  (w) ,  occluded  (o) ,  or  quasi¬ 
stationary  (qs)  fronts  passed  or  were  in  the  vicinity  of 
Flatland.  Note  that  every  active  period  is  coupled  with  a  front. 
While  this  coupling  suggests  the  interesting  notion  that  fronts 
are  a  source  of  gravity  waves,  it  hinders  the  prospects  for 
retrieving  large-scale  vertical  motions  from  the  Flatland  data 
because  the  geophysical  "noise"  apparently  increases  just  at 
those  times  when  we  hope  to  find  a  large-scale  signal. 

Next  we^.will  examine  two  cases  in  detail  to  illustrate  two 
points:  first,  there  is  a  frontal  passage  "signature"  which  is 

found  in  the  vertical  velocity  time  series  in  some  cases;  second, 
matching  the  signal  at  one  station  with  a  system  which  may  be 
moving  at  an  unsteady  rate  (as  fronts  and  storms  often  move)  and 
which  may  be  decaying  or  developing  at  a  rate  which  is  only 
poorly  defined  by  the  12-hourly  radiosonde  observation  schedule 
is  an  extremely  complex  problem.  Weather  systems  are  fully 
three-dimensional  phenomena  which  move  and  change  with  time.  The 
Flatland  radar  provides  a  continuous  record  of  the  vertical 
motions  over  a  single  place.  The  morphology  of  the  low-frequency 
signal  in  vertical  velocity  varies  from  case  to  case,  depending 
upon  the  rate  and  consistency  of  the  movement  of  the  weather 
system  and  upon  its  rate  of  development  or  stage  of  maturity.  In 
some  cases  the  vertical  velocity  pattern  seen  by  the  radar  is 
relatively  clear  and  the  synoptic  feature  easily  discerned,  while 
in  other  cases  it  is  less  well  defined.  We  illustrate  each  of 
these  points  with  specific  examples. 

First,  we  consider  a  relatively  clear  frontal  passage 
signature.  On  5  March  a  weak  low  was  over  western  Illinois  at  06 
UTC  with  a  warm  front  moving  northward  across  the  Flatland  area. 
The  front  passed  near  05  UTC  based  on  surface  observations  at 
Champaign  airport,  although  the  National  Weather  Service  analysis 
indicated  it  was  weak  and  nearly  dissipated  by  12  UTC  (Figure  3)  . 
The  time  series  of  vertical  velocities  display  upward  motion 
ahead  of  the  surface  frontal  passage  and  downward  motions  behind 
it  below  5  km,  as  seen  in  Figure  2  and  as  shown  more  clearly  by 
the  analysis  of  hourly  mean  velocities  in  Figure  4.  This  pattern 
was  noticed  by  Larsen  and  Rottger  (1982)  in  connection  with  the 
passage  of  a  warm  front,  and  we  have  seen  it  in  other  cases  of 


moving  fronts.  While  it  seems  unlikely  that  the  regions  of 
upward  and  downward  motion,  located  only  about  4  hours  (150  km  if 
the  front  moved  at  10  m  s-1  )  apart,  are  due  to  the  front  itself 
as  envisioned  in  the  Norwegian  frontal  model,  we  suggest  it  is  a 
typical  frontal  signature.  Perhaps  it  is  due  to  a  gravity  wave 
launched  by  the  front  and  propagating  along  with  it,  similar  to 
the  events  discussed  by  Gall,  et  al.(1988). 

A  front  is  in  many  respects  a  mesoscale  feature,  yet  in 
order  to  compute  vertical  motions  aloft  from  conventional  data  we 
must  rely  on  the  radiosonde  network  which  reflects  only  larger 
scales.  For  example,'  around  Flatland  the  triangle  of  radiosonde 
stations  at  Peoria,  Salem  and  Dayton  has  an  average  leg-length  of 
379  km;  which  defines  the  smallest  horizontal  scale  observable 
with  this  network  over  Flatland.  The  vertical  motions  computed 
using  the  continuity  equation  applied  to  the  rawinsonde  data  over 
this  triangle  for  5/00  and  6/00  UTC;  using  Ekman  pumping  as  a 
lower  boundary  condition  following  Nastrom,  et  al.(1985),  are 
compared  with  radar  time-averages  in  Figure  5.  The  radar  data 
represent  6-hour  averages  and  the  error  bars  are  the  standard 
error  of  the  mean.  At  5/00  UTC  the  agreement  appears 
satisfactory  and  at  6/00  UTC  it  appears  excellent.  However,  the 
period  centered  at  6/00  UTC  is  depicted  by  plus-signs,  and  that 
centered  at  6/03  UTC  by  triangles.  The  point  of  this  case  is 
that  the  rawinsonde  based  values  are  mean  values  over  a  very 
large  area  (  of  course,  the  radiosonde  data  contain  errors  and 
noise  which  we  have  ignored  here) ,  and  the  radar  time  mean  is  for 
only  one  position  inside  the  triangle,  and  that  even  though  the 
radar  is  near  the  center  of  the  triangle  the  time-mean  and  the 
space-mean  do  not  necessarily  match.  This  lack  of  agreement  does 
not  imply  any  error  in  either  the  radar  data  or  the  radiosonde 
data,  but  rather  points  to  the  frequent  mismatch  of  the  scales  of 
motion  they  represent.  In  the  next  case  we  will  use  NMC  data  to 
further  illustrate  this  notion  of  "shooting  at  a  moving  target". 

On  March  18th  an  occluded  front  approached  Flatland  from  the 
southwest,  as  illustrated  in  Figure  6.  The  mean  vertical  motions 
during  this  period  were  upward  as  seen  in  Figure  2  and  also  in 
Figure  7.  Green,  et  al.(1988),  used  proxy  indicators  of  vertical 
motions,  such  as  clouds  and  precipitation  from  surface  reports 
and  widespread  echoes  from  storm  surveillance^  radars,  to 
corroborate  the  upward  motions  observed  by  the  radar.  Also,  they 
noted  that  the  12 -hour  forecasts  from  the  NMC  model  showed  that  a 
region  of  large  vertical  motion  apparently  passing  Flatland 
during  this  period.  Computation  of  diagnostic  vertical  motions 
specific  to  the  Flatland  location  suitable  for  comparison  with 
time-mean  radar  observations  was  outside  the  scope  of  their 
study. 

We  have  computed  vertical  motions  over  Flatland  based  on  the 
NMC  gridded  analyses  using  the  kinematic  method,  the  adiabatic 
method,  and  the  adiabatic  quasi-geostrophic  omega-equation 
method.  Details  follow  Nastrom,  et  al.(1985).  The  results  for 
18/12  UTC  are  compared  with  9-hour  averages  from  the  radar  in 
Figure  8.  There  is  general  agreement  among  all  curves  that  the 
motion  is  upward  at  several  cm  s-1,  although  at  any  height  the 
comparison  is  not  exact.  We  have  included  radar  curves  for  two 


overlapping  time  periods  as  the  radar  is  representative  of  only 
a  portion  of  one  NMC  grid  square  and  the  storm  is  moving  in  this 
square  throughout  the  averaging  period.  The  large  difference 
between  the  radar  means  and  the  indirect  values  at  lowest  heights 
may  reflect  a  difference  in  the  scales  each  data  type  represents, 
or  it  may  be  due  to  limitations  of  the  indirect  methods,  such  as 
the  neglect  of  latent  heat  release  and  the  use  of  linear 
interpolation . 

3.  Conclusion 

The  Flatland  radar  experiment  has  provided  copious  data  on 
the  vertical  motions  over  flat,  mid-continental  temperate 
latitude  terrain,  free  from  the  effects  of  mountains. 

Preliminary  studies  of  this  data  set  have  begun  to  shed  light  on 
the  nature  of  the  motion  field  in  the  absence  of  nearby  sources 
of  orographic  effects.  In  particular,  qualitative  comparison  of 
the  time  series  of  vertical  motions  with  those  times  when  surface 
fronts  passed  or  were  near  Flatland  indicates  an  association  with 
active  periods  throughout  the  troposphere  and  sometimes  into  the 
stratosphere  of  large  lateral  extent.  Further,  the  passage  of 
fronts  is  often  accompanied  in  the  lower  troposphere  with  a 
region  of  upward  motion  ahead  of  and  downward  motion  behind  the 
front,  creating  a  signature  of  frontal  passage. 

Although  these  effects  related  to  synoptic-scale  features 
are  easily  observable,  the  observation  of  the  large-scale 
vertical  motion  itself  in  the  midst  of  larger  amplitude  signals 
due  to  small-scale  motions  is  difficult.  The  absence  of 
orographic  noise  in  the  signal  has  greatly  facilitated  the 
understanding  of  other  aspects  of  the  signal,  such  as  the  effect 
of  doppler-shifted  gravity  waves,  and  this  understanding  may 
contribute  to  the  eventual  understanding  of  the  large  scale 
signal  as  well. 

Verification  of  the  Taylor  hypothesis  through  comparison  of 
temporal  and  spatial  averages  is  also  difficult.  In  particular, 
it  has  become  clear  that  the  velocity  and  rate  of  decay  or 
development  of  a  given  system  will  determine  the  horizontal  scale 
to  which  a  temporal  average  at  the  radar  site  can  be  properly 
compared.  Only  in  certain  cases  will  this  scale  match  that  of 
the  radiosonde  network. 

Each  of  these  questions  has  important  implications  for  the 
future  application  of  MST  radars.  It  appears  that  a  very 
promising  avenue  for  future  progress  is  using  multiple  radars  to 
assess  spatial  correlation  features. 
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Figure  1.  Average  of  34  spectra  of  vertical  velocity  over  48 
hour  data  periods  at  Flatland  during  March  through  September, 
1987,  in  area  preserving  coordinates. 
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Figure  5.  Vertical  motions  computed  using  the  continuity 
equation  and  horizontal  winds  from  the  radiosonde  station 
triangle  around  Flatland  compared  with  6-hourly  mean  radar 
motions.  See  text. 


Figure  6.  Synoptic  weather  map  for  18  March  1987,  12  UTC. 
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Figure  8.  Vertical  motions  computed  from  NMC  gridded  analyses 
compared  with  9-hourly  mean  radar  motions.  See  text. 
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A  simple  model  for  the  enhanced  frequency  spectrum  of  vertical 
velocity  based  on  tilting  of  atmospheric  layers  by  lee  waves 
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This  paper  is  concerned  with  the  interpretation  of  the  enhanced 
frequency  spectrum  of  vertical  velocity  that  has  been  observed  at  many  loca¬ 
tions  under  strong  wind  conditions  in  the  vicinity  of  rough  terrain.  These 
enhanced  frequency  spectra  are  considerably  more  energetic  than  the  vertical 
velocity  spectra  observed  under  low-wind  conditions,  and  their  shape  resembles 
the  shape  of  the  oblique  spectra  observed  by  the  same  radars.  In  this  paper 
we  present  the  composite  spectra  of  horizontal  and  vertical  velocity  observed 
over  a  range  of  wind  conditions  by  the  Platteville,  Colorado  radar,  located  in 
the  lee  of  the  Colorado  Rockies.  The  vertical  velocity  spectra  show  clearly 
the  increase  in  magnitude  and  change  in  spectral  slope  that  occurs  as  the 
troposphere  becomes  increasingly  disturbed.  We  compare  the  observations  with 
a  simple  model  that  relates  the  magnitude  of  the  vertical  velocity  spectrum  to 
the  horizontal  velocity  spectrum  and  the  effective  tilt  of  atmospheric  layers. 

Our  results  show  that  it  is  possible  to  simulate  the  observed  vertical 
velocity  spectrum  with  effective  tilts  in  the  range  of  2.5  to  10  degrees. 

Recent  observations  of  vertical  velocity  spectra  from  the  Flatland  radar, 
located  in  flat  terrain  in  central  Illinois,  do  not  exhibit  this  behavior, 
which  leads  us  to  conclude  that  lee  waves  are  responsible  for  the  enhancement 
in  the  Colorado  spectra. 
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1.  INTRODUCTION 

Much  has  been  learned  in  recent  years  concerning  the  nature  of  the  ob¬ 
served  spectrum  of  mesoscale  variability  in  the  free  atmosphere.  Major 
progress  has  been  made,  for  example,  through  spectral  analysis  of  time  series 
of  Doppler  radar  observations  of  horizontal  and  vertical  velocities.  The 
frequency  spectra  determined  from  radar  observations  have  been  compared  with 
spectral  models  for  internal  waves  and  quasi- two-dimensional  turbulence 
[VanZandt,  1982,  1985;  Scheffler  and  Liu,  1985;  Gage  and  Nastrom,  1985a, b] . 

While  there  is  still  considerable  debate  concerning  the  nature  of  the 
spectrum  of  horizontal  velocities,  there  is  general  agreement  that  the  verti¬ 
cal  velocity  spectra  in  undisturbed  conditions  are  almost  entirely  due  to 
internal  waves.  This  view  is  supported  by  the  climatological  study  of  the 
frequency  spectra  of  vertical  motions  reported  by  Ecklund  et  al .  [1986]  and 
the  detailed  analysis  of  vertical  velocity  observations  from  the  Flatland 
radar  [Green  et  al.,  1988;  VanZandt  et  al.,  1989]. 

Although  the  vertical  velocity  spectrum  seems  to  result  from  internal 
waves,  the  horizontal  velocity  spectrum  may  be  expected  to  have  contributions 
from  both  internal  waves  and  quasi -two -dimensional  turbulence  as  discussed  in 
Gage  and  Nastrom  [1985a, b;  1986]  and  Gage  [1989],  Recently,  this  subject  has 
received  increased  attention  from  the  oceanographic  community  in  the  recent 
work  of  Muller  et  al.  [1988],  who  estimate  the  relative  contributions  of 
internal  waves  and  vortical  modes  to  the  ocean  current  spectra  observed  during 
IWEX.  The  coexistence  of  waves  and  quasi-two-dimensional  turbulence  is  an¬ 
ticipated  in  the  analysis  of  stratified  turbulence  contained  in  Riley  et  al. 
[1981]  and  Lilly  [1983],  The  subject  of  turbulence  in  stratified  fluids  has 
recently  been  reviewed  by  Hopfinger  [1987]. 

This  paper  is  concerned  with  the  explanation  of  the  enhanced  frequency 
spectra  of  vertical  motions  that  are  often  observed  near  mountains  under 


strong  wind  conditions.  We  hypothesize  that  these  enhanced  spectra  of  verti¬ 
cal  motions  are  due  to  a  contamination  of  the  vertical  motions  by  a  component 
of  the  quasi-horizontal  motions  that  occur  on  isentropic  surfaces.  According 
to  this  hypothesis,  when  isentropic  surfaces  are  tilted  as  they  are  by  lee 
waves,  the  vertical  beam  of  the  Doppler  radar  observes  a  component  of  the 
quasi-horizontal  motions.  This  hypothesis  is  tested  by  employing  a  simple 
model  that  relates  the  observed  vertical  velocity  spectrum  to  the  observed 
horizontal  velocity  spectrum.  Observations  taken  by  the  Platteville,  Colorado 
radar  are  used  for  this  purpose. 

2.  THE  ENHANCED  FREQUENCY  SPECTRUM  OF  VERTICAL  MOTIONS 

A  clear  illustration  of  the  enhanced  frequency  spectrum  of  vertical 
motions  is  shown  in  Figure  1.  These  spectra  were  reported  by  Ecklund  et  al. 
[1985]  and  were  taken  in  southern  France  during  the  ALPEX  experiment.  The 
enhanced  spectra  contrast  dramatically  with  the  relatively  flat  spectra  of 
vertical  velocity  observed  under  quiet  conditions  shown  in  Figures  1  and  2. 
Furthermore,  the  quiet- time  spectra  observed  during  ALPEX  closely  resemble 
vertical  velocity  spectra  observed  elsewhere  under  similar  conditions  [Ecklund 
et  al.,  1986].  They  also  resemble  the  Flatland  vertical  velocity  spectra 
reported  by  VanZandt  et  al.  [1989]  observed  under  most  conditions. 

The  possibility  that  the  enhanced  frequency  spectrum  of  vertical 
motions  might  be  due  to  Doppler  shifting  of  an  internal  wave  spectrum  by 
strong  winds  needs  to  be  considered.  Both  Scheffler  and  Liu  [1986]  and  Fritts 
and  VanZandt  [1987]  have  analyzed  the  effect  of  mean  winds  on  a  spectrum  of 
internal  waves.  Recently,  VanZandt  et  al.  [1989]  have  shown  that  the  rela¬ 
tively  small  changes  of  vertical  velocity  spectra  that  are  observed  at 
Flatland  can  be  explained  by  Doppler  shifting. 

The  fact  that  enhanced  vertical  velocity  spectra  are  not  observed  at 
Flatland  is  a  strong  argument  that  the  enhanced  vertical  velocity  spectra 
observed  at  other  locations  is  due  to  the  influence  of  terrain  and  not  to 


Doppler  shifting.  In  the  remainder  of  this  paper,  we  examine  the  enhanced 
vertical  velocity  spectra  that  are  observed  at  Platteville,  Colorado  in  the 
lee  of  the  Rocky  Mountains. 

3.  FREQUENCY  SPECTRA  OF  VERTICAL  AND  HORIZONTAL  MOTIONS 
OBSERVED  AT  PLATTEVILLE,  COLORADO 

The  Platteville  radar  has  been  in  continuous  operation  since  the  early 
1980's.  It  was  originally  constructed  as  a  prototype  for  the  Poker  Flat  MST 
radar  [Ecklund  et  al.,  1979]  and  has  been  operated  in  recent  years  by  NOAA's 
Wave  Propagation  Laboratory.  The  spectra  presented  in  Figure  3  and  Figure  4 
are  composite  spectra  averaged  by  season  and  stratified  by  the  standard  devia¬ 
tion  of  vertical  velocity.  At  Platteville,  there  is  a  close  relationship 
between  the  background  wind  speed  and  the  vertical  velocity  variance  as  evi¬ 
dent  in  Figure  3  and  Figure  4. 

Composite  vertical  velocity  spectra  for  the  winter  season  at 
Platteville,  Colorado  are  shown  in  Figure  3.  These  spectra  show  clearly  the 
enhancement  in  spectral  magnitude  and  systematic  variation  of  spectral  slope 
that  accompanies  increasing  winds  at  this  location.  The  dashed  curve  indi¬ 
cates  a  spectral  slope  of  -5/3.  With  increasing  wind  speed  the  observed 
vertical  velocity  spectra  approach  the  -5/3  spectral  slope.  Only  those  time 
periods  when  both  horizontal  and  vertical  wind  components  were  measured  are 
used  here,  although  the  results  are  very  similar  when  all  available  data  are 
used. 

Corresponding  composite  zonal  wind  spectra  at  5.8  km  for  the  winter 
season  as  stratified  by  vertical  velocity  variance  are  shown  in  Figure  4. 

These  spectra  also  show  a  systematic  but  less  pronounced  variation  with  back¬ 
ground  wind  speed.  Note  that  the  magnitude  of  the  horizontal  velocity  spectra 
are  in  all  cases  considerably  larger  than  the  magnitude  of  the  corresponding 
vertical  velocity  spectra.  Again,  the  dashed  line  shows  the  -5/3  slope  for 
comparison  with  the  observed  spectra. 
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4.  A  SIMPLE  MODEL  FOR  THE  ENHANCED  VERTICAL  VELOCITY  SPECTRUM 


The  simple  model  used  here  to  explain  the  enhanced  vertical  velocity 
spectrum  is  based  on  the  idea  that  in  stably  stratified  flows  the  velocity 
field  is  comprised  of  both  internal  waves  and  potential  vorticity  modes  [Riley 
et  al. ,  1981;  Lilly,  1983;  Muller,  1984;  Herring  and  Metais,  1989],  The 
potential  vorticity  modes  are  comprised  of  quasi-horizontal  eddies  that  follow 
very  closely  isentropic  surfaces  in  the  stratified  fluid.  Accordingly,  when 
the  isentropic  surfaces  are  horizontal  and  undisturbed,  vertical  motion  will 
be  due  only  to  internal  waves.  The  horizontal  motion  field,  however,  will  be 
comprised  of  both  internal  waves  and  vortical  modes.  When  isentropic  surfaces 
are  tilted,  the  vortical  modes  will  also  be  tilted  to  conserve  potential 
vorticity  so  that  there  will  now  be  a  vertical  component  to  the  vortical 
motion . 

The  consequences  of  the  model  on  observations  by  a  vertically  directed 
radar  are  simple  and  straightforward.  When  isentropic  surfaces  are  flat,  the 
radar  will  observe  only  the  internal  gravity  wave  field  and  vortical  motions 
will  not  be  observed.  However,  when  isentropic  surfaces  are  tilted,  as  they 
may  be  by  mountain  lee  waves,  the  radar  will  observe  a  component  of  the  vorti¬ 
cal  motion. 

The  process  can  be  quantified  as  follows.  If  the  effective  tilting 

angle  is  S,  the  observed  vertical  velocity  spectrum  $  (w)  is  related  to  the 
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Thus  if  6  were  known  it  would  be  possible  to  estimate  the  vertical  velocity 
spectrum  from  the  observed  horizontal  velocity  spectrum  since  the  internal 
wave  spectrum  is  fairly  well  known.  From  the  Platteville  observations,  we  can 
deduce  the  effective  tilting  angle  that  satisfies  (1). 
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5.  VERTICAL  VELOCITY  SPECTRA  DEDUCED  FROM  THE  SIMPLE  MODEL 

Vertical  velocity  spectra  can  be  deduced  from  the  observed  horizontal 
velocity  spectrum  by  first  determining  the  effective  tilting  angle.  The 
effective  tilting  angle  is  determined  by  the  value  of  5  required  to  reduce  the 
horizontal  velocity  spectrum  to  the  observed  magnitude  of  the  vertical 
velocity  spectrum.  For  example,  Figure  5  shows  that  S  -  10°  brings  the  ob¬ 
served  horizontal  velocity  spectrum  for  U  -  17  ms  1  into  near  coincidence  with 
the  observed  vertical  velocity  spectrum.  Note  that  only  when  the  effective 
tilt  angle  is  less  than  about  5°  does  the  internal  wave  spectrum  become  impor¬ 
tant  in  (1)  for  these  observations. 

Model  vertical  velocity  spectra  for  the  various  background  horizontal 
velocities  in  Figure  4  are  shown  in  Figure  6.  Comparison  of  the  two  sets  of 
spectra  shows  excellent  agreement.  The  effective  tilt  angles  that  are  consis¬ 
tent  with  the  model  fall  in  the  range  of  2.5°  -  10°.  This  magnitude  of  tilt 
is  easily  produced  by  lee  waves  that  possess  horizontal  wavelengths  of  order 
ten  kilometers  and  vertical  displacements  ranging  from  several  hundred  meters 
to  greater  than  one  kilometer  [Gage,  1986].  Of  course,  in  the  absence  of 
mountain  lee  waves  tilting  would  not  be  sufficient  to  produce  the  enhanced 
vertical  velocity  spectrum.  On  the  synoptic  scale  isentropic  surfaces  are 
generally  tilted  less  than  a  few  tenths  of  a  degree.  Propagating  internal 
waves  can  cause  tilting  of  perhaps  a  few  degrees  which  could  cause  a  slight 
enhancement  in  the  observed  magnitude  of  vertical  velocity  spectra  at  very  low 
frequencies  even  over  flat  terrain. 


6.  CONCLUSIONS 


A  simple  model  has  been  used  Co  simulate  the  enhanced  frequency  spectra 
of  vertical  motion  observed  at  Platteville,  Colorado  under  disturbed  condi¬ 
tions.  If  it  is  hypothesized  that  the  horizontal  velocity  spectrum  is 
primarily  due  to  quasi-horizontal  motions  associated  with  potential  vorticity 
modes ,  the  magnitude  and  shape  of  the  observed  vertical  velocity  spectrum  can 
be  explained  by  tilting  of  isentropic  surfaces  due  to  lee  waves. 

While  the  simple  model  proposed  here  is  reasonably  consist 
observations  at  Platteville,  more  stringent  tests  are  required  before  too  much 
confidence  can  be  placed  in  the  result.  For  example,  independent  estimates  of 
the  slope  of  isentropic  surfaces  associated  with  lee  waves  expected  under 
various  background  wind  conditions  would  strengthen  the  analysis. 
Alternatively,  observations  ta cen  simultaneously  at  more  than  one  oblique 
zenith  angle  would  help  with  the  evaluation  of  how  much  of  the  horizontal 
velocity  spectra  is  associated  with  vortical  modes  and  how  much  is  due  to 


internal  waves. 
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FIGURE  CAPTIONS 


Figure  1.  The  frequency  spectrum  of  vertical  velocity  observed  in  southern  France 
during  ALPEX  for  active  and  quiet  days.  (After  Ecklund  et  al . , 

[1985]). 


Figure  2. 

Figure  3. 


Figure  4. 


Figure  5a. 

Figure  5b. 
Figure  5c. 
Figure  5d. 
Figure  6 . 


The  frequency  spectrum  of  vertical  velocity  observed  in  southern  France 
during  ALPEX  for  quiet  days.  (After  Ecklund  et  al .  ,  [1985]). 

The  frequency  spectrum  of  vertical  velocity  observed  at  Plattevillle , 

CO  during  the  winter  season  at  5.8  km.  Spectra  shown  are  composite 
spectra  averaged  over  several  years  and  stratified  by  vertical  velocity 
variance . 

The  frequency  spectra  of  zonal  wind  speed  observed  at  Platteville,  CO 
during  the  winter  season  at  5.8  km.  The  spectra  represent  the  same 
period  of  observation  as  the  spectra  in  Figure  3  and  are  stratified  in 
a  similar  manner. 

An  illustration  of  the  deduction  of  the  vertical  velocity  spectrum  from 
the  observed  horizontal  velocity  spectrum  using  (1).  Also  shown  is  the 
observed  vertical  velocity  spectrum  for  u  -  17  ms  ^  for  comparison. 
Same  as  for  Figure  5a  except  for  u  -  14  ms 

Same  as  for  Figure  5a  except  for  u  -  9  ms  . 

Same  as  for  Figure  5a  except  for  u  -  4  ms  ^ . 

Frequency  spectra  of  vertical  velocity  derived  from  the  observed 
horizontal  velocity  according  to  (1) .  These  spectra  should  be  compared 
with  the  observed  spectra  in  Figure  3. 
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ABSTRACT 


Whenever  the  time  series  of  a  variable  is  not  sampled 
continuously,  a  fraction  of  its  variance  is  missed.  This 
fraction  contributes  to  the  uncertainty  of  the  mean  value 
determined  from  the  available  sample.  In  the  case  of 
vertical  velocity,  where  a  relatively  large  fraction  of  the 
variance  is  due  to  high  frequency  motions,  the  effect  of 
even  small  gaps  in  the  sampling  is  very  important. 
Measurements  from  the  Flatland  VHF  radar  are  used  to  study 
this  effect.  Quasi-continuous  data  are  sampled  with  several 
strategies,  and  it  is  found  that  the  strategy  with  the 
smallest  gap  between  data  samples  gives  the  smallest 
uncertainty  of  the  mean  vertical  velocity.  It  is  found  that 
the  time  between  statistically  independent  samples  of 
vertical  velocity  at  Flatland  is  5.1  minutes  in  the 
troposphere. 
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1.  INTRODUCTION 


Clear-air  Doppler  radars  can  provide  profiles  of  the 
vertical  velocity  over  a  station  on  a  nearly  continuous 
basis,  and  at  many  research  sites  these  data  are  recorded  at 
intervals  on  the  order  of  a  minute.  In  some  cases  it  is 
necessary  to  use  data  samples  spaced  at  greater  intervals 
than  at  the  reporting  frequency  and  in  recent  years  there 
have  been  several  studies  which  have  used  time-means  of  the 
vertical  velocity.  For  example,  Nastrom,  et  al.(1985), 
compared  time-means  of  vertical  velocity  from  MST  radar 
measurements  with  large-scale  values  based  on  radiosonde 
data;  and  Larsen,  et  al. (1988) ,  compared  them  with  ECMWF 
model  forecasts  of  vertical  velocity.  Cornish  (1988) 
mentioned  the  difficulty  of  obtaining  long-term  mean  values 
at  the  Arecibo  observatory,  and  Balsley,  et  al.  (1988) ,  in  a 
much  more  extensive  study,  used  time  mean  vertical 
velocities  to  study  tropical  motions  at  7°N.  Strauch,  et 
al.  (1987),  studied  the  importance  of  including  the  effects 
of  large  vertical  motions  in  the  measurements  of  the 
horizontal  wind  components  from  obliquely  directed  radar 
beams.  For  purposes  such  as  these  it  is  important  to  use 
the  sampling  strategy  which  yields  the  mean  value  of 
vertical  velocity  with  the  smallest  possible  uncertainty 
under  given  experimental  constraints  and  to  estimate  the 
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magnitude  of  the  uncertainty  of  the  mean  once  a  sampling 
strategy  is  selected. 

Because  the  spectrum  of  atmospheric  vertical 
velocity  covers  a  complete  range  of  frequencies,  some 
variance  will  be  missed  whenever  the  temporal  sampling  is 
not  continuous.  This  unresolved,  or  missed,  variance 
contributes  to  the  uncertainty  of  mean  values  based  on  the 
available  observations,  along  with  other  effects  such  as 
aliasing  of  the  observed  frequency  spectrum.  In  this  paper 
we  examine  the  effects  of  various  temporal  sampling 
strategies  on  estimates  of  the  time-mean  vertical  velocity 
from  MST  radar  measurements. 

2 .  DATA 

Data  used  for  this  study  are  vertical  velocities  from 
the  Flatland  VHF  radar.  The  radar  operating  parameters  and 
examples  of  the  data  are  given  by  Green,  et  al.  (1988) . 
Briefly,  during  the  period  of  record  used  here,  March 
through  May,  1987,  the  radar  was  operated  with  one  beam 
which  was  pointed  vertically.  Doppler  spectra  were  averaged 
for  about  2.5  minutes  and  then  recorded  on  tape.  Mean 
velocities  at  2.5-minute  intervals  were  computed  from  these 
average  spectra.  The  available  data  are  thus  quasi- 
continuous  in  time  as  only  a  few  seconds  gap  was  missed  each 
data  processing  and  recording  cycle.  In  practice,  some 
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observations  are  missing  due  to  insufficient  signal-to-noise 
ratio  of  the  doppler  spectrum,  or  other  problems  such  as 
contamination  by  reflections  from  aircraft.  Also,  there  are 
occasional  short  gaps  due  to  equipment  malfunction.  For  the 
analyses  of  frequency  spectra  and  the  autocorrelation 
function  given  here  a  segment  of  data  was  not  used  at  a 
given  height  if  a  gap  longer  than  10  minutes  occurred. 


3.  STATISTICAL  FEATURES  OF  THE  VERTICAL  VELOCITY 

The  fluctuations  of  vertical  velocity  at  Flatland  have 
the  most  energy  at  high  frequencies,  as  can  be  seen  in  the 
frequency  spectrum  (Figure  1)  which  shows  that  most  of  the 
energy  falls  between  about  10  and  30  minutes  period.  The 
coordinates  in  Figure  1  are  area  preserving,  i.e.,  the 
variance  in  a  given  frequency  band  is  proportional  to  the 
area  under  the  curve.  From  comparisons  with  a  model  of  a 
spectrum  of  gravity  waves  subjected  to  doppler  shifting  by 
horizontal  winds.  Van  Zandt,  et  al.  (1989),  conclude  that 
the  observed  spectrum  of  vertical  velocity  at  Flatland  is 
due  mostly  to  gravity  waves.  This  is  in  contrast  to  the 
spectra  of  vertical  velocity  at  locations  in  or  near 
mountains  which  show  the  effects  of  standing  lee  wave 
activity  (Ecklund,  et  al.,  1986;  Gage  and  Nastrom,  1989). 

Because  a  large  portion  of  the  variance  at  Flatland  is 
due  to  relatively  high  frequency  variations  we  should  expect 
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that  a  relatively  high  sampling  rate  will  be  required  to 
avoid  large  uncertainties  in  mean  values  computed  from  the 
observations.  In  Figure  1  we  note  that  the  peak  variance  is 
found  at  about  12  minutes  period,  and  in  order  to  avoid 
serious  sampling  problems  we  might  expect  that  observations 
should  be  taken  more  often  than  about  half  this  period, 
i.e.,  every  6  minutes  or  so.  Another,  more  rigorous, 
estimate  of  the  time  between  independent  samples  can  be 
obtained  from  analysis  of  the  autocorrelation  function  as 
discussed  next. 

Figure  2  shows  the  autocorrelation  function,  R(t) ,  at 
Flatland  at  5  km  for  lags  less  than  about  15  minutes  and, 
for  comparison,  the  R(t)  functions  at  Poker  Flat,  Alaska, 
from  Nastrom  and  Gage  (1983) .  R(t)  can  be  used  to  estimate 
the  effective  time  between  independent  samples  (T)  using  the 
relation  (Leith,  1973) 


where  t  is  the  lag.  An  important  application  of  T  is  in 
computing  statistical  significance  levels  (Mitchell,  et  al., 
1966)  from  time  series  of  observed  data.  Nastrom  and  Gage 
(1983)  found  that  at  Poker  Flat  T  ranged  from  about  15 
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minutes  for  the  curves  labeled  "winter  1"  and  "summer", 
which  represent  active  periods,  to  about  9  minutes  for  the 
curve  labeled  "winter  2",  which  represents  quiet  periods. 
There  were  not  enough  cases  to  prepare  a  summer  curve  for 
quiet  periods,  and  only  the  curve  for  active  times  is  shown. 
The  curve  labeled  "Flatland"  was  prepared  by  averaging  the 
autocorrelation  functions  from  each  6-hour  period  available. 
The  data  within  each  6-hour  period  were  linearly 
interpolated  to  153-second  intervals  and  residuals  from  a 
linear  trend  line  fit  by  least  squares  were  then  used  to 
compute  R(t) .  Note  that  the  average  curve  at  Flatland 
decreases  more  rapidly  than  all  of  the  curves  from  Poker 
Flat,  suggesting  that  the  vertical  velocity  at  Flatland  has 
a  relatively  short  "memory".  The  semi-logarithmic  plots  in 
Figure  2  show  that  an  appropriate  analytic  model  after  about 
lag  5-minutes  is  R(t)  =  c  exp(-vt),  where  v  defines  the 
rate  of  decrease.  At  shorter  lags  the  observed  curves  fall 
above  this  model,  suggesting  that  the  observed 
autocorrelation  function  is  a  mixed  first-order  moving 
average,  first-order  autoregressive  process  as  discussed  in 
more  detail  by  Nastrom  and  Gage  (1983)  .  Using  (1)  with  this 
model  of  R(t)  we  estimate  that  T  at  Flatland  is  5.1  minutes, 
which  is  less  than  the  estimate  for  any  of  the  curves  at 
Poker  Flat.  Also,  we  note  that  this  estimate  (5.1  minutes) 
is  near  that  of  6  minutes  inferred  above  from  Figure  1. 

This  result  does  not  mean  that  observations  should  be  taken 
only  each  5.1  minutes  at  Flatland,  but  rather  that,  on  the 
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average,  independent  information  is  gained  only  at  this 
rate.  In  fact,  considerable  uncertainty  of  the  mean  is 
found  when  only  every-other  observation  is  used  (i.e.,  with 
data  spaced  about  5.1  minutes),  as  discussed  next. 

4.  EFFECTS  OF  SAMPLING  STRATEGY  ON  MEAN  VERTICAL  VELOCITY 


The  analysis  approach  we  will  use  will  be  to  form 
averages  over  one-  , two-,  three-,  and  six-hour  periods  using 
the  complete  set  of  observations.  These  averages  will  be 
used  as  the  standards  for  comparison.  Next,  the  data  set 
will  be  degraded  to  simulate  various  sampling  strategies  and 
the  mean  values  thus  obtained  will  be  compared  with  the 
"true"  mean  values  based  on  complete  data.  This  procedure 
will  give  estimates  of  the  uncertainty  of  measurement  due  to 
temporal  sampling  strategy.  There  are,  of  course,  other 
sources  of  uncertainty,  such  as  system  limitations  and  the 
effects  of  spatial  averaging,  which  might  be  important  in 
practice  and  which  are  beyond  the  scope  of  this  study. 

Four  schemes  were  used  to  illustrate  the  effects  of  a 
strategy  which  uses  incomplete  sampling.  Three  of  them  use 
one-half  of  the  data  but  in  varying  sized  segments:  the 
first  uses  alternate  observations,  the  second  uses  data  for 
every-other  10  minutes,  and  the  third  uses  data  for  every- 
other  30  minutes.  Each  of  these  could  simulate  a  strategy 
applied  in  a  field  experiment  where,  for  example,  the  radar 
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was  used  to  gather  vertical  data  half  of  the  time  and 
oblique  data  half  of  the  time.  The  fourth  strategy  uses 
data  for  the  first  12  minutes  of  each  hour,  similar  to  the 
observation  cycle  of  the  data  used  by  Larsen,  et  al.  (1988). 
The  results  from  these  strategies  will  be  labeled  "alt  obs", 
"0-10  min",  "0-30  min",  and  "12  min"  in  Figures  3  and  4  and 
in  Table  1. 

RMS  differences  between  mean  values  based  on  these  four 
sampling  strategies  and  the  "true"  mean  computed  from  all 
available  data  are  given  in  Figure  3  as  functions  of 
averaging  period.  Two  patterns  are  evident  in  Figure  3: 
first,  the  RMS  differences  decrease  with  increasing 
averaging  time.  This  decrease  is  consistent  with  the  notion 
that  the  degraded  observation  cycles  are  sufficient  to 
resolve  long-period  variations  reflected  in  the  left-hand 
tail  of  the  spectrum  in  Figure  1,  and  that  the  variance  due 
to  longer  periods  thus  contributes  to  the  uncertainty  of  the 
mean  only  for  shorter  averaging  periods. 

The  second  pattern  seen  in  Figure  3  is  that  larger  gaps 
between  data  samples  lead  to  larger  uncertainty  of  the 
estimates  of  the  mean.  This  effect  reflects  the  increased 
spectral  leakage  caused  by  longer  gaps,  an  effect  discussed 
by  Baer  and  Tribbia  (1976).  They  found  that  there  is 
reduced  spectral  fidelity  for  all  frequencies  higher  than 
the  frequency  corresponding  to  the  longest  gap  in  a  data 


sequence.  The  reduced  spectral  fidelity  in  the  present 
context  leads  to  greater  uncertainty  of  the  estimate  of  the 
mean  value. 

Finally,  we  note  that  the  uncertainty  of  the  estimate 
of  the  mean  is  relatively  very  large  when  only  12  minutes  of 
data  per  hour  are  used.  For  a  one-hour  averaging  period  the 
RMS  difference  from  the  true  mean  is  nearly  9  cm  s-1  for 
this  strategy,  which  is  over  half  as  large  as  the  standard 
deviation  of  all  data  (about  14  cm  s-1) . 

Other  statistical  quantities  can  be  used  to  illustrate 
the  patterns  seen  in  Figure  3.  For  example,  Figure  4  shows 
the  mean  deviation,  defined  as  the  average  absolute  value  of 
the  difference  between  the  true  mean  and  that  from  the 
various  sampling  strategies,  in  a  format  similar  to  Figure 
3 .  The  patterns  are  consistent  with  those  in  Figure  3 , 
although  the  numerical  values  of  this  measure  of  uncertainty 
are  only  about  two-thirds  as  large  as  the  RMS  values. 

Another  measure  of  the  capability  of  one  estimate  of  the 
mean  to  track  the  true  mean  is  the  linear  correlation 
coefficient,  r.  Table  1  gives  the  values  of  r  over  N  pairs 
of  means  with  the  true  mean  for  the  same  averaging  times  and 
sampling  strategies  used  in  Figures  3  and  4.  The  percentage 
of  variance  in  one  variable  explained  by  another  is  given  by 
100  r2 .  Applying  the  results  in  Table  1  shows  that  for  all 
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averaging  periods  the  strategy  of  using  alternate  2. 5-minute 
data  intervals  accounts  for  90  percent  or  more  of  the 
variance  of  the  true  mean  in  all  cases  here,  while  using 
every-other  10-minute  period  accounts  for  over  80  percent  of 
the  variance  for  averaging  periods  of  2  hours  or  more.  At 
the  extreme,  using  only  the  first  12  minutes  of  each  hour 
the  estimates  of  the  means  account  for  less  than  50  percent 
of  the  variance  of  the  true  mean  except  for  the  6-hour 
averaging  period.  These  results  show  that  the  uncertainty 
of  estimates  of  the  mean  vertical  velocity  is  minimized  when 
the  average  length  of  gaps  in  the  data  is  minimized. 


5.  CONCLUSIONS. 

The  following  conclusions  are  based  on  analysis  of 
vertical  velocities  measured  by  the  Flatland  MST  radar 
during  the  spring  of  1987: 

1.  The  autocorrelation  function  of  vertical  velocity 
resembles  that  expected  for  a  mixed  first-order  moving- 
average,  first-order  autoregressive  process. 

2.  The  average  time  between  independent  samples  in  the 
data  from  Flatland  used  here  is  about  5  minutes. 
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3.  The  uncertainty  of  the  mean  value  of  vertical 
velocity  decreases  as  the  averaging  period  used  to  compute 
the  mean  increases  regardless  of  the  sampling  strategy  used. 

4.  The  uncertainty  of  the  mean  value  of  vertical 
velocity  increases  as  the  length  of  gaps  in  the  data 
increases.  Because  of  this  effect,  every  effort  should  be 
made  to  minimize  gaps  in  samples  of  vertical  velocity  in 
future  field  programs. 

These  results  are  concerned  only  with  temporal  sampling, 
and  in  practice  other  effects  such  as  instrument  errors  may 
require  attention.  In  particular,  these  results  do  not 
include  any  sampling  uncertainties  which  might  arise  from 
spatial  variability.  It  remains  as  a  future  experiment  to 
compare  averages  from  nearby  stations  as  functions  of 
averaging  time,  sampling  strategy,  and  other  variables  to 
study  such  effects. 
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AUTOCORRELATION 


Figure  2.  Autocorrelation  function  of  vertical  velocity  at 
Flatland  and  at  Poker  Flat.  See  text. 
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INTRODUCTION 

An  experimental  campaign  was  conducted  in  June  1988  at  the  Flatland  VHF 
clear-air  Doppler  Radar  site,  which  is  located  in  very  flat  terrain  far  removed 
from  mountains  near  Champaign-Urbana,  Illinois,  to  measure  height  profiles  of 
tha  refractivity  turbulence  structure  parameter  C2  and  related  turbulent 
parameters.  Three  different  techniques  were  used  to  measure  profiles  of  C2: 
the  Flatland  radar  (GREEN  et  al. ,  1988),  a  stellar  scintillometer  (OCHS  etnal., 
1977),  and  thermosonde  balloon  flights  (BROWN  et  al.,  1982).  Both  the  radar 
and  scintillometer  measured  C2  remotely  by  sensing  the  effect  of  the  fluctua¬ 
tions  of  refractive  index  on  ?he  propagation  of  electromagnetic  waves.  At  radio 
wavelengths  the  humidity  and  its  gradient  as  well  as  the  temperature  gradients 
contribute  significantly  to  the  refractive  index,  whereas  at  optical 
wavelengths  only  the  temperature  fluctuations  are  important.  In  contrast  to 
the  remote  sensors,  the  thermosonde  measures  the  temperature  fluctuations 
directly  by  measuring  the  RMS  temperature  difference  between  two  very  fast 
temperature  sensors  separated  by  a  meter.  The  balloon  instrument  also  made 
measurements  of  the  height  profile  of  the  standard  thermodynamic  parameters 
(pressure,  temperature,  and  humidity),  plus  wind  speed  and  direction.  Model 
estimates  of  C2  were  calculated  from  the  thermodynamic  and  wind  measurements 
using  the  numerical  methods  described  by  WARNOCK  AND  VANZANDT  (198 3).  In 
addition,  optical  measurements  were  made  of  the  transverse  coherence  length 
(EATON  et  al . ,  1988),  and  of  the  isoplanatic  angle  (EATON  et  al.,  1985).  Both 
these  parameters  depend  on  a  weighted  integrated  value  of  C2  through  the  atmos¬ 
phere. 


Comparisons  among  these  measurements  taken  simultaneously  in  simple  topog¬ 
raphy  provide  a  unique  opportunity  to  compare  these  different  measurement 
techniques  and  to  contrast  these  measurements  with  previous  comparisons  made  in 
rough  terrain  (e.g.,  GOOD  et  al.,  1982;  GREEN  et  al. ,  1984;  EATON  et  al., 

1988).  In  this  paper  we  present  some  preliminary  results,  which  emphasize  the 
radar  measurements. 


EXPERIMENTAL  SET  UP 


Thermosonde  System 


The  thermosonde  system  consisted  of  two  parts,  which  were  mounted  below  an 
ascending  balloon  (BROWN  et  al .  ,  1982).  One  part  was  a  standard  VIZ  digital 
microsonde,  the  other  was  a  micro- thermo  bridge  thermosonde  linked  to  the 
microsonde.  The  thermosonde  measured  the  RMS  temperature  fluctuations  between 
two  unheated  fine  wire  tungsten  probes  separated  horizontally  by  one  meter. 

The  noise  level  of  the  instrument  is  about  0.002°  C,  and  the  data  were  recorded 
every  four  seconds  giving  about  a  20-meter  height  resolution.  The  microsonde 
makes  excellent  height-resolution  measurements  of  the  standard  thermodynamic 
parameters  (pressure,  temperature,  and  relative  humidity),  and  wind  speed  and 
direction.  The  thermodynamic  data  were  recorded  every  four  seconds,  which  gave 
about  a  20  m  height  resolution.  The  wind  speed  and  direction  were  determined 
by  using  the  Loran-C  navigator  system,  and  the  data  were  recorded  every  ten 
seconds,  which  gave  about  a  50  m  height  resolution.  A  1200-gram  meteorological 
balloon  was  used  to  lift  the  instrument  package  with  an  ascent  rate  of  about 
5  m/s.  The  package  was  suspended  from  90  to  180  meters  below  the  balloon  to 
ensure  that  the  turbulent  wake  from  the  two  meter  balloon  did  not  affect  the 
measurements . 

Flatland  Radar 


The  Flatland  radar  (GREEN  et  al.,  1988)  is  located  about  8  km  west  of  the 

Champaign-Urbana  Airport  [40.05°  N,  88.38°  E,  212  m  above  mean  sea  level 

(MSL) ] .  This  clear-air  Doppler  radar  (also  called  wind  profiler  or  ST  radar) 

operates  at  a  frequency  of  49.8  MHz  (6.02  m  wavelength),  with  peak  and  average 

power  of  about  10  kw  and  150  watts,  respectively.  The  pulse  length  and  the 

range  resolution  are  selectable  from  150  to  2400  m.  All  the  data  reported  in 

this  paper  were  taken  with  a  pulse  length  of  1.5  km  and  over-sampled  with  a 

range  resolution  of  750  m.  In  this  experiment,  the  antenna  beam  was  tilted 

twenty  degrees  off  the  vertical  toward  the  east  or  south  directions.  With  this 

tilt  angle,  echoes  due  to  specular  scattering  were  essentially  eliminated,  so 

that  unambiguous  C2  measurements  were  obtained, 
n 

Optical  Systems 

Three  optical  systems  were  used  in  this  study:  a  stellar  scintillometer, 
an  isoplanometer  (EATON  et  al . ,  1985),  and  a  transverse  coherence  length  system 
(EATON  et  al.,  1988).  The  isoplanometer  measures  the  isoplanatic  angle  8  , 
which  is  the  maximum  angular  extent  of  an  extended  object  that  can  be  viewed 
through  turbulence,  and  the  transverse  coherence  system  measures  the  transverse 
coherence  length  r  ,  which  is  related  to  the  spread  of  a  star  image  due  to 
turbulence.  The  scintillometer,  described  below,  measures  a  height  profile  of 
C2 ,  whereas  both  the  isoplanometer  and  transverse  coherence  systems  measure  an 
optical  quantity  which  is  related  to  the  weighted  integrated  value  of  C2 
through  the  atmosphere,  i.e.,  8q  (radians)  and  r^  (meters)  are 
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where  k  is  the  wavenumber  of  the  light  and  z  is  the  height  above  ground. 


S tellar  Sc ir.t il lome ter  Model  II 


The  stellar  scintillometer  Model  II  was  developed  by  the  Wave  Propagation 
Laboratory  of  NOAA  (OCHS  et  al . ,  1977).  To  operate  the  system,  a  star  of 
second  magnitude  or  brighter  and  within  45°  of  the  zenith  is  selected.  The 
system  is  sensitive  to  spatial  wavelengths  ranging  from  5  to  15  cm,  and 
measurer  seven  different  height  regions  of  optical  turbulence  ranging  from  2.2 
to  18.5  km  above  ground  level  (AGL) .  The  height  weighting  functions  for  these 
seven  heights  are  broad,  and  are  broadest  at  the  highest  altitudes. 

Surface  Measurements 


An  instrumented  tower  was  installed  near  the  telescope  domes.  In  addition 
to  the  standard  meteorological  measurements,  the  temperature  structure 
parameters  C2  and  the  solar  radiation  were  measured.  The  temperature,  dew 
point  temperature,  wind  speed  and  direction,  and  C2  measurements  were  made  at 
both  one  and  four  meters  above  the  surface,  and  the  pressure  and  solar  radia¬ 
tion  were  each  measured  at  one  height,  one  and  four  meters,  respectively. 

Model  Estimates  of  C2 

- n 

Model  calculations  of  C2  were  made  from  the  microsonde  upper  air  data; 
i.e.,  pressure,  temperature , humidity,  and  wind  speed  and  direction.  The  basic 
model  concepts  are  given  by  VANZANDT  et  al .  (1981).  Since  then,  the  model  has 
been  extended  and  revised  considerably.  We  used  the  latest  version  described 
by  WAR-NOCK  et  al .  (1985),  and  used  the  numerical  techniques  described  by 
WARNOCK  AND  VANZANDT  (1985),  to  evaluate  the  model  estimates.  To  compare  the 
model  estimates  with  results  from  previous  comparisons  with  the  Sunset  and 
Stapleton  radars,  which  are  located  in  and  near  the  mountains,  respectively 
(WARNOCK  et  al .  1986;  1988),  we  used  identical  values  of  all  the  model 
parameters  and  constants  except  one:  the  constant  in  the  equation  giving  the 
distribution  of  wind  shear.  This  parameter  quantifies  the  shear  environment  in 
the  range  of  scales  important  in  the  onset  of  turbulence  flow;  therefore,  the 
evaluation  of  this  parameter  allows  the  shear  environment  at  these  scales  above 
the  Flatland  radar  to  be  contrasted  with  its  value  above  rough  terrain. 

Description  of  Campaign 

This  experimental  campaign  was  conducted  at  the  Flatland  radar  site  from  6 
to  15  June  1988.  The  scintillometer  and  isoplanometer  systems  shared  a  single 
telescope;  the  scintillometer  operated  at  night  and  the  isoplanometer  operated 
during  the  day  and  in  twilight.  Thermosonde  balloon  launches  were  usually  made 
after  noontime  and  before  midnight  local  time.  A  typical  schedule  was  to 
launch  a  few  packages  in  the  afternoon  and  a  few  after  dark.  The  radar 
operated  during  clear  and  cloudy  conditions,  whereas  the  optical  systems 
operated  during  clear  sky  conditions.  The  optical  telescopes  were  located 
about  forty  meters  northeast  of  the  center  of  the  radar  antenna,  and  the 
thermosonde  launch  site  was  about  sixty  meters  east  of  the  antenna  center. 

PRELIMINARY  RESULTS 

Measurements  made  by  both  the  radar  and  scintillometer  remote  sensors  are 
average  values  over  both  time  and  space.  Five -minute  data  averages  were  used 
in  this  paper  by  both  systems,  and  the  height  profiles  measured  by  both  systems 
are  relatively  smooth.  In  contrast,  each  thermosonde  in  situ  measurement  is  an 
average  over  four  seconds,  which  gives  about  a  20-meter  height  resolution.  The 
four-second  thermosonde  data  displays  many  very  tA/f/ 'large  peaks  in  C2.  ■> - 

Frequently,  the  measured  C2  values  fall  to  the  observing  noise  level.  ^Thus ,  to 
compare  the  thermosor.de  profile  with  the  others  we  smoothed  the  thermosonde 
data.  We  used  a  Gaussian  filter  with  a  -  0.5  km  and  truncated  the  filter  at 
±2.5  a . 


Figures  1  through  3  show  height  profiles  of  the  C2  measurements  from  the 
radar,  thermosonde,  and  scintillometer  together  with  two  model  estimate 
profiles.  All  the  radar  data  used  in  this  preliminary  study  were  taken  with 
the  antenna  pointed  20°  towards  the  east.  The  thermosonde  and  model  profiles 
have  been  smoothed  with  the  same  Gaussian  filter.  One  model  profile  gives  the 
total  C2  including  the  humidity  terms;  this  model,  called  model  (radar),  is 
compared  to  the  radar  data.  The  other  model,  called  model  (dry),  omits  the 
humidity  terms;  it  is  compared  to  the  scintillometer  and  thermosonde  profiles. 
Note  that  the  dry  and  radar  model  profiles  merge  together  at  about  ten  km  and 
are  identical  at  higher  altitudes. 

Figure  1  shows  the  first  example  in  this  data  set  that  has  data  from  all 
three  instruments.  Recall  that  the  thermosonde,  scintillometer,  and  dry  model 
profiles  form  one  set  of  profiles,  and  that  these  three  profiles  are  to  be 
compared  with  each  other;  whereas  the  radar  and  radar  model  are  a  separate  set 

of  profiles.  The  agreement  among  the  profile  in  each  set  is  very  good  to 

excellent.  The  radar  data  are  on  the  model  (radar)  curve  except  for  the  lowes 
point  at  A. 85  km;  at  that  height  the  model  (radar)  is  smaller  than  radar 

measurement.  The  next  example,  shown  in  Figure  2,  is  a  daytime  example,  so 

there  are  r.o  scintillometer  data.  Figure  3  shows  the  next  nighttime  example. 
In  both  cases  the  model  (radar)  is  smaller  than  the  radar  data  for  the  lowest 
two  to  three  range  gates.  The  data  from  the  other  range  gates  fit  the  model 
(radar)  well. 


LOG  CN2 

Figure  1.  Height  profile  of  Log  C2  for  02-03  hours, 

8  June  1988  'll.  The  large  solid  dot  is  the  median  of 
the  radar  data  for  the  hour,  and  the  horizontal  bars 
give  the  extreme  values  of  the  radar  measurements  for 
the  hour.  The  solid  line  is  the  thermosonde  profile; 
the  balloon  was  launched  at  02:04  UT .  The  hatched  area 
gives  the  range  of  scintillometer  measurements  ever  the 
hour.  The  long -short  dashed  line  is  the  model  (radar) 
profile,  and  t’^o  dashed  line  is  model  (dry),  which  is  to 
be  compared  w i *'h  the  thermosonde  and  scintillometer  profiles. 
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Figure  3.  Same  as  for  Figure  1  except  for  00-01  hours, 

11  June  1988  UT ;  the  balloon  was  launched  at  00:44  UT. 


Figure  2  is  not  a  typical  example  of  the  thermosonde  measurements  made 
during  the  day.  For  flights  later  in  the  afternoon,  the  thermosonde  data  were 
much  larger  than  the  nighttime  measurements.  This  day/night  effect  is  not  yet 
understood. 

All  observational  methods  used  in  this  study  make  several  fundamental 
assumptions  to  derive  a  C2  value  from  the  raw  data.  The  most  important  assump¬ 
tion  is  that  the  mixing  is  due  to  turbulent  flows,  and,  further,  that  the 
turbulent  flows  are  homogeneous  and  isotropic  and  that  the  observing  scales  are 
in  the  inertial  subrange.  Since  the  measurements  and  model  estimates  used  in 
this  preliminary  study  are  generally  consistent,  these  simplifying  assumptions 
and  analysis  in  terns  of  the  C2  structure  parameter  are  generally  useful.  The 
relative  importance  of  convective  mixing,  viscous  damping  at  scales  smaller 
than  the  innerscale,  and  anisotropic  turbulent  fluctuations  will  require  fur¬ 
ther  research. 

SUMMARY  AND  CONCLUSIONS 

An  experimental  campaign  was  conducted  at  the  Flatland  clear-air  VHF  radar 
site  located  near  Champaign-Urbana ,  Illinois,  in  June  1988  to  measure  height 
profiles  of  the  refractivity  turbulence  structure  parameter  C2  and  related 
turbulent  parameters.  This  Flatland  site  was  chosen  because  it  is  located  in 
very  flat  terrain  far  removed  from  mountains,  so  that  orographic  effects  are 
minimized.  Three  different  techniques  were  used  to  measure  the  height  profiles 
of  C2 .  The  50  MHz  clear-air  Doppler  Flatland  radar  and  a  stellar  scintil¬ 
lometer  measured  the  profile  remotely,  and  high  resolution  in  situ  measurements 
of  C2  were  obtained  from  over  20  thermosonde  balloon  flight.  The  balloon 
instruments  also  measured  the  standard  thermodynamic  and  wind  data  with  excel¬ 
lent  height  resolution. 

Model  estimates  were  calculated  from  the  standard  balloon  data  and  com¬ 
pared  with  the  measurements.  Because  the  radar  measurements  are  sensitive  to 
humidity  and  its  gradient,  whereas  the  thermosonde  and  scintillometer  are  not, 
two  model  profiles  were  calculated.  One  model  profile  included  the  humidity 
terms,  called  model  (radar),  and  the  other,  called  model  (dry),  did  not. 

During  the  nighttime,  all  the  measurements  and  model  profiles  are 
generally  consistent.  There  are  two  exceptions  that  occur  systematically 
through  the  data  set.  One  difference  is  that  the  model  (radar)  values  of  C2 
near  five  kilometers  are  consistently  lower  than  the  values  measured  by  the 
radar.  Another  difference  is  that  the  scintillometer  measurements  at  about  14 
km  are  always  smaller  than  the  thermosonde  and  model  estimates,  and  are  near 
the  instrumental  noise  level. 

In  all  the  nighttime  cases,  identical  values  of  all  the  model  parameters 
and  constants  were  used  in  the  model  calculations.  Furthermore,  these  values 
were  identical  to  those  used  in  previous  studies  using  the  Sunset  and  Stapleton 
radars  except  for  one  very  important  parameter.  This  parameter  is  the  constant 
in  the  distribution  function  of  wind  shears.  Its  value  at  Sunset  and  Stapleton 
was  50%  greater  than  its  value  at  Flatland.  This  suggests  that  the  wind  shears 
at  scales  of  the  fine  structure  are  smaller  over  Flatland  than  over  mountainous 
terrain . 
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1.  INTRODUCTION 

The  capability  of  ST  radars  to  provide 
direct  measurements  of  the  momentum  flux 
presents  an  opportunity  for  major  advances 
in  atmospheric  science.  Some  of  the 
issues,  needs,  and  techniques  surrounding 
the  measurement  of  momentum  flux  are 
discussed  by  Fritts  (1984)  and  Schoeberl 
(1984),  among  others.  Tne  vertical 
momentum  flux  was  first  measured  using  the 
two-beam  technique  Dy  Vincent  and  Reid 
(1983),  who  used  mesospneric  observations. 
They  al so  presented  the  body  force  due  to 
the  vertical  divergence  of  the  vertical 
momentum  flux.  Cornisn  and  Larsen  (19d4) 
examined  the  vertical  momentum  flux  at  14.5 
km  over  AreciDo,  but  did  not  look  into  the 
momentum  flux  divergence.  It  appears  that 
no  other  results  from  the  troposphere  and 
stratosphere  are  available. 

The  purpose  of  this  paper  is  to 
present  preliminary  results  of  momentum 
flux  and  momentum  flux  divergence 
calculations  made  using  data  from  the 
Sunset  radar.  In  an  attempt  to  illustrate 
changing  conditions,  we  present  results 
from  a  day  wnen  the  background  wind  speeds 
aloft  changed  abruptly.  The  results  given 
here  are  examples  taken  from  a  more 
complete  study  of  momentum  flux  values 
under  a  variety  of  synoptic  weather 
conditions. 

3.  DATA  AND  METHOD 

The  Sunset  Radar  (Green,  et  al.,  1985) 
is  located  in  a  narrow  mountain  canyon 
15  km  west  of  Boulder,  Colorado,  and  is 
just  east  of  the  continental  divide.  The 
array  antenna  of  this  VHF  (ST)  pulsed 
Doppler  radar  can  be  steered  in  the 
east-west  or  north-south  vertical  plane. 
During  the  experiment  reported  here,  five 
antenna  beam  positions  were  used:  vertical 
and  15°  to  the  east,  west,  north,  and 
south.  Three  consecutive  observations  were 


made  at  each  beam  position  at  90-second 
intervals  and  then  the  beam  was  moved.  A 
full  cycle  could  thus  be  made  each  20 
minutes.  The  relative  locations  of  the 
radar  volumes  are  illustrated  in  Figure  1. 

The  synoptic  situation  early  in  the 
day  chosen,  January  28,  1985,  was 
characterized  by  light  winds  throughout  the 
height  region  sampled  by  the  radar  (from 
about  4  to  14  km).  After  about  15  UT  a 
weak  jet  stream  moved  over  Sunset  with 
winds  at  10  km  over  the  radar  increasing 
from  10  ms*1  to  over  30  ms'1  by  lb  UT . 

Winds  at  all  levels  above  4  km  were 
primarily  from  the  west;  surface  wind  data 
were  not  available. 


Anlenno 


Figure  1.  Schematic  depiction  of  the  radar 
beams  positions  at  Sunset. 
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The  upper  air  soundings  at  Denver, 
about  50  km  southeast,  show  the  tropopause 
rose  from  about  9  km  to  10  km  through  the 
day.  A  strong  Inversion  whicn  was  present 
in  the  lower  troposphere,  near  650  mb, 
early  In  the  period  dissipated  after 
1200  UT.  Further,  the  temperature  profile 
in  the  lower  stratosphere  late  in  the 
period  was  characterized  by  wave-like 
perturbations. 

The  method  of  Vincent  and  Reid  (1983) 
was  useo  for  tnis  calculation.  The 

vertical  momentum  flux,  u'w1  or  v‘w' ,  is 
found  from 

u'w'  (z)  =  [ V2 ( at , R )  -  V2 ( -a ,R ) 3  /  (2  sin  2a) 

where  V  is  the  radial  velocity,  a  is  tne 
antenna  beam  angle,  R  Is  the  range  and  z  is 
the  altitude.  Altitudes  from  4  to  14  km 
were  sampled  at  1  km  intervals.  Decreasing 
signal -to-noise  sometimes  rendered  the 
upper  level  useless.  The  calculation  was 
made  at  a  given  height  only  If  all  beam 
positions  were  sampled  every  cycle  during  a 
full  hour.  Averaging  periods  of  one  hour 
and  three  hours  were  used;  the  results  are 
similar  except  that  the  three-hourly  values 
are  obviously  smootned,  thus,  we  have 
chosen  to  present  the  one-hour  results  in 
order  to  capture  as  much  detail  of  the 
effects  of  transients  as  possible. 

3.0  RESULTS 

Figure  2  shows  the  hourly  momentum 
flux  from  00-18  UT.  Values  for  the 

meridional  component  (v'w’l  are 
typically  a  few  tens  of  cm2/s2  from 
about  00-09  UT,  and  then  increase  in 
magnitude  to  on  the  order  of  one  mzs*2 
after  09  UT.  The  most  notable  event  is 

found  below  5  km  at  IQ-12  UT,  when  v'w' 
reached  over  -20  nrs'S  The  values  for 

the  zonal  component  (u  w  )  are  usually 
less  than  a  few  n^s-2  before  09  UT , 
although  a  few  magnitudes  exceed  10 
m2s*z.  These  results  can  be  compared 
with  the  aircraft  measurements  given  oy 
Lilly  and  Kennedy  (1973).  They  report 


horizontally  averaged  momentum  flux  values 
of  about  8  dynes/cm2  (which  corresponds 
with  2  m2/s2  for  mean  density  of  5  x 
10-9  gcnr3);  although  their  traces  of 
integrated  momentum  flux  snow  large 
variability,  indicating  local  values  range 
far  from  the  mean.  At  10-12  UT  a  maximum 

is  found  in  u'w'  at  the  same  location 

l  i 

where  a  minimum  was  found  in  v  w  . 

After  12  UT  the  values  of  u'w' 
generally  exceed  10  m2/s2  above  8  km, 
with  a  local  minimum  found  at  11  km  at 
18  UT. 

Figure  3  shows  the  vertical  flux 
divergence  of  meridional  and  zonal 
momentum.  These  results  were  computed  from 
the  data  in  Figure  2  by  taking  differences 
across  layers  1  km  apart.  No  smoothing  has 
been  applied  in  an  effort  to  preserve  as 
much  detail  in  the  results  as  possible. 

The  units  used,  10'3  ms-2,  correspond 
to  3.6  m/s/hour  (e.g.,  the  contour  labelled 
10X10-3  ms*2  is  the  same  as  36 
m/s/hour).  In  Figure  3,  large  values  of 
zonal  and  meridional  momentum 
flux  divergence  occur  near  6  km  at 
10-12  UT.  Otherwise,  the  meridional  values 
are  nearly  all  less  than  about  3  units. 
Alternating  periods  of  large  zonal  values 
are  found  above  10  km  at  14-13  UT.  After 
about  11  UT  the  contours  of  zero  zonal 
momentum  flux  divergence  slope  upward  to 
the  right,  while  at  the  same  time  the 
contours  of  zero  meridional  f'ux  slope 
downward  to  the  right.  The  significance  of 
this  pattern,  if  any,  is  not  yet  clear. 

4.  CONCLUSION 

We  have  presented  preliminary  results 
of  the  momentum  flux  and  flux  divergence 
during  a  transient  episode,  as  a  jet  stream 
moved  over  the  radar.  The  zonal  and 
meridional  momentum  flux  and  flux 
divergences  displayed  remarkable  continuity 
with  altitude  in  time,  increasing  in 
intensity  as  lee  waves  and  other  gravity 
wave  activity  developed  while  the  jet 
stream  approached.  The  momentum  flux 
values  observed  compare  favorably  with 
aircraft  measurements  made  over  similar 
topography,  at  least  during  the  early  part 


Fiqure  2.  Time-height  section  of  the  vertical  momentum  flux  at  Sunset  on 
January  28,  1985.  Units:  m2/s2.  (a)  meridional  (b)  zonal. 
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Figure  3.  Time-height  section  of  the  vertical  divergence  of  the 
vertical  momentum  flux  at  Sunset  on  Januar>  28,  1985. 
Units:  10'3  m/s^.  (a)  meridional  (b)  zonal. 


of  the  day.  The  accelerations  due  to  the 
momentum  flux  divergence  seem  rather  large 
at  first  glance,  especially  for  tne  late 
part  of  tne  day.  However,  we  note  that 
there  may  be  compensating  forces  due  to 
effects  not  considered  here,  such  as 
transverse  circulations  or,  likely  more 
important,  scales  of  motion  too  small  to  be 
resolved  by  these  data. 
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1.  INTRODUCTION 

In  recent  years  much  progress  has  been 
made  in  determining  the  spectrum  of  mesoscale 
atmospheric  motions.  The  frequency  spectra  of 
vertical  and  horizontal  velocities  have  been 
determined  in  the  free  atmosphere  by  means  of  the 
nearly  continuous  measurement  of  radial  velocity 
by  wind-profiling  Doppler  radar  (Balsley  and  Car¬ 
ter,  1982;  Gage  and  Nastrom,  1985;  Eckiund  et  al . , 
1986).  In  addition,  wind  measurements  by  com¬ 
mercial  aircraft  collected  during  the  NASA  Clobal 
Atmospheric  Sampling  Program  (GASP)  have  been 
analyzed  co  yield  wavenumber  spectra  in  the  upper 
troposphere  and  lower  stratosphere  that  cover 
scales  ranging  from  a  few  km  to  10,000  km  (Nastrom 
and  Cage,  1985). 

The  synthesis  of  the  observed  acmospheric 
spectra  into  a  dynamical  framework  has  proceeded 
along  several  lines.  Two  conceptual  models  have 
been  proposed  to  explain  the  nature  of  the  spectra. 
The  first  is  analogous  to  the  Carrett-Mur.k  spectrum 
of  internal  waves  in  the  ocean.  The  second  relies 
bn  a  combination  of  Che  internal  wave  spectrum 
and  a  spectrum  of  quasi-horizontal  motions  presum¬ 
ably  associated  with  the  vorcicicy  bearing  mode  of 
fluid  motions.  This  vortical  mode  possesses  verti¬ 
cal  variation  and  is  often  referred  to  as  quasi- 
two-dimensional  or  stratified  turbulence.  In  the 
following  sections  we  shall  briefly  review  the 
current  status  of  these  two  competing  models  of 
mesoscale  atmospheric  spectra  and  present  some  new 
results  of  an  analysis  of  the  dependence  of  fre-  1 
quency  spectra  of  horizontal  velocity  upon  back-  ! 
ground  wind  speed  that  may  provide  a  basis  for 
evaluating  the  relative  contributions  of  waves  and 
turbulence  to  the  spectrum  of  horizontal  motions 
in  the  atmosphere.  j 

!  2.  BACKGROUND 


magnitude  and  shape  of  the  spectrum  of  atmospheric 
|Vertical  motions,  the  VanZandt  model  did  not  take 
them  into  account.  Subsequently,  as  summarized 
In  Eckiund  et  al.  (1986),  it  now  appears  that  a 
jTather  flat,  nearly  universal  spectrum  of  vertical 
velocity  exists  in  the  atmosphere  at  least  under  i 
(light  wind  conditions.  The  possibility  that  this  1 
(vertical  velocity  spectrum  is  nearly  universal  has 
been  reinforced  by  recent  results  of  the  Flatland 
radar  (Green  et  al.,  1988)  tliat  show  that  in  the  ab¬ 
sence  of  significant  orography  flat  vertical  velocity 
frequency  spectra  are  obtained  under  a  wide  range 
jof  atmospheric  conditions.  However,  the  magnitude 
of  this  observed  vertical  velocity  spectrum  is  not 
(Consistent  with  the  VanZandt  (1982)  model  spectrum 
as  was  pointed  out  by  Cage  and  Nastrom  (1985). 


The  original  VanZandt  model  spectrum  was 
Formalized  to  fit  the  observed  atmospheric  horizon¬ 
tal  velocity  spectrum.  If  the  observed  vertical 
Velocity  spectrum  is  assumed  to  be  an  internal  wave 
spectrum,  it  is  possible  co  calculate  the  horizon¬ 
tal  velocity  spectrum  of  internal  waves  by  employ¬ 
ing  the  polarization  relations  and  dispersion  j 
relation  for  internal  waves.  The  spectrum  of  hori¬ 
zontal  motions  due  to  internal  waves  using  this  I 
approach  is  shown  in  Fig.  1.  Also  shown  in  Fig.  l! 
are  observed  horizontal  velocity  spectra.  In  these 
examples,  the  observed  horizontal  velocity  spectra1 
contain  more  energy  than  the  model  spectrum  at  all 
frequencies.  While  this  suggests  tliat  other  proces¬ 
ses  besides  internal  waves  may  be  responsible  for  i 
the  spectrum  of  horizontal  motions,  it  must  be 
kept  in  mind  that  up  to  this  point  no  attempt  has  I 
been  made  to  account  for  Doppler  shifting  effects. ( 
Recently  Scheffler  and  Liu  (1986)  and  Fritts  and 
VanZandt  (1987)  have  been  able  to  model  the  influ-! 
ence  of  mean  wind  speeds  on  internal  wave  spectra  I 
and  their  results  show  that  the  effect  of  Doppler  | 
Shifting  on  atmospheric  internal  waves  is  not  i 
negligible.  ; 


The  idea  thac  the  spectrum  of  mesoscale 
atmospheric  motions  might  be  due  to  an  incoherent 
spectrum  of  internal  wave  motions  analogous  to  the 
Carrett-Munk  spectrum  in  the  ocean  appears  to  have 
originated  with  the  work  of  Dewan  (1979)  and 
VanZandt  (1982).  VanZandt  examined  the  observed 
Atmospheric  spectra  of  horizontal  motions  and 
Ahowed  that  it  was  possible  to  construct  a  univer¬ 
sal  model  of  t lie  atmospheric  spectrum  in  the 
ppirJt  of  Carrett-Munk  tliat  fit  the  atmospheric 
pbservatlons  quite  well.  However,  since  little 
nformation  was  available  at  the  time  on  the 


■informs 


The  possibility  that  the  observed  atmo¬ 
spheric  spectrum  of  horizontal  motions  may  be  j 

dominated  by  processes  other  than  internal  waves  1 
was  considered  by  Lilly  (1983)  and  Cage  and  Nastrom 
(l985a,  1985b).  These  authors  suggested  that  the 
horizontal  velocity  spectrum  may  be  a  manifestation 
qf  the  vortica1  mode  of  fluid  motions  also  known  as 
Stratified  turDulence  or  quasi-two-dimensional  j 
tjurbulence.  Lilly  based  his  analysis  on  earlier  | 
work  by  Riley  et  al .  (1981)  that  showed  that  at  low 
ffroude  numbers  it  was  possible  through  a  scaling 
analysis  to  separate  the  equations  of  motion  into 


Um 


fig.  1.  Comparison  of  observed  frequency  spectra 
of  horizontal  motions  with  the  frequency  spectrum 
of  horizontal  buoyancy  wave  motions  consistent  with 
a  model  tropospheric  vertical  velocity  spectrum. 


a  set  of  equations  governing  wave  motions  and  an¬ 
other  set  of  equations  governing  turbulunce. 
Physically,  at  low  Kroude  number  in  a  stratified 
fluid  it  should  be  possible  to  decompose  the  motion 
field  into  waves  and  stratified  turbulence.  The 
possibility  that  the  observed  atmospheric  velocity 
spectra  contain  contributions  from  waves  and  quasi- 
two-dimensional  stratified  turbulence  has  been 
considered  by  Cage  and  Nastrom  (1985a,  1985b). 

j  Both  points  of  view  discussed  above  accept 
.the  idea  that  the  spectrum  of  vertical  motions  is 
largely  determined  by  a  spectrum  of  internal  waves. 
Only  in  the  importance  attached  to  quasi-two- 
dimensional  turbulence  do  the  points  of  view 
differ.  Doppler  shifting  effects  may  provide  a 
way  to  differentiate  waves  and  turbulence.  The 
calculations  of  Scheffler  and  Liu  (1986)  and  Fritts 
and  VanZandt  (1987)  proviJe  a  means  co  quantify 
the  Doppler  shifting  effect  of  a  mean  velocity  on 
a  spectrum  of  internal  waves.  The  Taylor-transform- 
ation  can  be  used  to  relate  frequency  spectra  of 
turbulence  to  the  wave  number  spectra  of  atmos¬ 
pheric  motions 

5.  'iilK  FIT' LOT  OF  A  MFAN  WIND  ON  INTERNAL 

;  WAVE  AND  TURBULENCE  SPECTRA 

According  to  the  analysis  of  Scheffler 
and  Hu  (1986)  and  Fritts  and  VanZandt  (1987)  the  I 
(influence  of  a  mpan  wind  on  the  spectrum  of  in¬ 
ternal  waves  can  be  parameterized  by  the  quantity  1 


iwhere  in,  is  a  characteristic  vertical  wavenumber, 

i  " 

jU  is  the  mean  wind  speed,  and  N  is  the  Bruni- 
'Vaisala  frequency.  The  influence  of  the  mean  wind 
o n  the  spectrum  of  horizontal  velocity  for  inter¬ 
nal  waves  is  illustrated  In  Fig.  2  whicii  is  taken 
!frora  Fritts  and  VanZandt  (1987).  Figure  2  snows 
'that  qualitatively  the  influence  of  the  mean  wind 
i  is  to  increase  the  spectral  amplitude  at  the  high 
frequency  end  of  the  spectrum  and  to  decrease  the 
spectral  amplitude  at  the  low  frequency  end  of  the 
spectrum.  The  net  effect  of  the  mean  wind  is  Co 
increase  (make  less  negative)  the  spectral  slope 
above  -  2  which  is  the  slope  of  die  model  spectrum 
without  Doppler  shifting.  Note  Chat  this  change 
in  spectral  slope  is  in  the  sense  required  to 
fit  die  observed  spectral  slope  in  Fig.  I.  Note 
also  chat  only  those  waves  with  a  component  of 
their  pnase  velocity  in  the  direction  of  die 
mean  wind  are  Doppler  shifted. 


W/f 


Fig.  2.  Doppler  shifting  effect  of  a  mean  wind 
on  a  model  buoyancy  wave  spectrum  of  horizontal 
velocity.  (After  Fritts  and  VanZandt,  1987.) 

The  Taylor  transformation  can  be  used  | 
to  determine  the  frequency  spectrum  of  turbulence  i 
seen  by  a  fixed  observer  when  a  known  wavenumber  ! 
spectrum  is  advected  past  at  a  given  velocity.  I 
The  frequency  spectrum  is  given  by  ! 


E(f)  -  2n  E(k)/U 

(L) 

f  =  Uk/27t 

(3) 

To  decermine  the  dependence  of  the  wind  speed, 
we  adopt  the  model  wavenumber  spectrum  contained 
in  Fig.  3.  This  model  wavenumber  spectrum  Is  a 
good  approximation  to  the  climatological  CASE 
spectrum  (Nastrom  and  Gage,  1985).  Using  Eqs.  (2) 
and  (3)  we  calculate  the  frequency  spectra  perti¬ 
nent  to  various  advection  velocities  as  illustrated 
in  Fig.  4.  Note  that  unlike  the  wave  spectra  in 
Fig.  2,  the  turbulence  spectra  in  Fig.  4  retain  j 
their  original  spectral  slope.  Their  spectral 
magnitude  Increases  uniformly  at  all  frequencies  | 
yith  increasing  advection  velocity.  Actually,  over 
mountainous  terrain,  the  GASP  spectra  are  modified; 
bomewhat  (Nastrom  ec  al.,  1987)  and  these  modific-i 
jitions  should  be  taken  into  account  in  a  more  | 

complete  analysis. 
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Fig.  3.  Model  horizontal  wavenumber  spectrum  of 
horizontal  velocity  based  on  GASP  spectra. 


Fig.  4.  Doppler  shifted  frequency  spectra  of 
horizontal  velocity  determined  by  Taylor  trans¬ 
formation  of  the  model  wavenumber  spectrum 
contained  in  Fig.  3. 

4.  THE  OBSERVED  DEPENDENCE  OF  THE  FREQUENCY 

SPECTRUM  OF  HORIZONTAL  VELOCITY  UPON 
WIND  SPEED  AT  PLATTEVILLE,  COLORADO 

The  50  MHz  Doppler  radar  located  at 
Platteville,  Colorado,  has  been  observing  horizon¬ 
tal  and  vertical  velocities  routinely  since  1981. 
•For  the  purposes  of  the  present  paper  accumulated 
data  from  1981-1984  were  analyzed  to  determine 
the  dependence  of  the  frequency  spectra  upon 
background  wind  speed.  The  frequency  spectrum  of 
:the  observed  wind0  were  calculated  for  3  hour  and  i 
j9  hour  periods  and  data  were  stratified  according 
jto  mean  horizontal  wind  speed.  Composite  spectra  ■, 


were  formed  for  each  wind  speed  interval.  Compo¬ 
site  3-hour  spectra  for  the  zonal  wind  are  shown 
in  Fig.  5.  These  spectra  show  a  general  increase 
in  spectral  amplitude  with  increasing  wind  speed. 
Spectral  shape  does  not  appear  to  be  greatly 
affected  by  increasing  mean  winds  over  the  range 
2  to  20  ms 
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Fig.  5.  Observed  3-hour  frequency  spectra  of 
horizontal  zonal  velocity  ac  Platteville,  Colorado 
during  summer  months  stratified  by  wind  speed. 

3.  COMPARISON  OF  THE  OBSERVED  DEPENDENCE  OF 

THE  HORIZONTAL  VELOCITY  SPECTRUM  ON  MEAN 
VELOCITY  WITH  THE  WAVE  AND  TURBULENCE 
MODELS 

In  this  section  we  compare  quantitatively 
the  dependence  of  the  observed  frequency  spectrum 
of  horizontal  velocity  upon  mean  wind  speed  with 
the  Doppler  shifting  effect  discussed  in  Section 
3.  This  is  accomplished  by  determining  the  change 
in  spectral  amplitude  with  changing  wind  speed. 
Since  we  are  concerned  primarily  with  the  relative 
magnitudes  of  the  spectra,  all  changes  in  spectral 
amplitude  are  determined  relative  to  the  spectral 
amplitude  of  tne  non-Doppler  shifted  (or  lowest 
wind  speed)  spectrum.  Changes  in  the  spectral 
amplitude  are  determined  as  a  function  of  wind 
speed  at  several  frequencies.  The  frequencies 
chosen  for  the  comparison  are  the  Brunt-Vaisala 
frequency  N,  a  frequency  equal  to  .IN,  and  a 
frequency  equal  to  .01N.  The  observed  frequency 
spectrum  does  not  extend  to  low  enougli  frequencies 
to  compare  with  the  wave  and  turbulence  models  at 
. 01 N  but  the  changes  In  spectral  amplitude  for  the 
models  have  been  Included  here  for  completeness 

:  In  order  to  determine  the  changes  in  j 

spectral  amplitude  for  the  wave  model  it  is 
necessary  to  assign  a  value  to  the  characteristic  | 


wavenumber  tn  .  We  have  tor  this  purpose  adopted  » 
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a  value  of  =  .75  x  1U  cptn  which  is  consistent 
with  values  anticipated  by  Fritts  and  YanZaadc 
(1987).  For  the  troposphere  pert.'nent  to  the 
summertime  Plutte.ville  spectra  contained  in  Fig.  5 
we  adopt  N  =  1.67  x  10  cps.  With  these  choices 
for  and  N,  3  =  .450. 

The  results  of  the  comparisons  are  con¬ 
tained  in  Fig-  6.  Altogether  there  are  six  curves 
plotted  in  Fig.  b.  Each,  curve  gives  the  spectral 
density  ratio  as  a  function  of  wind  speed.  Curves 
with  positive  (negative)  slope  have  increasing 
(decreasing.)  spectral  amplitude  with  increasing 
wind  speed.  Two  curves  are  plotted  for  the 
observed  change  in  spectral  amplitude  correspond¬ 
ing  to  frequencies  of  N  and  .IN,  respectively. 
These  c.uc.es  show  a  modest  increase  in  slope  with 
decreasing  frequency.  The  turbulence  model  result 
does  not  depend  on  frequency  so  that  only  a  single 
curve  is  drawn.  This  curve  has  a  positive  siope 
intermediate  to  the  slopes  of  the  two  curves  for 
the  observed  spectra.  Three  curves  are  plotted 
for  the  wave  model.  These  curves  show  the  depen¬ 
dence  upon  wind  speed  of  .he  model  wave  spectrum 
for  frequencies  of  N,  .IN  and  .01N,  respectively. 
The  slopes  for  me  three  curves  decrease  rather 
markedly  vi.h  decreasing  frequency.  Tills  is  the 
opposite  sense  of  the  more  modest  change  noted 
previously  for  me  observed  spectra. 
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1  Fig.  6.  Comparison  of  the  observed  dependence  of 
|  the  horizontal  velocity  spectrum  on  wind  speed 
i  with  wave  and  turbulence  models. 

:  6.  CONCLUDING  REMARKS 

In  this  paper  we  have  considered  the 
dynamical  processes  that  may  be  responsible  for 
i  the  observed  mesoscale  atmospheric  wind  spectra. 

I  We  have  addressed  the  issue  of  whether  the  observ- 
|  ed  spectrum  of  horizontal  velocity  is  due  primur- 
|  LLy  to  waves  or  turbulence  by  comparing  the 
!  dependence  of  the  observed  horizontal  velocity 
j  spectra  on  wind  speed  with  the  Doppler  shifting 
effect  anticipated  for  a  model  wave  spectrum  and 
for  a  model  turbulence  spectrum.  The  results 
i  show  that  the  observed  spectra  do  not  follow 
j  either  the  turbulence  model  or  the  wave  model 
|  very  closely.  However,  the  turbulence  model  seems 
|  to  fit  the  observations  more  closely  than  does 
I  -'“C  wave  rno-JcJ  . 


I  The  wave  and  turbulence  models  employed 

here  must  be  considered  as  approximate  and  die 
results  shown  are  necessarily  of  a  statistical 
nature.  The  observations  are  limited  to  a  single 
iocation  and  may  not  be  typical  of  other  locations 
For  these  reasons  the  results  presented  here  must 
be  regarded  as  preliminary. 

i 
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1.  1 NIRCOL  CTION 

Clear-air  Doppler  radars,  also  called 
wind  profilers,  have  the  almost  unique  capability 
to  profile  the  vertical  air  motion  over  the  radar 
or.  a  nearly  continuous  basis.  Clark  et  al.  (198o; 
point  out  that  5w  3z,  the  vertical  divergence 
component,  may  easily  be  obtained  from  these  mea¬ 
surements  through  simple  height  differencing.  This 
technique  is  sensitive  to  all  scales  of  motion  and 
should  be  coopuimentary  to  values  of  horizontal 
velocity  divergence  determined  from  the  NWS  rawin- 
sonde  network  cr  those  obtained  from  the  future 
clear-air  Doppler  radar  netuorks  (see  Chadwick  and 
Hassel,  1987>,  using  techniques  similar  to  those  of 
Zamora  et  al.  (1987). 

Vertical  divergence  observations  seem¬ 
ingly  associated  with  the  passage  of  a  front  are 
cescribed  in  Clark  et  al.  (198b).  However,  Chese 
cbser cations  were  made  at  the  Sunset  radar  in  the 
foothills  of  the  Colorado  Rocky  Mountains  and  were 
complicated  by  orographic  effects  (such  as  lee 
waves)  ar.d  the  convective  activity  accompanying 
the  front  . 

Here  we  present  a  case  study  of  the 
vertical  motion  and  vertical  divergence  observed 
during  the  passage  of  a  cold  front  under  more 
favorable  conditions.  The  data  presented  are  from 
the  Fiatland  radar  (Green  et  al.,  1988)  which  is 
located  in  the  extremely  flat  terrain  near 
Champaign,  Illinois,  far  from  any  mountains,  and 
there  was  little  convective  activity  associated 
with  the  frontal  passage. 

2.  CASE  STUDY 

The  NUS  500  mbar  charts  at  1200  UT  April 
29  and  30  show  strong  north-westerly  flow,  increas¬ 
ing  north  of  the  radar,  and  associated  with  a 
center  of  low  pressure  to  the  north-east.  This 
condition  is  often  associated  with  subsidance  and 
clear  skies.  A.  cold  front  patted  over  the  Flat’.and 
radar  at  about  0200  UT  or.  the  30th  of  April.  Thie 
front  advanced  steadily  from  the  North  at  about 
JOkm/hr  (8  m/s).  The  NWS  surface  analysis  charts 
show  the  front  390km  north  of  the  radar  at  1200  UT 
April  29,  over  the  Fiatland  radar  around  0200  UT, 
end  330km  south  of  the  radar  at  1200  UT  April  30, 
1987.  There  was  no  precipitation  aesociated  with 
this  system  as  it  paased  over  the  radar.  The 
weather  observer  at  Willard  airport,  8km  east  of 
the  radar,  reported  clear  skies  except  for  scat¬ 
tered  clouds  at  3ka  (10,000  feet)  from  10-5-1545  UT 


on  the  29th,  and  again  from  2345-02-5  UT  on  the 
30th.  The  observer  was  off  duty  between  0-45  _T 
until  1045  UT  the  next  day,  tut  at  1145  UT 
scattered  clouds  at  3km  were  again  noted,  while 
at  12-5  UT  there  were  clear  skies.  The  basically 
clear  skies  are  consistent  with  the  expectation 
cf  subsidance  associated  with  the  upper  level  n.ortr.- 
westerly  flow. 
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April  29-30,  1987 

gure  1.  Vertical  velocity  time  series  from  the 
■latland  radar  covering  April  29-30,  1987.  The 
horizontal  axis  is  marked  in  hours  UT,  while  the 
vertical  axis  shows  the  height  of  each  time  series 
in  km  MSL.  The  velocity  scale  is  shown  at  the 
lower  right. 

The  vertical  velocity  data  for  the 
48-hour  period  cantered  on  the  frontal  passage  is 
presented  at  full  time  resolution  in  Figure  1. 
Focussing  our  attention  on  the  lower  heights  the 
traces  show  variability  roughly  within  a  range  of 
plus  or  minus  20cm/s.  An  exception  to  this  occurs 
between  the  hours  of  1800  to  2400  UT  on  the  29th. 
During  this  more  active  period,  w  ranges  between 
plus  or  minus  about  50cm/s,  and  close  examination 
with  a  straight  edge  shows  that  the  excursions 
lasting  a  few  minutes  or  more  are  basically  in 
phase  at  sll  heights  below  7  or  6  km.  From  the 
airport  observations  described  above,  we  know  that 
the  sky  was  clear  during  this  period,  so  that  the 
oscillations  are  unlikely  to  be  due  to  convective 
activity.  A  more  likely  explanation  is  That  these 
features  are  shear  generated  gravity  waves  launched 
from  the  jet  stream  to  the  north.  The  large 
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amplitude  upward  excursion  between 
Co-V'  and  0100  "T  on  the  30th,  however, 

„s  seer,  only  at  the  lowest  height,  2.1km  MSI.  It 
may  well  be  associated  with  a  weak,  pre-frontal 
cumulus  and  is  not  inconsistent  with  the  s.attered 
cleans  observes  during  this  time  at  Willard  airport. 

For  the  detection  of  large-scale  features 
in  the  flow,  we  clearly  need  to  remove  local  scale 
features  and  signals  from  propagating  waves.  In 
Figure  2  we  show  the  results  of  simple  running  mean 
filtering.  The  top  panel  shows  the  results  of 
contouring  1  hour  means  of  w.  The  center  panel 
in  Figure  2  shows  3w  Fz,  determined  from  the 
difference  of  the  w  values  at  adjacent  heights. 

To  the  extent  that  the  vertical  divergence  is 
compensated  by  the  horizontal  velocity  convergence, 
we  can  say  that  the  light  solid  lines  represent 
horizontal  velocity  convergence,  while  the  dashed 
represent  horizontal  divergence.  The  bottom  panel 
was  obtained  from  the  - w  0 z  data  in  the  middle 
panel  by  performing  a  3  height  running  mean  in  a 
simple  attempt  to  reduce  the  amplitude  of  smali- 
vertioul  scale  features. 


Figure  2.  The  top  panel  shews  contours  of  hourly 
mean  w,  in  10cm  s  intervals,  where  light  solid 
lines  represent  upward  motion,  dashed  iines 
downward,  and  the  thick  solid  line  no  vertical 
motion.  The  middle  and  lower  panels  shows  aw/oz 

-4  - 1 

in  10  s  contour  intervals,  where  the  light 
solid  lines  represent  vertical  velocity  divergence , 
the  dashed  lines  represent  vertical  convergence , 
and  the  solid  line  depicts  the  transition.  The 
lower  panel  is  the  same  as  the  middle  panel,  except 
chac  the  data  has  been  smoothed  by  a  3-height 
running  mean. 


evident  in  the  middle  panel  of  Figure  2. 

4 .  SUMMARY 

The  vertical  velocity  divergence  term  is 
easily  measurable  with,  clear-air  Doppler  radars. 

The  divergence  pattern  observed  depends  ot,  the 
dominating  signal  of  the  various  atmospheric  events 
occurring  simultaneously  above  the  radar.  For 
example,  in  the  case  stud'  presented  here,  propa¬ 
gating  gravity  waves  seemed  t  clearly  dominate 
the  data  for  a  six  hour  in.erval.  During  other 
r  ere  source  of  the  signal  seems  associated 

with  the  pattern  of  upper  atmospheric  flow. 
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DISCUSSION 


The  hourly  mean  vertical  velocity  con¬ 
tours  in  Figure  2  show  basically  downward  motion, 
consistent  with  the  pattern  of  flow  shown  in  the 
NUS  500  xbar  analysis.  Superimposed  are  periods 
of  upward  motion,  which  are  not  yet  accounted  for 
in  "ur  study.  The  front  itself  is  apparently 
unobserved  in  these  lower  tropcsphei -c 
observations,  indicating  that  it  must  be  a  shallow 
feature,  perhaps  confined  to  the  boundary  layer. 
This  is  consistent  with  the  fact  that  its  south- 
ware  progression  is  not  reflected  In  the  500  mfcar 
wind  pattern.  The  pattern  of  vertical  divergence 
shows  the  tendency  of  the  atmosphere  to  arrange 
itself  in  layers  of  alternating  divergence  and 
convergence.  This  Is  especially  noticeable  after 
0200  L'T  on  che  30th.  The  decreased  vertical  scaie 
of  this  layering  associated  with  the  active  period 
believed  die  to  trapped  gravity  waves  is  also 
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1.  INTRODUCTION 

Since  the  development  of  the  ciear-air 
Doppler  radar  technique  (also  called  the  wind¬ 
profiling  or  MST-radar  technique)  at  Jicamarca, 
Peru  (Woodman  and  Guillen,  1974),  and  Sunset, 
Colorado  (Green  et  al . ,  1975),  it  has  been  applied 
to  a  wide  range  of  meteorological  problems  (see, 
e.g.,  Liu  and  Kate,  1985).  Despite  this  rapid 
progress,  research  on  some  important  problems  has 
been  frustrated  by  the  fact  that  most  clear-air 
Doppler  radars  are  near  mountains.  The  resulting 
orographic  effects  act  as  geophysical  noise  on 
observations  of  other  processes.  These  effects 
are  especially  serious  for  studies  of  the  vertical 
component  of  motion.  For  example,  Ecklund  et  al. 
(1982)  found  that  when  the  wind  flowed  over  the 
mountains,  the  variance  of  the  vertical  velocity 
was  strongly  correlated  with  the  wind  speed. 
Nastrom  et  al.  (1985)  found  that  they  could  ex¬ 
tract  the  small  synoptic-scale  vertical  velocity 
only  when  the  horizontal  wind  was  not  from  the 
direction  of  nearby  mountains.  Following  their 
suggestion,  we  have  constructed  a  new  clear-air 
Doppler  radar,  called  the  Flatland  radar,  in 
very  flat  terrain  near  Champaign-l'rbana ,  Illinois. 
We  find  that  the  vertical  velocity  field  over  very 
flat  terrain  is  Indeed  quite  different  from  that 
near  rough  terrain,  and  we  present  observations 
that  suggest  that  the  vercical  velocity  due  to 
other  processes,  such  as  synoptic-scale  motions 
and  gravity  waves,  can  be  studied  by  claar-air 
Doppler  radars  in  very  flat  terrain. 

2.  EXPERIMENTAL  DESIGN 

The  Flatland  Radar  is  lo.ated  at  40.5°N, 
88.4°u,  212m  above  mean  sea  level  (MSL),  about 
8km  west  of  the  Champa ign-Urbana  Airport.  The 
radar  operates  at  a  frequency  of  49.8MHz  (wave¬ 
length,  b.02m),  with  a  pulse  length  and  the  range 
resolution  of  750o.  The  3dB,  two-way  beamwidth 
is  3.2°.  To  minimize  contamination  of  vertical 
velocity  measurements  by  horizontal  winds,  the 
antenaa  was  carefully  leveled  to  within  0.02° 
from  the  vertical.  The  Doppler  spectra  have  a 

a  velocity  resolution  of  5cms  and  an  unaliased 
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velocity  range  of  ±  3.2  ms  .  The  Flatland  radar 
has  been  measuring  the  vertical  velocity  every 
153  seconds  almost  continuously  since  March  2, 
1987.  The  second  phase  of  the  Flatland  radar, 
with  steerable  oblique  beams  to  measure  both 
components  of  the  horizontal  wind,  will  be 
implemented  in  early  1988.  A  more  detailed 
description  of  the  radar  and  some  preliminary- 
results  are  given  in  Green  et  al.,  1988. 
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Figure  1.  Radar  vertical  velocities  averaged  ovtr 
fifteen-minute  periods:  (a)  Flatland  radar  during 
March  1987;  (b)  Platteville  radar  during  March  1981. 
(Panels  (b)  and  (c)  are  from  Ecklund  et  al.,  1982.) 

3.  RESULTS 

Figure  1  presents  time  series  of  15-minute 
averages  of  the  vertical  velocity  in  each  range 
gate  of  three  radars,  each  located  in  a  different 
kind  of  terrain.  Panel  (a)  is  from  the  Flatland 
radar  and  panels  (b)  and  (c)  are  from  tue  Flatte- 
ville  and  Sunset  'adars  in  Colorado,  80  and  16kr 
east  of  rhe  crest  of  the  Front  Range  (■>■  4000m  MSL1, 
respectively.  In  the  Flatland  time  series  the 
smaller  upper  height  limit  and  the  data  gaps 
around  10km  are  thought  to  be  due  to  a  smaller 
signal-to-noise  ratio. 
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In  the  Colorado  time  series  in  Figures 
i  i  b '  and  iu'  there  is  a  striking  alcernacion  of 
"active"  pence?  with  large  variance  and  "quiet" 
periods  with  relatively  small  variance.  Ecklund 
et  al.  (1962)  showed  that  the  variance  is  highly 
correlated  with  the  strength  of  the  5km  zonal 
wind  flowing  over  the  Front  Range,  and  Chey  con¬ 
cluded  that  the  active  periods  are  mostly  due  to 
mountain  waves.  In  cor.trasc,  the  variance  in  the 
Flatland  time  series  in  Figure  1(a)  is  nearly 
always  small,  comparable  with  that  during  the 
quiet  periods  at  Piatteviile  and  even  smaller 
than  the  quiet  periods  at  Sunset. 

Mountains  have  similar  effects  on  short- 
period  vertical  motions.  Figure  2  shows  frequency 
spectra  during  spring  1967  from  the  5.2km  range 
gate  at  Flatland,  plotted  wit.i  thick  curves,  to¬ 
gether  with  spectra  taken  in  southern  France 
during  ALP EX  (Ecklund  et  al.,  1985),  plotted  with 
thin  curves.  The  spectra  are  stratified  into 
quiet  days,  when  the  lower-tropospheric  winds 
were  light,  less  than  5m  s,  and  active  days, 
when  the  winds  were  greater  than  -  20  ms.  On 
quiet  days  both  the  ALPEX  and  Flatland  spectra 
are  flat  to  periods  just  less  than  the  buoyancy 
frequency  at  about  10  minutes.  The  Flatland 
active-days  spectrum  is  similar,  but  flatter  and 
silently  raised.  But  a:  ALPEX  the  strong  winds 
were  northerly  mistral  winds  that  passed  over 
nearoy  low  mountains,  and  the  ALPEX  active-day 
spectra  are  much  steeper,  with  a  slope  as  negative 
as  -5 '3  (indicated  by  the  thick  straight  line) 
and  with  amplitudes  much  larger  at  ail  frequencies. 


Pct*uU  (inmulCS } 

Figure  2.  Frequency  spectra  of  vertical  velocity 
fluctuations.  The  two  chick  curves  are  from 
Illinois  (Flatland)  and  the  thin  curves,  from 
southern  France  (ALPEX).  The  heavy  straight  line 
labeled  F  J  1-  an  approximation  to  the  ALPEX 
active-days  spectra.  The  Flatland  spectra  are 
from  the  5.2km  range  gate  and  are  the  average  of 
13  and  9  spectra,  respectively.  The  ALPEX  spectra 
are  the  average  of  four  750m  range  gates  centered 
from  3.85  to  6.10km. 


Doppler  shifting  of  an  lntrinisic  gravity  wave 
spectrum  by  the  background  wind  (Fritts  and 
VanZandc,  1987).  This  suggests  that  under  these 
conditions  the  vertical  motions  are  predominantly 
due  to  propagating  gravity  waves,  with  only  small 
contributions  from  other  processes. 

4.  CONCLUSIONS 

These  results  show  tha-  vertical  motions 
near  rough  terrain  are  often  dom.naced  by  crograpr.it 
effects,  at  all  periods  ranging  from  minutes  to 
many  hours.  The  absence  of  such  effects  over 
very  flat  terrain  suggests  that  clear-air  Doppler 
radars  can  be  used  to  study  vertical  velocities  due 
to  other  processes,  including  synoptic-scale 
motions  and  propagating  gravity  waves. 
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Moreover,  the  slight  change  in  shape  of 
the  Flatland  spectra  with  Increasing  wind  speed  Is 
not  inconsistent  with  the  change  predicted  due  to 
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1.  INTRODUCTION 

The  capability  of  clear-air  doppler 
radara  to  aonttor  the  air  velocity  vector 
over  e  station  on  a  continuous  basis 
enables  studies  of  a  variety  of  features 
not  adequately  sampled  by  other  currently 
available  or  past  measurement  system.'  (such 
as  balloons),  In  particular,  studies  of 
atmospheric  gravity  waves  using 
conventional  data  sources  have  been  impeded 
due  to  their  relatively  short  time  scales, 
from  a  few  alnutes  to  a  few  hours,  and 
many  questions  remain  unanswered.  Vertical 
velocity  data  from  radar  sites  located  in 
or  near  mountains  all  show  alternating 
quiet  and  active  periods  of  vertical 
velocity  variance  (e.g.,  Ecklund,  et  al., 
1982);  the  quiet  periods  are  believed  to 
reflect  a  gravity  wave  spectrum,  and  the 
active  periods  are  caused  by  the  effect  of 
standing  lee  waves  launched  by  the 
atmospheric  flow  over  rough  terrain.  The 
patterns  of  velocity  observed  during  lee 
wave  events  are  highly  localized  and  It  is 
not  clear  that  general  conclusions 
regarding  mesoscale  dynamics  can  be  drawn 
from  these  episodes.  In  contrast  to  those 
results,  the  vertical  velocities  measured 
by  the  Flatland  radar  in  Illinois  represent 
a  spectrum  of  freely  propagating  gravity 
waves  subject  to  doppler  shifting  by  the 
background  wind  (VanZandt,  et  al.,  1989). 


Processes  which  have  been  suggested  as 
sources  of  gravity  waves  include  flow  over 
orography,  unstable  wind  shears,  and 
convection.  There  have  been  a  few  case 
studies  which  show  the  variance  of  the 
vertical  velocity  Is  effected  by  frontal 
passages  and  near  thunderstorms  and  Jet 
streams  (e.g.,  Gage  and  Balsley,  1978; 
Roster,  et  al.,  1986),  but  in  general  it 
fas  been  difficult  to  separate  the  effects 
of  these  processes  from  lee  wav*  effects  in 
radar  data  obtained  near  mountainous  areas. 
T.ve  purpose  of  this  paper  Is  to  describe 
■Tr  atmospheric  processes  which  oause 
enhanced  variance  of  the  vertical  velocity 
at  Flatland  where  lee  waves  are  not  a 
problem. 


2.  DATA 

The  Flatland  radar  is  located  in  v*r> 
flat  terrain  near  Champaign,  Illinois,  as 
described  by  Green,  et  al.,  IBM.  In  its 
final  configuration  the  radar  will  provide 
horizontal  and  vertloal  wind  measurements 
from  several  beaea.  The  dete  available  for 
this  study  are  from  the  first  phase  of  the 
radar  ayatee  and  lnotuda  only  the  vertical 
velocity  messursmsnti.  Date  are  reoerded 
at  ties  Interval*  of  approximately  2 
alnutes  et  vertical  levels  separated  by 
0.75  ka  between  1.4  end  IB. 6  km.  An 
example  of  a  typical  tie*  eerie*  of 
vortical  velocity  ee  seen  by  the  Fletlend 
radar  was  prtaantad  in  Green  et  el.,  1968. 
Close  lnepaetion  of  such  redar  data  reveals 
periods  when  the  variance  of  vertical 
velocity  is  very  eeell  end  other  tlees  when 
the  verlano#  le  relatively  large.  W*  hav* 
studied  the  period  from  March,  1987, 
through  May,  1988,  and  find  the  "quiet" 
period*  correspond  to  days  when  no 
significant  weather  eyetema  were  located 
over  east  centre)  Illinois  end  the  region 
was  often  doelnated  by  high  pressure. 
Conversely  the  "active"  periods  can  ell  be 
Identified  with  eoee  eynoptio  or  eeeoscele 
phenomenon.  Although  th#  difference 
between  "quiet"  end  "active"  periods  seen 
at  Flatland  Is  siuch  smeller  then  that  at 
radar  sites  located  in  or  near  mountains 
where  le*  wave  effects  are  important,  their 
s i gnl f lclence  for  atmospheric  processes  may 
be  greater.  It  is  during  the  active 
periods  that  enhanced  gravity  wav*  activity 
leads  to  significant  flux#*  of  momentum  and 
energy,  end  th#  divergence*  of  thee# 
variable*  hav#  direct  bearing  on  th#  large- 
scale  flow.  Also,  1 r revere ib I e  mixing 
process##  associated  with  episode*  of  wave 
transience  lead  to  exchange  of  paesiv* 
tracers,  especially  in  th#  presence  of 
strong  gradients  euoh  as  near  upper  level 
fronts  end  tropopaue#  fold*.  Thus,  given 
th#  relatively  Imrgm  portion  of  the  earth’s 
surface  effected  by  proeeesmo  eooh  as  seen 
at  Fletlmnd  et  any  given  time,  their  net 
influenoe  on  the  atmoephero  may  bm  large. 

All  Instance#  of  high  variance  were 
found  to  be  associated  with  a  particular 
event  such  as  *  frontal  passage  or  a 

thur.dei  storm.  There  were  no  active  periods 


*13 


sr::  cat ed  i r  the  radar  data  not  aaaociatad 

w  ;  t  n  a  significant  weather  ivanti  likewise, 
: "spactlon  of  tha  daily  weather  aaps  and 
satallita  pictures  showed  that  all  waathar 
evanti  lad  to  a  aignatura  in  tha  radar  data 
during  thla  pariod.  In  ordar  to  illuatrata 
thaaa  signatures  wa  will  naxt  praaant 
specific  caaes  whan  tha  varianca  waa 
incraasad,  ralating  char ac t a r i a t 1 ca  of  tha 
data  collected  by  tha  Flatland  radar  to  the 
concurrent  ataoapheric  eonditiona. 

3.  CASE  STUDIES 

S  March  1987.  Dry  cold  front  with  light 
winds  aloft.  At  9/03  UTC  a  cold  front  was 
stretched  acroaa  northern  Illinois  between 
Chicago  and  Urbana.  Behind  tha  front  a 
massive  high  pressure  dose  of  cold  air 
cantered  in  southern  Manitoba  puahad 
rapidly  southward  creating  an  Arctic 
outbreak.  By  9/09  UTC  tha  front  had  passed 
Urbana  and  was  located  over  tha  southern 
tip  of  Illinois.  Surface  observations  from 
tha  Chaapaign  airport  indicated  frontal 
passage  at  approx laata 1 y  9/05  UTC. 

However,  no  precipitation  or  eonvaction  was 
observed  and  only  alddle  and  high  clouds 
were  reported.  At  500  ab  a  low  presaure 
trough  was  situated  over  the  southeastern 
U.S.,  and  winds  through  the  troposphere 
were  generally  weak  and  froa  the  north  over 
most  of  Illinois. 

Inspection  of  tha  Flatland  vertical 
velocities  for  this  day  (Figure  1)  shows 
that  at  about  05  UTC  there  was  a  sudden 
increase  of  variance  at  all  levels  in  the 
troposphere.  It  is  evident  that  this 
disturbance  in  the  vertical  velocity  field 
corresponds  to  the  ties  the  cold  front 
passed  the  radar  site.  An  Interesting 
feature  in  this  case  is  that  near  11  ka. 
Just  below  the  tropopause,  the  signal-to- 
noise  ratio  of  radar  echos  was  small  before 
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Fit.  1.  Flatland  vertical  velocities  (cm 
s'*)  for  9  March  1967.  Ties  resolution  is 
roughly  every  2.5  ainutes. 


fronts,  passage  and  rose  sharply  after  the 
large  negative  spike  in  vertical  velocity 
below  4  ka  near  06  UTC.  Although  the 
relationship  of  this  phenoaanan  to  the 
dynamics  of  the  frontal  systea  is  not  well 
understood,  it  is  clearly  an  upper-level 
signature  of  the  Increased  turbulence  on 
the  3 - e  scale. 


The  amplitude  of  the  oscillation  in 
the  troposphere  appears  to  deap  out  with 
time  after  9/09  UTC.  Spectral  analysis  of 
the  data  between  06  and  09  UTC  (Figure  2> 
shows  the  predoainent  period  wae  B.6 
minutes  in  the  troposphere.  The  Brunt- 
Vaisala  period  between  700  and  900  ab  was 
scout  9.0  ainutes  based  on  radloeonde 
soundings  at  Peoria  and  Salea.  We  note  that 
these  spectra  closely  reseable  those 
expected  for  a  spectrum  of  Internal  gravity 
waves  in  ca I  a  or  very  light  wind 
conditions.  In  fact,  the  slightly  positive 
slope  of  the  spectra  is  consistent  in  the 
gravity  wave  aode 1  with  the  -5/3  power  law 
often  observed  for  the  horizontal  velocity 
spectra. 
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Fig.  2.  Spectral  analysis  of  vertical 
velocities  for  06-09  UTC  9  March  1987. 


The  close  association  of  the 
oscillations  in  Figure  1  with  the  frontal 
passage  is  in  contrast  with  the  results  of 
Brodhua,  et  al.  (1976),  who  reported  that 
wave  amplitudes  at  the  surface,  as  seen  by 
microbarographs,  began  to  increase  4-5 
hours  before  the  passage  of  cold  fronts. 
Gedzeltnan  (1963)  also  reports  that  wave 
amplitudes  are  larger  than  average  when 
extratropical  cyclones  are  approaching.  Ue 
believe  that  the  unusual  factor  in  the 
present  case  is  that  the  winds  aloft  were 
very  I i ght . 

14-15  Apr i I  1967,  Cold  front  with 
pr ec i p 1 ta 1  ton .  At  14/12  UTC  a  low  at  the 
surface  was  located  over  northwestern 
Missouri,  with  a  wara  front  extending 
southeastward  across  central  Illinois  and 
southern  Indiana  and  with  a  cold  front 
curving  south  across  Missouri  (Figure  3a). 
Rainfal  I  was  recorded  throughout  the 
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*  r  oposphe  r  ■.  c  heights  is  reminiscent  of  the 

3  larch  case  discussed  earlier,  end 
spectral  analysis  again  showed  a  dominant 
period  of  occl Nation  Just  longer  than  the 
B runt - Va i sa 1  a  period. 
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29  Hay  1987.  This  case  illustrates  a  less 
draeatic  source  of  increased  variance, 
naeely,  dry  convective  currents  caused  by 
afternoon  heating  of  the  surface.  Surface 
reports  froe  Champaign  indicatad 
temperatures  in  excess  of  30°  C  and 
scattered  clouds  with  bases  at  1.5  km, 
rising  to  S  km  after  22  UTC.  Vertical 
velocities  for  this  day  (Figure  5)  showed 
lr.creaeed  variance  between  IS  and  22  UTC 
(12  to  16  local  standard  time1, 
particularly  below  about  u  km.  Secauae 
there  was  no  organized  activity  occurring 
at  this  time  at  the  su rface  or  in  the  lower 
troposphere  we  conclude  that  the  enhanced 
vertical  sat  ion  seen  here  was  related  to 
surface  heating.  Inspection  of  the  data 
for  all  tuiaer  reveals  this  type  of 
enhanced  variance  is  frequently  seen  in  the 
afternoon  hours  on  warm,  sunny  days. 

The  effects  of  afternoon  convection 
are  generalized  in  Figure  6  which  shows  the 
diurnal  variation  of  the  standard  deviation 
over  all  data  from  June-August  1987.  In 
the  aid  -  and  lower  troposphere  there  is 
largest  activity  from  12-19  LT.  Greatest 
values  are  below  about  3  km,  with  only  a 
slight  maximum  at  higher  levels.  The 
increase  at  low  altitudes  is  likely 
associated  with  convective  plumes,  while 
the  increase  at  upper  levels  say  bs  due  to 
the  convective  waves  launched  at  low  levels 
discussed  by  Kuettner,  et  al.  (1987). 


4.  CONCLUSIONS 

Several  case  studies  of  synoptic 
weather  events  and  Flatland  radar  data  have 
been  presented.  These  were  chosen  to 
illustrate  broad  features  although  it  must 
be  noted  that  the  time  series  of  vertical 
velocity  are  different  in  each  ease  despite 
a  similarity  In  synoptic  situations. 
However,  several  generalizations  can  be 
made:  1)  All  weather  events  resulted  in  a 

marked  increase  in  the  variance  of  vertical 
velocity,  including  thunderstorms,  fronts, 
synoptic  seals  cyclones,  upper  level 
disturbances,  and  convection  currents 
caused  by  surface  hestlng.  (2)  There  do 
not  appear  to  be  any  instances  of  spurious 
or  unexplained  bursts  of  activity  in  the 
radar  data.  All  easas  of  increased 
variance  were  found  to  be  related  to 
synoptic  or  rassoscals  events.  (3) 
Preliminary  results  indicate  that 
tropospheric  gravity  waves  associated  with 
fronts  and  thunderstorms  sometimes  lead  to 
quas 1 -monochromat 1 c  signatures  in  the  radar 
data . 

Results  from  tbs  Flatland  VHF  radar 
are  consistent  with  end  expend  upon  past 
data  from  surface-based  observations.  The 
high  temporal  and  spatial  resolution  of  the 
radar  data  and  the  ability  to  measure  all 
three  wind  components  simultaneously  offer 


Fig.  6.  Diurnal  variation  of  the  standard 
deviation  of  vertical  velocity  for  ail  data 
collected  from  June  -  August  1987.  The 
ordinate  is  offset  by  5  cm  s-1  for  each 
level  in  the  vertical. 
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1.  INTRODUCTION 

The  measurement  of  the  refractivity 
structure  parameter  C2  is  of  fundamental  im¬ 
portance  to  the  understanding  of  the  propagation 
of  radio  and  optical  waves  through  the  atmos¬ 
phere.  It  is,  therefore,  of  considerable 
interest  to  the  design  and  siting  of  many  types 
of  remote  sensors  and  astronomical  telescopes 
and  to  the  study  of  atmospheric  turbulence. 

Until  recently,  most  of  the  clear-air  radars 
were  located  near  mountains,  so  that  previous 
studies  of  C2  by  clear-air  radars  an  1  com¬ 
parisons  between  radar  measurements  and  other 
methods  of  measuring  C2  were  all  made  in  rough 
terrain  (e.g.,  Good  et  al.,  1982;  Green  et  al., 
1984;  Eaton  et  al . ,  1988).  The  resulting 
orographic  effects  may  alter  the  turbulence  in¬ 
tensity  near  the  mountains  as  well  as  complicate 
the  interpretation  of  the  results.  Therefore, 
an  experimental  campaign  was  conducted  in  June 
1988  at  the  Flatland  clear-air  VHF  radar  site 
located  in  very  flat  terrain  far  removed  from 
mountains  near  Champaign-Urbana ,  Illinois,  to 
measure  height  profiles  of  the  refractivity  tur¬ 
bulence  structure  parameter  C2  and  related 
turbulent  parameters. 

Three  different  techniques  were  used  to 
measure  height  profiles  of  C2 .  The  Flatland 
clear-air  radar  (Green  et  al."  1988),  and  a 
stellar  scintillometer  (Ochs  et  al . ,  1977), 
measured  the  profile  remotely,  and  in  situ 
measurements  were  obtained  from  over  20  thermo¬ 
sonde  balloon  flights  (Brown  et  al . ,  1982).  The 
balloon  instrument  also  made  measurements  of  the 
height  profile  of  the  standard  thermodynamic 
parameters  (pressure,  temperature,  and 
humidity),  plus  wind  speed  and  direction.  Model 
estimates  of  C2  were  calculated  from  the  ther¬ 
modynamic  and  wind  measurements  using  the 
numerical  methods  described  by  Warnock  and 
VanZandt  (1985).  In  addition,  optical  measure¬ 
ments  were  made  of  the  transverse  coherence 
length  r0  (Eaton  et  al..  1988),  and  of  the 
isoplanatic  angle  60  (Ecton  et  al.,  1985).  Both 
these  parameters  depend  in  a  weighted  integrated 
value  of  C2  through  the  atmosphere. 

Both  the  scintillometer  and  radar 
measurements  of  C2  are  obtained  remotely  by 
sensing  the  effect  of  the  fluctuations  of 
refractive  index  on  the  propagation  of 
electromagnetic  waves;  the  scintillometer 
measures  the  optical  scintillation  from  a  star, 
“hereas  the  radar  measures  backscatter  from 
these  fluctuations.  At  radio  wavelengths  the 
humidity  and  its  gradient  as  well  as  the  tem¬ 
perature  gradients  contribute  significantly  to 
the  refractive  index;  whereas  at  optical 


wavelengths  only  the  temperature  fluctuations 
are  important.  In  contrast  to  the  remote  sen¬ 
sors,  a  thermosonde  (Brown  et  al . ,  1982) 
measures  the  temperature  fluctuations  directly 
by  measuring  the  RMS  temperature  difference  be¬ 
tween  two  very  fast  temperature  sensors 
separated  by  a  meter.  Comparisons  among  these 
measurements  taken  simultaneously  in  simple 
topography,  provide  an  unique  opportunity  to 
compare  these  different  measurement  techniques. 
In  this  paper  we  present  some  preliminary 
results,  which  emphasize  the  radar  measurements. 

2.  EXPERIMENTAL  SET  UP 

2 . 1  Flatland  Radar 

The  Flatland  radar  (Green  et  al.,  1988) 
is  located  about  8  km  west  of  the  Champaign- 
Urbana  airport  [40.05*N,  88.38’E,  212m  above 
mean  sea  level  (MSL) ] .  This  clear-air  Doppler 
radar  (also  called  wind  profiler,  or  ST  radar) 
operates  at  a  frequency  of  49.8  MHz  (6.02  m 
wavelength)  with  peak  and  average  power  of  about 
10  kw  and  150  watts,  respectively.  All  the  data 
reported  In  this  paper  were  taken  with  a  pulse 
length  of  1.5  km  and  over  sampled  with  a  range 
resolution  of  750  m. 

The  antenna  consists  of  two  collocated 
dipole  arrays  covering  a  57  x  57  m  area,  with  a 
3  dB,  two-way  beamwidth  of  3.2  degrees.  Each  of 
the  arrays  is  composed  of  strings  of  coaxial- 
collinear  dipoles.  In  this  experiment  the 
antenna  beam  was  tilted  twenty  degrees  off  the 
vertical  toward  the  east  or  south  directions. 
With  this  tilt  angle,  echoes  due  to  specular 
scattering  were  essentially  eliminated,  so  that 
unambiguous  measurements  were  obtained. 

The  signal  received  by  the  radar  was 
sampled,  coherently  filtered,  and  Fourier  trans¬ 
formed  to  produce  Doppler  power  spectra  in  real 
time  at  the  radar  site.  The  Doppler  spectra 
have  128  points  with  a  velocity  resolution  of  5 
cra/s  and  an  unaliased  radial  ±3.2  m/s.  Three  of 
these  spectra  were  averaged  together  and  re¬ 
corded  on  magnetic  tape.  Several  steps  were 
required  to  derive  C2  from  the  power  spectra. 
First,  the  region  in  ?he  spectrum  that  contained 
the  signal  had  to  be  identified  (e.g.,  Clark  and 
Carter,  1980).  Next,  the  signal • to-noise  was 
obtained  from  the  area  under  the  signal  of  the 
power  spectrum.  We  calibrated  the  signal-to- 
noise  by  calculating  the  noise  power  from  radio 
astronomy  cosmic  noise  charts  and  the  system 
temperature.  Finally,  we  calculated  the  radar 
reflectivity  from  the  radar  equation  using  the 
measured  system  parameters  and  recorded  trans¬ 
mitted  power  (e.g.,  Green  et  al.,  1979). 
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2 . 2  Stellar  Sc i r t i 1 1 o^e ;e r  Model  IT 

The  stellar  scint i 1 lomerer  Model  II  was 
developed  by  the  Wave  Propagation  Laboratory  of 
NOAA  (Ochs  et  al.,  1977).  To  operate  the  sys¬ 
tem,  a  star  of  second  magnitude  or  brighter  and 
within  u5“  of  the  zenith  is  selected.  The  sys¬ 
tem  uses  a  Schmidt -Cassegrain  telescope 
35.5  cm  in  diameter  to  focus  starlight  from  the 
star  into  an  attached  instrument  package  that 
sequentially  measures  the  scintillation  inten¬ 
sity  at  differ  nt  spatial  wavelengths  ranging 
from  5  to  15  cm.  Since  particular  spatial 
wavelengths  observed  in  the  scintillation  pat¬ 
tern  originate  from  specific  altitude  regions  in 
the  atmosphere,  the  surface -based  measurements 
can  be  combined  with  weighting  functions  to  ob¬ 
tain  C2  height  profiles.  The  system  can 
measure  seven  different  height  regions  of  opti¬ 
cal  turbulence  ranging  from  2.2  to  18.5  km  above 
ground  level  (AGL) .  The  height  weighting  func¬ 
tions  for  these  seven  heights  are  broad,  and  are 
broadest  at  the  highest  altitudes. 

2 . 3  Thermosonde  System 

The  thermosonde  system  consisted  of  two 
parts,  which  were  mounted  below  an  ascending 
balloon  (Brown  et  al . ,  1982).  One  part  was  a 
standard  VIZ  digital  microsonde,  the  other  was  a 
micro-thermo  bridge  linked  to  the  microsonde.  A 
1200  gram  meteorolog-' cal  balloon  was  used  to 
lift  the  instrument  package  with  an  ascent  rate 
of  about  5  m/s.  The  package  was  suspended  from 
90  to  180  meters  below  the  balloon  to  ensure 
that  the  turbulent  wake  from  the  two  meter  bal¬ 
loon  did  not  affect  the  measurements. 

2.3.1  Thermosonde 

The  thermosonde  measured  the  RMS  tem¬ 
perature  fluctuations  between  two  unheated  fine 
wire  tungsten  probes  separated  horizontally  by 
one  meter.  The  minimum  detectable  RMS  tempera¬ 
ture  fluctuation  was  about  0.002°C,  and  the  data 
were  recorded  every  four  seconds,  giving  about  a 
20-meter  height  resolution.  C2  was  calculated 
from  the  observed  temperature  ?luctuations  and 
the  pressure  and  temperature  observed  by  the 
microsonde . 

2.3.2  Standard  upper  air  data: 

The  microsonde  makes  excellent  height- 
resolution  measurements  of  the  standard  thermo¬ 
dynamic  parameters  (pressure,  temperature,  and 
relative  humidity)  and  wind  speed  and  direction. 
The  thermodynamic  data  was  recorded  every  four 
seconds,  which  gave  about  a  20  height  resolution 
The  wind  speed  and  direction  were  determined  by 
using  the  Loran-C  navigator  system,  and  the  data 
were  recorded  every  10  seconds,  which  gave  about 
a  50  i  height  resolution.  These  data  are  useful 
in  describing  the  meteorological  conditions 
during  the  campaign,  and  are  the  input  for  the 

C2  model  calculations  described  below, 
n 

2. A  Model  Estimates  of _ ££ 

Model  calculations  of  C2  were  made 
from  the  standard  upper  air  data;  i.e.,  pres¬ 
sure,  temperature,  numidity,  and  wind  speed  and 
direction.  These  data  are  described  in  section 

2.3.2  above.  The  basic  model  concepts  are  given 
by  VanZandt  et  al .  (1981).  Since  then,  the 
model  has  been  extended  and  revised  considerably 
We  used  the  latest  version  described  by  Warnock 
et  al.  (1985),  and  used  the  numerical  techniques 


described  by  Warnock  and  VanZandt  (1985)  co 
evaluate  the  model  estimates.  To  compare  the 
model  estimates  with  results  from  previous  com¬ 
parisons  with  the  Sunset  and  Stapleton  radars 
(Warnock  et  al . ,  1986;  Warnock  et  al . ,  1988), 
which  are  located  in  and  near  the  mountains, 
respectively,  we  used  identical  values  of  all 
the  model  parameters  and  constants  except  one: 
the  constant  in  the  equation  giving  the  dis¬ 
tribution  of  wind  shears.  This  parameter 
quantifies  the  shear  environment  in  the  range  of 
scales  important  in  the  onset  of  turbulence 
flow;  therefore,  the  evaluation  this 
parameter  allows  the  Flatland  shear  environment 
at  these  scales  to  be  contrasted  with  the  moun¬ 
tainous  environment. 

2 . 5  Description  of  Campaign 

This  experimental  campaign  was  con¬ 
ducted  at  the  Flatland  radar  site  from  6-15  June 
1988.  The  radar  and  r„  measurements  were  made 
as  continuously  as  possible  throughout  both  the 
daytime  and  nighttime.  The  scintillometer  and 
isoplanometer  systems  shared  a  single  telescope; 
the  scintillometer  operated  at  night  and  the 
isoplanometer  operated  curing  the  day  and  in 
twilight.  Thermosonde  balloon  launches  were 
usually  made  after  noontime  ar.d  before  midnight 
local  time.  A  typical  schedule  was  to  launch  a 
few  packages  in  the  afternoon  and  a  few  after 
dark.  The  radar  operated  during  clear  and 
cloudy  conditions,  whereas  the  optical  systems 
operated  during  clear  sky  conditions.  The  opti¬ 
cal  telescopes  were  located  about  forty  meters 
northeast  of  the  center  of  the  radar  antenna, 
and  the  thermosonde  launch  site  was  about  sixty 
meters  east  of  the  antenna  center. 

3.  PRELIMINARY  RESULTS 

As  described  in  section  2.1  above  we 
used  the  cosmic  noise  power  received  by  the 
radar  system  to  calibrate  the  observed  signal- 
to-noise  ratio  to  derive  C2 .  At  the  low  ranges 
there  are  several  practical  problems  with  this 
procedure  such  as  receiver  recovery  after  the 
transmitted  pulse.  Clark  et  al .  (1988)  dis¬ 
cussed  this  problem  for  the  data  taken  with  the 
Flatland  system  operating  during  this  campaign. 
They  found  that  the  cosmic  noise  calibration  may 
be  used  for  the  sixth  range  gate  (4.85  km)  and 
above .  Thus ,  the  radar  measurements  used  in 
this  paper  begin  at  4.85  km  and  continue  until 
the  signal-to-noise  ratio  becomes  too  small  to 
measure  accurately. 

Measurements  made  by  both  the  radar  and 
scintillometer  remote  sensors  are  average  values 
over  both  time  and  space.  Five-minute  data 
averages  were  used  in  this  paper  by  both  sys¬ 
tems,  and  these  profiles  are  relatively  smooth. 
In  contrast,  each  therraosonde  in  situ  measure¬ 
ment  is  an  average  over  four  seconds  which  gives 
about  a  20-meter  height  resolution.  The  four- 
second  thermosonde  data  displays  many  very  thin 
large  peaks  in  C2 .  Frequently,  the  measured 
C2  values  fall  to  the  observing  noise  level  be¬ 
tween  peaks.  Thus,  to  compare  the  thermosonde 
profiles  with  the  others,  we  smoothed  the 
thermosonde  data.  We  used  a  Gaussian  filter 
with  a  -  0.5  km  and  truncated  the  filter  at  ±5 


Figures  1  through  3  show  height 
profiles  of  the  C2  measurements  from  the 
radar,  thermosonde,  and  scintillometer  together 
with  two  model  estimate  profiles.  All  the  radar 
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data  used  in  this  preliminary  study  were  taken 
with  the  antenna  pointed  20°  toward  the  east. 

The  thermosonde  and  model  profiles  have  been 
smoothed  with  the  same  Gaussian  filter.  One 
model  profile  gives  the  total  C2  including  the 
humidity  terms;  this  model,  called  model 
' radar),  is  compared  to  the  radar  data.  The 
other  model,  called  model  (dry),  omits  the 
humidity  terms;  it  is  compared  to  the  scintil¬ 
lometer  and  thermosonde  profiles.  Note  that  the 
dry  and  radar  model  profiles  merge  together  at 
about  ten  km  and  are  identical  at  higher  al¬ 
titudes. 

Figure  1  shows  the  first  example  in 
this  data  set  that  has  data  from  all  chree  in¬ 
struments.  Recall  that  the  thermosonde, 
scintillometer ,  and  dry  model  profiles  form  one 
set  of  profiles,  and  that  these  three  profiles 
are  to  be  compared  with  each  other;  whereas  the 
radar  and  radar  model  are  a  separate  set  of 
profiles.  The  agreement  among  the  profiles  in 
each  set  is  very  good  to  excellent.  The  radar 
data  are  on  the  model  (radar)  curve  except  for 
the  lowest  point  at  4.85  km;  at  that  height  the 
model  (radar)  is  smaller  than  radar  measurement. 
The  next  example,  shown  in  Figure  2,  is  a  day¬ 
time  example,  so  there  are  no  scintillometer 
data.  Figure  3  shows  the  next  nighttime  ex¬ 
ample.  In  both  cases  the  model  (radar)  is 
smaller  than  the  radar  data  for  the  lowest  two 
to  three  range  gates.  The  data  from  the  other 
range  gates  fit  the  model  (radar)  well. 

1988/08/08  02:04  UT 


LOG  CN2 

Fig.  1  Height  profile  of  Log  C2  for  02-03 
hours,  8  June  1988  UT.  Each  large  solid  dot  is 
the  median  of  the  radar  data  for  the  hour,  and 
the  horizontal  bars  give  the  extreme  value  of 
the  radar  measurements  for  the  hour.  The  solid 
line  is  the  thermosonde  profile;  the  balloon  was 
launched  at  0204  UT.  The  hatched  area  gives  the 
range  of  scintillometer  measurements  over  the 
hour.  The  long-short  dashed  line  is  the  model 
(radar)  profile,  and  the  dashed  line  is  model 
(dry),  which  is  compared  with  the  thermosonde 
and  scintillometer  profiles. 

Figure  2  is  not  a  typical  example  of 
the  thermosonde  measurements  made  during  the 
dav.  For  flights  later  in  the  afternoon,  the 
thermosonde  data  were  much  larger  than  the 
nighttime  measurements.  This  day/night  effect 
is  not  yet  understood. 
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Fig.  2  Same  as  for  Fig.  1  except  for  18-19 
hours,  10  June  1988  UT;  the  thermosonde  balloon 
was  launched  at  17:12  UT.  Because  this  was  a 
daytime  flight,  there  are  no  scintillometer 
data . 
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Fig.  3.  Same  as  for  Fig.  1  except  for  00-01 
hours,  11  June  1988,  UT;  the  balloon  was 
launched  at  00:44  UT. 


All  observational  methods  used  in  this 
study  make  several  fundamental  assumptions  to 
derive  a  C2  value  from  the  raw  data.  The  most 
important  assumption  is  that  the  mixing  is  due 
to  turbulent  flow,  and,  further,  that  the  tur¬ 
bulent  flow  is  homogeneous  and  isotropic  and 
that  the  observing  scales  are  in  the  inertial 
subrange.  Since  the  measurements  and  model  es¬ 
timates  used  in  this  preliminary  study  are 
generally  very  consistent,  these  simplifying  as¬ 
sumptions  and  analysis  in  terms  of  the  C2 
structure  parameter  are  useful.  The  relative 
Importance  of  convective  mixing,  viscous  damping 
at  scales  smaller  than  the  innerscale,  and 
anisotropic  turbulent  fluctuations  in  explaining 
the  times  of  disagreement  will  require  further 
research . 
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4.  SUMMARY  AND  CONCLUSIONS 

An  experimental  campaign  was  conducted 
at  the  Flatland  clear-air  VHF  radar  site  located 
near  Champaign-Urbana ,  Illinois,  in  June  1988, 
to  measure  height  profiles  of  the  refractivity 
turbulence  structure  parameter  C2  and  related 
turbulent  parameters.  This  Flatland  site  was 
chosen  because  ic  is  locaced  in  very  flat  ter¬ 
rain  far  removed  from  mountains,  so  that 
orographic  effects  are  minimized.  Three  dif¬ 
ferent  techniques  were  used  to  measure  the 
height  profiles  of  C2  .  The  50  MHz  pulse 
Doppler  Flatland  radar  and  a  stellar  scintil¬ 
lometer  measured  the  profile  remotely,  and  high 
resolution  in  situ  measurements  of  C2  were  ob¬ 
tained  from  over  20  thermosonde  balloon  flights. 
The  balloon  instruments  also  measured  the  stan¬ 
dard  thermodynamic  and  wind  daCa  with  excellent 
height  resolution. 

Model  estimates  were  calculated  from 
the  standard  balloon  data,  and  compared  with  the 
measurements.  Because  the  radar  measurements 
are  sensitive  to  humidity  and  its  gradient, 
whereas  the  therraosonde  and  scintillometer  are 
not,  two  model  profiles  were  calcu'ated.  One 
model  profile  included  the  humidity  terms, 
called  model  (radar),  and  the  other,  called 
model  (dry),  did  not. 

During  the  nighttime,  all  the  measure¬ 
ments  and  model  profiles  are  generally 
consistent.  There  are  two  exceptions  that  occur 
systematically  through  the  data  set.  One  dif¬ 
ference  is  that  the  model  (radar)  values  of  C2 
near  five  kilometers  are  consistently  lower  than 
the  values  measured  by  the  radar.  Another  dif¬ 
ference  is  that  the  scintillometer  measurements 
at  about  14  km  are  always  smaller  than  the 
thermosonde  and  model  estimates,  and  are  near 
the  instrumental  noise  level. 

In  all  the  nighttime  cases,  identical 
values  of  all  the  model  parameters  and  constants 
were  used  in  the  calc ul  i-1 .  Furtheic v. e , 
these  values  were  identical  to  those  used  in 
previous  studies  using  the  Sunset  and  Stapleton 
radars  except  for  one  important  parameter.  This 
parameter  is  the  constant  in  the  distribution 
function  of  wind  shears.  Its  value  at  Sunset 
and  Stapleton  was  50%  greater  than  ics  value  at 
Flatland.  This  suggests  that  the  wind  shears  at 
scales  of  the  fine  structure  are  smaller  over 
Flatland  than  over  mountainous  terrain. 
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1.  INTRODUCTION 

The  clear-air  Doppler  radar  technique 
(also  called  the  wind-profiling  or  MST  radar 
technique)  has  been  applied  to  a  wide  range  of 
meteorological  problems  since  its  development  by 
Woodman  and  Guille'n  ( 1 97 At )  and  Green  et  al . 
(1975).  Nevertheless,  research  on  some  impor¬ 
tant  problems  has  been  frustrated  by  the 
location  near  mountains  of  most  such  radars. 

The  resulting  geophysical  noise  has  been  espe¬ 
cially  serious  for  studies  of  the  vertical 
velocity  w. 

When  the  background  wind  was  light 
Rdttger  (1981)  and  Ecklund  et  al.  (1986)  found 
that  frequency  spectra  of  w  were  flat  at 
frequencies  less  than  the  buoyancy  frequency  and 
broke  to  a  steep  negative  slope  at  higher  fre¬ 
quencies.  Very  similar  w  spectra  are  observed 
in  the  ocean,  which  Garrett  an d  Munk  (1972, 

1975)  modelled  as  a  spectrum  of  Internal  gravity 
waves.  Thus,  we  can  be  confident  that  the 
light-wind  spectra  are  due  to  gravity  waves. 

But  when  the  background  wind  speed  was 
large  at  stations  near  mountains,  Ecklund  et  al. 
found  that  the  spectral  amplitude  increased  and 
the  shape  changed  su  ^  to  become  quite  incon¬ 
sistent  with  the  light-wind  spectra,  as  shown  in 
Figure  1.  They  suggested  that  these  changes 
might  be  due  to  mountain  waves.  Also,  Nastrom 
et  al.  (1985)  found  that  large-scale  vertical 
velor; t ies  due  to  syncptic-scale  motions  could 
be  extracted  from  the  data  when  the  wind  blew 
toward  the  mountains  from  the  plains,  but  not 
when  the  wind  blew  over  the  nearby  Rocky 
Mountains  in  Colorado. 

It  appeared,  therefore,  that  a  clear-air 
Doppler  radar  located  in  flat  terrain  might  be 
used  to  study  vertical  gravity  wave  motions  and 
large-scale  vertical  motions.  For  this  reason 
we  have  constructed  the  Flatland  radar  near 
Champal gn-Ur bana ,  Illinois.  In  this  paper  we 
show  that  the  behavior  of  w  over  flat  terrain  is 
Indeed  quite  different  from  that  near  mountains. 


2.  DATA  AND  METHODS 

The  Flatland  radar  is  located  about  8km 
west  of  C! ampa ign-Urfcana ,  lil'-ois,  at  VO.S'N, 
88.9^,  212m  above  sea  level.  It  operates  at  a 
frequency  of  99.8MHz.  The  antenna  is  a 
60m  «  60m  array  of  coaxial-coil  inear  dipoles, 
with  a  two-way,  half-power  to  half-power  beam- 
width  of  3.2°.  The  data  analyzed  here  were 


taken  using  750m  pulse  lengths  and  range  gates 
centered  from  1.9  to  19.9km,  but  useful  data 
were  usually  obtained  only  from  2.2  to  16.9km, 
with  often  a  region  of  missing  echoes  Just  below 
the  tropcpcuse.  The  Doppler  spectra  have  128 
points  with  a  velocity  resolution  of  5cm/s  and 
an  unaliased  velocity  range  of  ±3-2m/s.  Further 
details,  including  examples  of  the  data,  are 
given  in  Green  et  al .  (1988). 


Fig.  1.  Frequency  spectra  of  w  in  the  tropo¬ 
sphere.  Tne  thin  curves  are  from  southern 
France  (Ecklund  et  al.,  1986)  and  the  two  thick 
curves  ar%  from  Flatland.  The  spectra  labeled 
QUIET  and  ACTIVE  were  obtained  when  the  wind 
speed  was  S  5  and  i  20m/s,  respectively. 

The  vertical  velocity  was  measured  about 
every  2  1/2min  from  2  March  1987  until  April 
1986,  with  brief  interruptions  due  to  power 
failures,  etc.  and  a  break  of  about  six  weeks  in 
the  summer  of  1987.  Power  spectra  were  derived 
from  time  series  of  two  different  lengths,  6h 
and  95h.  If  more  than  three  successive  observa¬ 
tions  were  missing  or  eliminated  by  the  quality 
control  procedures  described  in  Green  et  al . 
(1988),  that  is,  if  there  was  a  gap  greater  than 
lOmln,  then  the  time  series  was  rejected.  The 
data  in  the  accepted  time  series  at  each  al¬ 
titude  were  then  linearly  Interpolated  to 
uniform  153s  intervals.  The  mean  and  a  linear 
trend  were  removed  and  power  spectra  were 
derived  by  Fourier  transform  of  the  residuals. 


717 


Routine  radiosr:,'.  profiles  of  horizontal 
wind  30d  temperatu'e  or  Peoria,  Illinois,  about 
*25-r  no"t  hwest  .  '  ,ie  Flat  land  radar,  were 
obtained  from  tne  National  Climatic  Data  Center. 


3 .  RESULTS:  6h  SPECTRA 

We  computed  6h  spectra  from  March  through 
May  centered  on  the  nominal  times  of  the  routine 
radiosonde  balloons,  0000  and  12002  ( 1 300  and 
0600  90 °W  time).  Because  the  number  of  in¬ 
dividual  6h  spectra  is  large,  we  are  able  to 
study  the  dependence  of  the  spectra  or.  altitude, 
wind  speed,  buoyancy  frequency,  etc. 

In  Figure  2  the  spectra  are  stratified  by 
altitude  from  2.2  to  19.2km.  For  each  mean 
spectrum  the  number  N  of  individual  spectra  that 
enter  the  average  and  the  variance  VAR  in 
(cm-'s)2  are  given  In  the  table  in  the  upper 
righthand  corner.  During  this  period  the 
average  tropopause  altitude  was  about  12km. 

There  is  no  discernible  altitude  variation 
within  the  troposphere  and  the  stratosphere ,  but 
at  long  periods  the  spectral  energy  density  in 
the  troposphere  is  about  a  factor  of  two  larger 
than  In  the  stratosphere . 

These  spectra  closely  resemble  the  QUIET 
spectra  at  ALPEX  in  Figure  1  ,  even  though  the 
mean  horizontal  wind  for  the  present  spectra  was 
about  15m/s,  and  they  are  quite  different  from 
mean  spectra  at  at  sites  near  mountains,  such  as 
Poker  Flat  (Bemra  et  al . ,  1985;  Ecklund  et  al., 

1 936;  and  SOUS*  (Kuo  et  al . ,  1985;  Larsen  et 
al.,  1988). 
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Fig.  2.  Frequency  spectra  stratified  by  al¬ 
titude.  The  mean  spectrum  at  15.6km  Is  nol3y 
because  it  Includes  only  9  spectra. 


We  hav-  also  sorted  the  spectra  at  3 ■  7 
and  5.2km  according  to  the  buoyancy  frequency  N 
over  the  layer  from  700nPi  to  500hPa  (about  3- 
to  5.5km)  at  Peoria.  There  was  not  any  detec¬ 
table  difference  between  the  mean  spectra 
corresponding  to  the  upper  and  lower  quart  lies, 
with  N  -  13-1  and  9  6xl0->  ra d/s,  respectively, 
a  ratio  of  1.9.  This  is  surprising,  since  in 
Fig.  2  the  ratio  of  1.8  between  the  stratosphere 
and  the  troposphere  causes  very  notl cable  dif- 
f  erence3 . 

This  behavior  is  considered  in  more 
detail  in  Figure  3.  where  the  spectra  at  3.7k-. 
are  plotted  as  a  function  of  horizontal  wind 
speed  WS  in  five  bins  with  boundaries  at  0,  3-6, 
5.5,  11,  and  22 m/s.  These  boundaries  were 
selected  so  that  each  bin  contains  at  least  -1'. 
spectra  and  so  that  the  boundaries  lie  near 
minima  in  the  histograms  of  wind  speeds.  The 
spectra  have  been  multiplied  by  WS2  to  separate 
them  vertically  and  divided  by  VAR  to  facilitate 
comparison  with  model  Doppler-shifted  spectra  ;r 
liie  next  section. 
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Fig.  3-  Observed  spectra  and  model  Doppler- 
shifted  spectra  stratified  by  wind  speed.  The 
observed  and  model  spectra  have  been  multiplied 
by  WS2  and  8’,  respectively,  and  divided  by  the 
variance . 

It  Is  evident  that  the  shape  of  the 
spectra  changes  systematically  with  Increasing 
WS.  As  WS  Increases,  the  slope  becomes  more 
negative  at  periods  longer  than  the  buoyancy 
period  1 /N  and  less  negative  at  shorter  periods, 
so  that  the  spectra  become  flatter..-  The  spectra 
at  other  altitudes  behave  in  much  the  same  way. 
It  will  be  shown  In  the  next  section  that  this 
behavior  is  consistent  with  Dopplei — shifted 
spectra  of  gravity  waves. 
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Fig.:--:  :.  shows  ".r  VAR  inc'eas-s  1::  ; 

power  f  WS.  Th:  s  behavior  will  be  discusses 
more  quantitatively  in  t!>?  next  section. 

4.  TJMPARISDN’  WITH  I T  H  L  DOPPLER-SHI  FTE~; 
GRAVITY  WAVE  SPECTRA 

To  '.;  rcr  :  .  swifted 
gr a v i t  v  wave  spectra  have  been  presented  by 
Schemer  ar.d  Liu  ::?$6  and  by  Fritts  and 
Var.Zar.dt  (195").  I r,  this  paper  we  have  chosen 
to  use  the  approach  of  Fritts  and  VanZandt , 
because  we  are  more  familiar  with  It.  Their 
model  has  been  program  numerically  in  order 

ws0'8 


WS  (m/s).  The  lower  set  of  points  are  the 
observed  variances  VAR.  The  upper  set  fitted  by 
the  straight  line  are  the  total  variances  es¬ 
timated  using  the  Doppler-shifting  model. 

to  be  able  to  Investigate  the  effect  of  both 
their  approximations  to  the  physios  of  gravity 
waves  and  to  the  assumed  intrinsic  spectrum. 

We  found  that  the  physical  approximations  had 
little  effect  on  the  res. Its.  On  the  other 
hand,  some  of  the  approximations  to  the  intrin¬ 
sic  spectrum  are  important. 

Fritts  ar.d  VanZandt  used  p  ■  2  in  order 
to  obtain  solutions  in  closed  form.  We  find 
that  in  c-der  for  the  calculated  spectra  to 
agree  witn  the  observed  spectra,  p  must  be 
nearer  p/3,  3S  is  commonly  observed. 

We  also  considered  two  different  extreme 
azimuthal  distributions.  First,  we  used  the 
approximation  of  Fritts  anj  VanZandt,  with  a 
fraction  a  of  the  wave  energy  propagating  in 
the  azimuth  of  the  background  wind  vector  and  a 
fraction  1  -  a  ,  in  the  opposite  direction. 
Alternatively,  we  assumed  that  a  fraction  a  was 
uniformly  distributed  in  the  semicircle  contain¬ 
ing  the  wind  vector  and  a  fraction  1  -  a  ,  in 
the  opposite  semicircle.  We  found  that  if  a  - 
1/2,  then  either  d  i  of.  -  i  but  i  or.  can  fit  the  ob¬ 
served  spectra,  but  that  if  is  not  near  1/2, 
then  the  mod°i  spectra  cannot  fit  the  observed 
spectra . 


The  model  Doppler-shifted  spectra  depend 
on  the  wind  speed  through  the  scaled  background 
wind  speed  B  -  WS/c#,  where  c,  is  the  charac¬ 
teristic  horizontal  phase  speed,  estimated  to  be 
-  5.5m/s  in  the  troposphere  (Fritts  and  Chou, 
1987)  . 

In  Figure  3  the  model  Doppler-shifted 
spectra  are  shown  by  the  curves  labeled  with  B. 
In  order  to  scale  them  In  the  same  way  as  the 
observed  spectra  they  have  been  multiplied  by  82 
and  divided  by  the  variance  over  the  same  fre¬ 
quency  range  as  the  observed  spectra.  It  is 
evident  that  the  shape  of  the  model  spectra 
changes  In  the  same  way  as  the  observed  spectra. 
The  results  at  other  heights  are  similar.  From 
this  we  conclude  that  the  observed  spectra  are 
Doppler-shifted  gravity  wave  spectra. 

It  must  be  kept  in  mind  that  the  model 
Doppler-shl fted  spectra  all  use  the  same  intrin¬ 
sic  spectrum.  In  reality,  the  intrinsic 
spectrum  mu3t  be  a  function  of  wind  speed  or 
actually  a  function  of  the  entire  velocity  and  N 
profile.  The  fact  that  the  Doppler-shifted 
spectra  fit  the  observed  spectra  at  several 
altitudes  and  over  a  wide  range  of  wind  speeds 
shows  that  the  Intrinsic  spectrum  Is  not  a 
strong  function  of  altitude  or  wind  speed. 
However,  the  Intrinsic  spectrum  could  vary 
somewhat  at  large  wind  speeds,  since  for  large  B 
the  model  spectra  are  rather  insensitive  to  the 
shape  of  the  Intrinsic  spectrum. 

VAR  for  the  observed  spectra  was  plotted 
versus  WS  in  Figure  9.  However,  VAR  is  only  a 
fraction  of  the  total  variance  in  the  gravity 
wave  field,  since  part  of  the  wave  energy  has 
been  Doppler-shifted  outside  the  observed  fre¬ 
quency  range,  especially  toward  higher 
frequencies.  The  total  variances  estimated 
using  the  model  are  also  plotted  In  Figure  9. 

For  WS  >  3-5m/s  the  total  variance  varies 
approximately  as  WS'-*,  or,  to  compare  variance 
with  variance,  as  (WS2)0''.  This  increase  must 
be  taken  into  account  when  comparing  spectra 
taken  at  different  wind  speeds. 

5.  RESULTS:  LOW-FREQUENCY  SPECTRUM 

The  mean  of  the  31*  available  individual 
spectra  from  the  ^ time  series  at  5.2km  is 
shown  as  a  usual  log-log  graph  in  Figure  5a  and 
as  an  area-preserving  spectrum  in  Figure  5b, 
where  the  area  under  the  spectrum  is  propor¬ 
tional  to  the  contribution  to  the  total  energy. 
It  is  clear  that  most  of  the  vertical  energy 
lies  near  the  buoyancy  period  and  that  a  sig¬ 
nificant  fraction  of  it  appears  to  lie  at 
periods  shorter  than  the  5.lmin  Nyquist  period. 

For  the  95h  spectrum  the  mean  WS  is 
estimated  to  be  about  lOm/s,  corresponding  to 
B  -  2.  The  model  spectrum  for  6*2,  shown  as 
the  curve  in  the  figures,  fits  the  observed 
spectrum  rather  well  for  periods  shorter  than 
about  6h,  as  expecteu  from  the  previous  dis¬ 
cussion,  but  it  falls  increasingly  below  the 
observed  spectrum  with  increasing  period.  The 
difference  between  the  observed  and  model 
spectra  may  be  largely  due  to  large-scale  verti¬ 
cal  motions.  This  possibility  will  be  explored 
In  future  research. 
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Fig.  5.  Mean  *1511  3pectrum.  a)  Plotted  as  a 
usual  log-log  graph.  The  curve  Is  a  model 
Doppler-3hlf ted  spectrum  with  8  «  2.  b)  Plotted 
as  an  area-preserving  graph.  To  express  the 
ordinate  In  (ra/s) ’/decade ,  multiply  by 
ln10/45  «  3600. 

6.  CONCLUSIONS 

We  f  !  'hat  ti i high-frequency  vertical 
velocity  i*  >  -;  a  are  quite  consistent  with 
Doppler-.  ;  ;  i  gravity-wave  spectra.  We  sug¬ 

gest  that  •  iow-frequency  spectra  are 
domlnateu  by  :„-~o-scale  vertical  motions. 
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